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ON A PROBLEM OF GEOCZE.* 


By Trpor Rano. 


1. Introduction. 
1.1. The Lebesgue area A(S) of a continuous surface § is defined by 


the formula 
A (S) gi. b. lim 


where the greatest lower bound is taken with respect to all sequences of poly- 
hedra II, converging to S in the Fréchet sense, and H(II,) is the area of 


Hi, in the elementary sense.’ 


1.2. Let us denote by A*(S) the quantity obtained by requiring, in 
‘the preceding definition, that the polyhedra II, be inscribed in S. By re- 
stricting the polyhedra, we certainly do not decrease the greatest lower bound. 
Hence A*(S) = A(S). Obviously then, we have surely A*(S) = A(S) if 
A(S) = ~. 


1.3. It may be argued that it would be more natural to consider A*(S), 
defined in terms of inscribed polyhedra, as the area of S. Still, A(S) is 
generally accepted as the area, probably because very little is known about 
A*(S). The problem of determining whether A*(S) = A(S) may properly 
be called the problem of Gedcze, because Geécze was the first to study this 
question in a substantial manner. Several of his many profound papers on 
the area of surfaces are devoted to this subject.2_ The purpose of the present 
paper is to improve the result obtained by Geécze concerning this problem. 


1.4. We shall be concerned with a surface S of the special form 
(1) S: f(x,y), (x,y) elo 


Where J, is the unit square 


and f(v,y) is continuous in For 0S x1, let L(x) denote the length 
of the are x = const., 0: yS1, z=—f(2,y). Similarly, let us denote, for 

* Received March 9, 1942; Presented to the American Mathematical Society at the 
Meeting in Chicago, April 1942. 

For general information concerning the Lebesgue area of a surface, see for 
example Radé [1], Chapter I. 

*See Gebcze [1], [2], [3]. 
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0=y=1, by M(y) the length of the arc y =const., 0 Sa=1, z=f(z,y), 
Geécze proved the equality A*(S) — A(S) under the following assumptions. 


(i) The functions L(x), M(y) are bounded. 


(ii) f(z,y) satisfies a uniform Lipschitz condition with respect to one 
of the variables, say x. That is, there exists a finite constant L, such that 


| (a, y)| S Lo | — | 
0S¢921. 


1.5. We shall extend this result by proving the following theorem. 


THEOREM. Suppose that the function f(x,y), defining the continuous F 
surface (1), is absolutely continuous with respect tox fora.e.yinO0SySlf 
(or, alternatively, that f(x,y) is absolutely continuous with respect to y for 
Then A*(S) =A(S). 


1.6. Let us recall that A*(S) = A(S) if A(S) =o (see 1.2). Inf 
proving the preceding theorem we can therefore assume that A(S) < o. But 
then, by a result of Gedcze,* the functions L(x), M(y) of 1.4. are summable. 
Hence, in proving our theorem, we shall really use the following assumptions. 


(i)* The functions L(x), M(y) are summable. 


(ii)* For ae. y in 0 Sy the function f(7,y) is absolutely con- 
tinuous with respect to 2. 


Stated in this manner, our assumptions may be more readily compared 


with those of Gedcze (see 1.4). 


1.7. Geédcze stated that he found a proof of A*(S) = A(S), for surfaces 
of the form (1), without any restriction upon f(z, y) beyond mere continuity. F 
He also stated that his proof depended upon a thorough study of Riemann-type 
approximations to Lebesgue integrals and also upon a topological study of § 
the sections z= const. of the surface z=f(z,y). While he did publish 
extensive and profound investigations in both of these subjects,® the present 
writer was unable to find in his published work either a general proof of the 
theorem A*(S) = A(S) or sufficient information to construct such a proof. 
The proof presented in this paper for the partial result stated in 1.5 makes 
essential use of the work of Geécze on Riemann-type approximations, but the 
author was unable to utilize the topological results of Geédcze on the sections 


See Gedcze [1]. 5 See Geédcze [2]. 
* See [4]. See Gedcze [2], [3]. 
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z= const. of the surface which he stated were essential for his proof. For 
this reason, it is unlikely that our proof is similar to the one Gedcze had in 
mind for the general theorem. 

1.8. The method followed in this paper, in proving the theorem of 1.5, 
is essentially a refinement of the method used by Geécze* to prove the result 
stated in 1.4. Various changes were necessary on account of the greater 


generality of our assumptions, and substantial simplifications were achieved 


by properly modifying some of the auxiliary quantities introduced by Gedécze 
(see especially 3.7). However, a definite effort was made to disturb as little 
as possible the original line of thought of Gedcze. 

1.9. The literature on the area of surfaces contains many beautiful 
results concerning approximations by inscribed polyhedra. As far as the 
author is aware, no result of a generality comparable to that of our theorem 
in 1.5. has been previously obtained. For this reason, we restrict ourselves 
to a reference to the work of Kempisty [1, 2], because he proved several results 
of great interest. beyond the mere equality A*(S) =A(S). On the other 
hand, his assumptions are far more restrictive than ours (cf. Kempisty [1]). 


2. On Riemann-type approximations. 


2.1. In studying the continuous surface S given by (1), we can and 
shall assume that A(S) < o (cf. 1.6). As a consequence, condition (i)* 
of 1.6 is satisfied. Furthermore, condition (ii)* in 1.6 will be assumed 
throughout the paper. 


2.2. X will designate a subdivision of the interval 0=2=1 of the 
form 


It will be convenient to describe the generic interval of X by writing 
e =x". Henceforth, 2’ will always refer to the left-hand end-point and 
2” to the right-hand end-point of a generic interval of XY. The following 
symbols will be used. || X || will designate the norm of X, that is the maximum 
length of the sub-intervals of Y, q(X) the number of intervals of X, and 
p(X) the parameter of regularity of X, that is the ratio of the longest to the 
shortest interval of ¥. Thus always p(X) 21. We shall also use sub- 
divisions Y of the interval 0 y <1 and the corresponding symbols y’, y”, 
|Y |, q(Y), p(Y). A pair of subdivisions XY, Y determine, in an obvious 
Manner, a subdivision into rectangles of the unit square OS a1, 
0SyZ1. This subdivision will be denoted by XY. 


7 See Geécze [1]. 
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2.3. Let (A) be an interval-function defined for all closed sub-intervals 
AinO0OS=2z=1. Then xy will denote the sum of the values of y(A) for 
the sub-intervals of Y. The symbol Vy will be used in a similar sense. 
If w(r) is a rectangle-function defined for all closed rectangles r: a’ S aS as 
y S=y=y’ in the unit square Jo, then XY, will denote the sum of the values 


of ¥(r) for the rectangles of the subdivision YY. 


2.4. Let us consider now the function L(x) of 1.4. Since L(x) is 
summable in 0 = #1, it is finite a.e. in this interval.. We agree once and 
for all to use only subdivisions Y such that L(2) is finite at every interior 
point of division of X. Of course, L(0) and L(1) may be infinite. For this 
reason, it may happen that in certain summations ‘Vy the terms depending 
upon 2’ = 0 or 2’ —1 are not finite. In such a case, we shall write X*,, 
instead of Xy, to indicate that such exceptional terms are to be discarded, and 
hence their influence should be studied separately. Similar conventions will 
apply to the symbols Y, Y* in relation to the function M(y) of 1.4. 


An important example is furnished by the summations 


X*L(z2’) (2” —2’), X*L (2) (2” — 2’), 
Y*M(y’)(y” —y’), Y*M(y”)(y”—y’). 


It will be a matter of great importance for us to know that by proper choice 


of X and Y these summations may be brought as close as we wish to 


1 J 
f, Morar, f M(y)dy 


respectively, and that furthermore Y, Y may be made to satisfy a variety 
of further conditions necessary in the sequel. The required facts will be 
summarized presently.® 

2.5. Let there be given two constants e, p such that O0<e<1, p> 
Then we have an 7 > 0, 0 < » < 1, such that the following statements hold. 
Given any measurable subset e, of such that |e. | < and any 
measurable subset e, of 0 = y= 1 such that | ey | < », there exists a sequence 
of subdivisions Y,¥’y of the unit square Jo, depending upon e, p, er, éy and the 


functions L(x), M(y), with the following properties. 
(a) 
(b) g(Xv)=—q(Yv), (v—=1,2,---). 
(c) p(Xv)<p, p(¥v)<p, 


| 0. 


* Except for slight changes in wording, these facts are proved in Geécze [3]. Sub 


sequently, these results of Geécze were extended in Hahn [1]. 
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(d) For each v, Xy,; is a refinement of X, and Y>,, is a refinement of Yy. 
(e) L(x) is finite at every interior point of division of Xv, and M(y) 
is finite at every interior point of division of Y>. 


(f) <6 L(2")(a”—2’) 


, 


—y’) — May | <6 
0 


1 
| (y” —y’) — J, M(y)dy| <«. 


2.6. In the sequel. « and p will be arbitrarily assigned and then kept 
fixed. Next, a corresponding », in the sense of 2.5, will be taken and then 
kept fixed. Then the sets e, and ey will be selected in the following manner. 
L(x) being summable, we can select e, to satisfy the following conditions. 

(2) The complement F, of with respect to OS 271, 


is a closed set. (3) L(x) is finite and continuous on F£,.° 


2.7. The selection of e, is more complicated, since we shall bring into 
play condition (ii)* of 1.6. Since M(y) is summable, we can first select a set 
ey, in 0OSy=1 to satisfy the following conditions. (1) | e’y| < 7/3. 


(2) The complement £’, of e’, with respect to 0S yX1 is a closed set. 
(3) M(y) is finite and continuous on #’,. Next, since f(x,y) is absolutely 
continuous with respect to 2 for a.e. y in 0 Sy 1, we can select a subset 

E,’ with the following properties. (1) |e” |< 7/3. (2) The set 
Ny’ = —e,” is closed. (3) For every ye Fy”, f(z, y) is absolutely con- 
tinuous with respect to z. Take now a ye L,”. Since then f(2,y) is abso- 
lutely continuous with respect to 2, for every interval x’ = 2S 2” the length 


sof the are y=y, is given by 


f [1 + 


Consider then any finite system of such ares, with lengths s,,s2,- + +, corre- 
sponding to the given ye H,” and to a system of non-overlapping intervals 


*See Saks [1], p. 72. 
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- - of total length = 1/k, where k is a positive 


integer. Let ox(y) designate the least upper bound of s,; + for all 
such systems, with y and & fixed. By the absolute continuity of the indefinite 
Lebesgue integral, it follows from (2) that for fixed ye HF,” 


wx(y) for k> 


Further, for fixed k, ox(y) is clearly a lower semi-continuous, and hence 
measurable, function of ye H,’”.. Hence we have a subset e,’” of F,’” with the 
following properties. (1) |< 7/3. (2) The set = —e,’” is 
closed. (3) wx(y) is continuous on F,’”, k=1,2,---. Thus we have, on 
the closed set H,’”, a decreasing sequence of continuous functions ox(y) con- 
verging to zero. By the theorem of Dini, the convergence of such a sequence 
is uniform.*® Hence if we define, for y > 0, 


= 1.u. b. wx(y) 


then 0 for 0. 


2.8. We put now ey—e,’+e/’+e,’". Then the complement of e,, 


, 


with respect to 0= y=1, is the set F,’”” which we shall denote by F,, and 


the following statements hold. 

(1) <> 

(2) E,, the complement of e, with respect to 0S y =1, is closed. 

(3) JM (y) is finite and continuous on Fy. 

(4) We have a function w(y), defined for y > 0, such that w(y) 
for y > 0, and such that the following holds. If yo ¢ F,, and if 7’ Sr=2,", 
are non-overlapping intervals in with total 
length S y, then the sum of the lengths of the arcs 


Y=Yoot Srsz,” z=f(2, Yo), 


is o(y). 


2.9. For given 0<¢<1 and p>1, and a corresponding 7» > 0, we 
have now selected the sets e,, e¢,, whose complements, with respect to0 = #1 
and 0 = y=1 respectively, we denoted by Ez, E,. Thereupon we select a 
sequence X,Y, as described in 2.5. Until further notice, ¢, p, e:, ey, Ex, Ey, 
w(y), Xv¥v will be kept fixed. In the sequel, we shall use the Landau symbol 
O in writing relations of the form A = O(B), with the following meaning: 


1? This is a well-known theorem of Dini. Sce any standard text on Real Variables. 


1/kS$ 
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there exists some finite constant T > 0, depending only upon the given con- 
stant p > 1 and the given function z = f(x,y), but independent of e¢, n, ez, ey; 
E,, Ey, o(W), Xv¥v, such that A STB. 


2.10. For the sequence X,Yy we have (see 2.5, (b)) g(Xv) =q(Y>). 
Let us designate this common value by gv. From the conditions (a), (b), (c) 
in 2.5 we infer easily the following relations for the generic intervals 
of Xv, respectively. 


7 =O qv"), 
(qv*) 


2.11. The sequence X,¥v gives rise to two sequences of functions 
l,(z), mv(y) defined as follows: 

ly(x) = — y’)? + (F(a, y”) — 

my(y) = Xy[ —- 2’)? + (f(2”, y) — f(2’, y))714, 06S 951. 
Thus l,(a2), mv(y) represent the lengths of polygons inscribed in the sections 
z= const., y= const. respectively of the surface (1). By condition (d) in 
2.5 we have therefore 
(3) L(t) 


and, by the geometrical meaning of l,(x), 
(4) L(r) for 


Since 1)(2) is obviously a continuous function of x, and since L(x) is con- 
tinuous on the closed set #,. (3) and (4) imply, by the theorem of Dini, that 
(x) converges to L(x) uniformly on £,."' Similarly, it follows that my(y) 
converges to M(y) uniformly on Fy. Therefore, if we define 


= max (L(x) —ly(x)) + max (M(y) —- m,(y)), 


yeEy 
then 
and 
(6) L(x) — = 


M(y) —mv(y) S &, ye Ey. 


2.12. Let d, designate the greatest one of the oscillations of f(z,y) in 


11 See 10° 


a” =O(1/qr), —y =O(1/q), — 2’) —y’) = O(1/o*), 
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the various rectangles of X,Y». From (a) in 2.5 it follows that d,-0 
for v—> 

2.13. By 2.6, 2.8 the functions L(x), M(y) are finite and continuous 
on the closed sets H,, Hy respectively. We can define therefore 


=max L(x), M = max M(y). 


reEz yeEy 
3. Auxiliary rectangle functions. 


-1. We consider the general ter ‘the sequence of s ivisions X)Y 
3.1. We consider the general term of the sequence of subdivisions YY, 
(see 2.5 to 2.11) and a generic rectangle 2’ Sa#@S2v",y Sysy’ of AY). 
The reader is requested to consult 2 whenever in doubt about the meaning 


of the various symbols used. 

3.2. We shall associate with the rectangle 
(7) 
the following four points of the surface (1). 


A: yoy’, 
B: y=y’, 


C: 


D: 


The symbols ABC, ACD will denote the areas, taken with the positive sign, 
of the triangles with the vertices A, B, C and A, (C, D respectively. The 
orthogonal projections of the points A, B, C, D upon the planes «xy, yz, 22 
respectively will be denoted by Aj, Bi, Ci, Di, i=1,2,3. The areas of the 
triangles AjBiCi, AiC;Di, 1 = 1, 2,3, each taken with a convenient sign to 
be explained in a moment, will be denoted by A;B;Cj, AiC;D;. Whatever 


signs we shall select for these areas, we shall have 
(8) ABC = ||, ACD = | rz ||, 
where 1, fr. are vectors, defined by their components as follows. 


(9) Yr, (A A,B,C;), 
r, = (A,C,D,, A,C2D2, AsC3D3). 


The quantity 


Py = X,Y,(ABC + ACD) 


represents then the area, in the elementary sense, of a polyhedron I, inscribed 


in the surface S. By condition (a) in 2.5, the sequence I, converges to the 


> () 


US 


he 
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surface (1) in the Fréchet sense, and hence we have, by the definition of 
A*(S) in 1.2, 
lim Py = A*(S) 2 A(S). 
vw 
3.38. We make the following agreements concerning the signs of 
AiBiCi, AiCyDi. 


A,B,C, = 3(2” —72’) --y’), A,C\D, = (a” —2’) (y’—y’), 
A,B,C, = $(y” — y’) (f(C1)— f(B1)), A202D. = — yy’) (f(D:) —f(A1)), 
A,B,C, = $(2” — 2’) (f(Bi1)—f(A1)), = — 2’) ). 


We use, for brevity, f(A,) instead of f(z’, y’) and so on. 


3.4. We define 


(x — 2’) (y’—y’), 
=3(y" —-y’) (f(D) + £(C1) —F(B1) —f(A1)), 
c= -—2’) (f(Bi) + —f(A1) —f(D)), 
( 


If (4,4,¢) denotes the vector with components 2,b,c, then we have (cf. 
3.2, 3.3) 
(a,b,c) 2%, n= + |.?° 


3.5. The inscribed polyhedron I, (see 3.2) can be represented by an 


equation of the form 


Hy: €Lo. 


We have then, for example, for 7 SxS2”", y=y” the formula 


” j —f B, 


2’). 


— 
Similarly, ¢v(2,y) is a linear function on each of the sides of the rectangle 
(7) of XY,¥y. From this remark we obtain the formulas 


y”’ 


b= (dv y) — dv(a’, y) )dy, 


c= (pv (2, — y’) ) da. 
3.6. We define 


B= — ay, 


J (f(x, — f(a, y’))da, 
t= (a? + + y?)%, 


** || x | designates the length of the vector r: 


Yy 
he 
22 
he ) 
to 
er 
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3.7. In defining B, y we deviate from Geécze ** who works with 


J f(x,y) —f(a’,y)| dy, 


” 


That is, we omit the absolute value signs appearing in the Gedcze ex- 
pressions B, y. The advantages of this modification were noticed in previous 
literature (cf. the last footnote on p. 626 in Radé [3]; see also Morrey [1]). 
The quantities 8B, y will permit us to dispense with a considerable amount of 
most ingenious but quite cumbersome geometrical detail in the method of 
Geédcze. We have to satisfy ourselves, however, that certain important 
properties of B, 7, 7 are shared by B, y, t. 

3.8. We shall need the following facts. 

(a) Lim XYr exists, for every sequence of subdivisions XY of Jo, pro- 
vided only that || X ||} 0, || Y || 0, and the limit is independent of the 
choice of the sequence XY. 

(b) This limit, to be denoted by 7, is = A(S).** 

(c) These statements remain valid for sequences of subdivisions of any 
sub-rectangle yy The corresponding limit will be 


denoted by 

(d) For every sub-rectangle we have 7 St. 

(e) XYt=T for every subdivision XY of Jp. 

(f) As a direct consequence of the preceding statements, we have 
XYr=T for every subdivision XY of Jo, and XY(t—vr)— 0 for every 
X || 0, || || 0. 


sequence of subdivisions YY such that | 
3.9. The proofs of these statements are precisely the same as in the 
case of the original Geécze expressions 8, 7,7. The reader is therefore referred 
to the excellent presentation in Saks [1], Chapter V, especially Theorem 3. 9, 
p. 168, or to Morrey [1]. 
3.10. We define, for the rectangle (7), 


(10) H = f(A,) + f(C1) —f(D). 


Gebcze [1]. 
14In fact, 7 = A(S) for every continuous surface (1). This theorem was sur- 
mised by Geécze who proved it however only in special cases. The general theorem 
was proved by the author (Radé [2]). See also the excellent presentation in Saks [1]. 
These authors used the original Geécze quantities, but the proofs remain valid, with 
obvious modifications, for the modified quantities f, , 7. 


f 
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Further, sx’, $2”, Sy, Sy” Will denote the lengths of the arcs 


ySySy’, 
z—fi(z,y), 
y=, z—f(z,y’”). 


Gz’; Tx", Gy’, oy” Will denote the lengths of the chords of these arcs. We also 
define 


02° = dy, = f "| y)— y)| dy, 
a 


by = | f(x, y’)— ov(a, y’)| dz, Oy" = f | f(a, dv(a, | dz. 
3.11. Put é=4(2’+ 2”). We define then 
p= (pv (x, y’) — f(a, y’)) dx + f (f(x, y’) — y’) )dé + 
x’ 
(f ( dv (2, ydx + J, (dv(z, — f(z, y”) ) dz. 
3.12. From (10), 3.3, (9) we obtain the identity 


= (0,4 — H, —}(2” 


3.13. We subdivide now the rectangle (7) into two rectangles by the 


, 


line c= §=—4(z7 Let us use the superscript 1 to refer to the 
2 ] 


rectangle Sa y SySy’ and the superscript 2 to refer to the 


rectangle €S2S2", ySySy’. We obtain then from 3.6, (10), 3.11 
the identity 

‘ 


3.14. Since A(S) < co by assumption, the partial derivatives df/dz. 
éf/dy exist a.e. in J, and are summable there (see Saks [1], p. 181). By 
the absolute continuity of the indefinite Lebesgue integral we have therefore, 
for y > 0, a function Q(y) such that the following holds. 


(a) 0 for yo 0. 


(b) If m,7,°-*+,7j,* + + is any system of non-overlapping rectangles 
in J, with total area < y, then 


Sf | Of /dx | drdy=Q(y), 
j 
rj 


ZS dedy ay), 
J 
rj 
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4. Inequalities. 

4.1. If a,, a, are vectors, in Euclidean three-space, then it follows from 
the definition of a, + a, that 

u=|ia, |, w—|a+a | 

are the lengths of the sides of a triangle whose area we shall denote by A. 

4.2. From the preceding remark the following inequalities can be 
derived immediately (the first three being, of course, well known). 

+ a2 || | + | ae |, 

(12) Ja a.+ a. + | a—a |, 
(13) Xa! + a | + | — +a, 
In the last inequality, a, X a. designates the vector product of aj, Qo. 

4.3. The proof of (13) is as follows. By a well-known elementary 
formula we have 

A? = p(p— u)(p—v) (p—w), 
where p= (u+v+w)/2. Thus p=u+v. Hence 
APs S 

(14) AS 
Since || a, X ae || = 2A, (13) is proved. 


4.4. Let there be given, in 2, = 272.2, a continuous function f(z). 


We assume that the length s of the arc 


(15) y=f(z) 


is finite. We denote by o the length of the chord of this arc. Let 


(16) 2, y= $(2) 


be the equation of the chord. We assert the inequalities 


(17) f(x) —o(2)| dx 2), 


1 


f f(x) dx S 2s°/*(s — a)”, 
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If we drop the assumption that s is finite we have the inequality 


(19) f(x) — o(x)| dx S — 2,)d, 


“1 


where d is the oscillation of f(z) ina SxS 


Proof. (19) is obvious. (17) and (18) can be derived either geo- 
metrically, by enclosing the area between the are (15) and its chord (16) 
in a parallelogram and in an ellipse respectively, or else by the following more 
analytic argument. Let F denote the maximum of | f(x) —¢(x)| in 


Then 


(20) f(r) —(r)| dx = (4%, —2,)F, 
< 


and we have a €, 2, é < 2, such that 

(21) | —o(é)| =F. 

Let us denote by P, R, S the points (2, f(21)), (v2, f(r2)), (é& 
(é£.¢(€)). From the triangle POR we see that F=—RS cannot exceed 
PR+ RQ. Hence 

(22) F=PR+ROSs. 


(20), (21), (22) prove (17). Next, observe that (2.—2,)F is twice the 
area of the triangle PQYR. Hence, by (14), 


(23) (a2 — <2(PR+ RQ)*2(PR + RO — PO). 
But PR + RQ <s, PQ=c. Hence (18) follows from (20) and (23). 

We shall list now, for easier reference, a number of inequalities for the 
auxiliary quantities defined in 3. The cross-references indicate the sections 
containing the information needed in the proofs. 

4.5. || == O(| H |/qv). See (9) in 8.2, (10) in 8.10, 3.12, 
2.10. 


4.7. |c—y|S6+6, |b—B|S0-+ 6. See 3.4, 3.5, 3.6, 
3.10. 


y |- See 3.4, 3.6, (11) in 4, 2. 


4.8. We have the relations 
(24) 6, = O(S2/qv), 


(25) 6. O[s.* 2 (8, 


and analogous relations for 6,”, 6y, 6". See 8.10, 4.4, 2.10. 
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4.9. The preceding estimates for 6, may be meaningless if x’ 0, 
because then possibly s,,= 0. For 2’ =0 we have the estimate 
6, = O(dv/qv) 
with analogous statements for 62”, Oy, Oy". See 8.10, (19) in 4.4, 2.10, 2.12. 


4.10. Let us return now to 8.138. We shall use the superscripts 1, 2 in 
the sense explained there. We introduce the vectors (see 3.6, 3.13) 


y) — (a, BO, y?), (a(?), y). 


Then (see 3.6, 3.13) 
(26) rao +o |. 


Ry the definition of the vector product we have 


(27) — ay | X HO 
By 3.13 we have 
27) and (28) yield 
(29) Ly? —y™ |S 


By 3.13, 3.6, 3.8, (26) we have 


(29), (13) in 4.2, (30), (31), 3.4, 2.10 yield now 
(32) | —- y) | —r)¥/?], 


4.11. |p|S6,+ 6. See 3.10. 


4.12. | H| = + qv(Oy + Oy") See 8.18, (32) in 
4.10, 4.11. 


4.13. ABC + ACD—n = O[qt?/?(t —1r)¥/? + Oy + 
See (8) in 3.2, 3.4, (12) in 4.2, 4.5, 4. 12. 


4.14. From the definition of H (see 3.10) it is obvious that 
| H| Ssy + sy. 


Hence (see references in 4. 13) 


ABC + ACD—n=O 
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4.15. This last estimate may become meaningless if or = 1, 
since then possibly sy == 0 or Sy” = ©. We derive therefore an alternative 


estimate. By 8.10 and 2.12 it is obvious that 
| H|—O(d,). 


12. 
| ; sv the same cross-references as in 4. 14 we obtain then the estimate 
ABC + ACD—n=O(dv/q). 
4.16. Let there be given now a positive constant K. For each rectangle 
rof X,¥y we define then a quantity Ax as follows. Ax = 1 if for the rectangle 
r all the inequalities 
Sr’, 82", Sy" S K/q, t= K/q’ 
hold. Otherwise A, —0 for that rectangle. We propose to show that 
(33) = O(1/K + 1/q), 
where a denotes the area of the generic rectangle r of Xv¥>. 
4.17. Let us consider first those rectangles of YyYy for which 
(34) > K/qv and 0. 
Let N(x’) designate the number of such rectangles for fixed 2’. Then (see 
1,4, 2.10, 2.4, 2.5) 
N(e’) = = O[ (qv?/K)L(a’) («” —2’)], 
= OL L(2)de + = 


Thus the number of all rectangles satisfying (34) is O(qv?/K) and hence 
their total area is O(1/K’). The sum of the areas of those rectangles of XyYv 
for which «’ 0 is O(1/qv) (cf. 2.10). Hence, the sum of the areas of 
those rectangles of for which > K/qv is O(1/K Applying 
the same reasoning to sz", sy, sy", it follows that the sum of the areas of those 
rectangles of Y,Y, for which the inequalities 


Se’, 82", Sy’, Sy" S K/qv 


do not hold simultaneously, is O(1/K + 1/qv). 
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4.18. Let now N designate the number of those rectangles of Y,Yy 
for which t > K/qv?.. Then (see 3.8) 
NK/q’? =T, 
and hence 
N = 0(q’/K). 
Thus the sum of the areas of these rectangles is O(1/K). Hence (see 4. 16, 
4.17) 
+ 1/q) + O(1/K) = O(1/K + 1/q), 
and (33) is proved. 
4.19. X* v8. Qvs 82", ‘dvs XyY *y8y"/qv == () 


Let us consider the first summation. By 2.4 to 2.10 


= O[LX* — 2’) = O[X*L (2’) (x” — 2’) ] = O(e). 


4.20. We have 
(35) NvY XvY y8y", qv = O[o(pM/K* ) |. 
y’€Ey 
sy’ > K'/?/qp > K'/?/qp 


Proof. It is sufficient to consider the first summation (see 2.6, 2.8, 
2.5, 2.13, 4.16 for notations). Let N(y’) designate the number of arcs sy, 
for fixed y’, that contribute to the first summation. Then by 2.18 


N(y’)K*® _ Mqv 
= I, N(y) 
qv 


For fixed y’, the sum of the intervals 2’ = «= 2” corresponding to these ares 
is therefore (see 2.5, (c) ) 
= N(y’) p/q pM/K*. 


By 2.8, (4) it follows that, for fixed y’, the sum of the lengths of these ares 
Sy is Sw(pM/K*). Since y’ takes on not more than qv values, (35) is 


established. 


4.21. 
== Ole wo(pM 1/K” k**/qv]. 


Proof. For notations see 2.4, 4.16, 2.5, 2.8, 2.18. It is sufficient to 


consider the first summation. We split it into three summations 3, =2, 3s 
as follows. %, is relative to those terms for which y’ee, (see 2.7). 32 i8.. 


relative to those terms for which 
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and sy > K*/q. 
3; is relative to those terms for which 
y eH, and K%/q. 
By 4.19, 3, = @(«). By 4.20, 3, = @(o(pM/K*)). Finally, by 2.10, 4. 16, 
3, S — Ax) K*/qv? = K*0[XVYv(1 — Ag) a] 
= K%0O[1/K + 1/qv] = O[1/K* + K*/q]. 
4.22, 92", Xv¥ Oy, XvY = O(e + dv). 


xv’ y’ €€y €ey 

See (24) in 4.8, 4.19. 4.9. 
4.23. Xv¥ 62", Xv¥vOy, XvY = O(e + dy + 
&'SK/qy 82"SK/qy sy" SK"/?/qv 

Proof. It is sufficient to consider the first summation. We split it into 
two summations where 3, is relative to those terms for which. 2’ ¢ ez, 
and is relative to those terms for which a’ « By 4.22, 3,—=O(e+ dy). 
By (25) in 4.8, the inequality of Schwartz, and (6) in 2.11, 

Rs Es 
82'SK/qp 


== O| (K* 2 qv? 2) qv 2(L(2’) ly 


O| J qv? 2. qv 2 O(K* 25,1 


4.24. We have 


Proof. We start with 
| ABC + ACD = | ABC + ACD—n|+|n—r|. 
By 4.6, 4.7, 4.23 
|n—1 | = Ole + dy + 
By 4.18, 4.23 


(88) | ABC + ACD—n | 
= O[XVY (t — + Ole + dy + 
By 4.16, 3.8, and the inequality of Schwartz, 
= (K*?/qv) (t — 
S (K*?/qv) y(t — 7) = (T — Xv¥ or)”. 
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(37), (38), (39) imply (36). 
4.25. We have 
(40) Xv¥v(1—Ax) ABC + ACD—7z| 
= O{0(1/K + 1/qv) + dv+e+ 0(pM/K*) + 1/K*? + K*?/q)}. 
Proof. By 4.6, 4.7 


By (33) in 4.16 
XvY,(1 O(1/K + 1/qv),; 


and hence, by 3.14, 3.6, 2.1, 
(42) Xv¥v(1—Aw) | SO[0(1/K + 1/q)]. 
From the definition of b (see 3.4) it is clear that 

|b) S3(y”—y) (sv + sy") = O(8y/qv + 


aad (see 2.12) 
|b | —y')dy = 


Hence, by 2.4, 4. 21, 
(43) XvYy(1 b | == O(dy)+ O[XY* (1 Ak) (sy/qv + Sy"/qv) | 


= O[dy + «+ 0(pM/K) + 1/K*¥? + K/2/q]. 
By 4.8, 4.9, 4.21, 


(44) Xp¥v(1—Ax) (By + —= O (dv) + (87/qv + 
= O[ dy + e+ o(pM/K) + 1/K¥? + 
By (41), (42), (43), (44), 


(45) = O{Q[0(1/K + 1/q)] 
+ dy o(pM/K*/*) + 1/K*? + 


By 4. 14, 4.15, 4.21 we obtain similarly 
(46) AXvYy(1 Ax) A BC A CD | 
= O[dy + + 1/K?? + 

(45), (46) imply (40). 

4.26. As a direct consequence of 4.24 and 4.25 we obtain 
| ABC + ACD | 

0{9[0(1/K + 1/qv)] + o(pM/K**) + + 8y + 4 
+ +. — 


/q 
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5. Proof of the theorem. 


5.1. Lemma. For the sequence XvY», described in 2.5 to 2.9 we have 


(47) lim | ABC + ACD —7|STe, 


where T is a constant that depends only upon the given surface (1) and the 
constant p of 2.5. 


Proof. In 4.26 keep K fixed and make vo. By 2.5, 2.12, (5) in 
2,11, (e) and (f) in 3.8 we obtain 


lim Xv¥, | ABC + ACD —r| = O{Q[0(1/K) ] + o(pM/K*¥?)+ 1/K"?}. 


This relation holds for every K >0. Hence, making K — 0, we obtain, 


by 3.14, 2.8, 


lim | ABC + ACD = O(e). 


By the meaning of the symbol O (see 2.9), this last relation proves (47). 
5.2. By 3.2, we have 
XvYv(ABC + ACD) = FE (IL), 
where /(II,) is the elementary area of an inscribed polyhedron 
Hy: 
such that dv(r,y) > f(x,y) uniformly in Jo. Since (cf. 3.8) 
| —T | =| XY¥,(ABC + ACD) — + — t)| 
= X,Y,| ABC + ACD—7|+ XvY,(t—7) 
the relation (47) implies that (cf. 3.8, (f)) 


(48) lim | | 


On the other hand, as noted above, 


(49) max | f(z, y) — @v(a, y)| > 0. 

(2,y) v->00 
We recall that «, p were assigned arbitrarily, and that T in (48) is independent 
of « and of the sequence Xy¥y. Hence, for any positive integer 7, we can 
select first p to satisfy 
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Once p has been so selected, the constant T in (48) is determined. Then we 
select « > 0 to satisfy 

Pe < 1/3}. 
After « has been so selected, we have a sequence Yy}y as described in 2.5 to 
2.9, and (48) and (49) hold. For v large enough, we have then, by (48) 
and (49), 
(50) | —T < 1/), max | f(x,y) — dv(@,y)| < 1/). 


We take such a y, and we designate the corresponding Xy., Yv, ¢v, Uy by 
X@, YY), ¢, 1% to stress their dependence upon the positive integer j. 
We obtain in this manner the following lemma. 

5.3. Lemma. We have a sequence of subdivisions X‘)Y"?, and a 


corresponding sequence of inscribed polyhedra 


WH: (a, y) elo, 


such that the following holds (ct. 2.5 and (50)). 


(a) (2,y4) of(a,y) uniformly in Io. 
j7x 

(b) g(X”) =q(Y’), 

(c) p(X") 1 for] > @. 

(d) 7. 


5.4. The proof of the theorem of 1.5 is now immediate. By the defi- 


nition of A*(S) in 1.2 we have 


(51) lim = A*(S). 

By 3.8, (b) we have 

(52) A(8). 

sy 1.2 
(53) A(S) = A*(8). 


(51), (52), (53) and (d) in 5.3 imply A(S) = A*(S). 


5.5. Infact we have proved more than the mere equality A(S) = A*(S). 
Indeed, we also proved that we have sequences of inscribed polyhedra as 
described in 5.8. The conditions stated there express the fact that the 
approximating inscribed polyhedron I) is such that its projection upon the 


zy plane consists of a number of triangles obtained by subdividing the unit 
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square J, first into rectangles that are practically congruent squares, and then 
drawing one diagonal in each of these rectangles. 
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INTERPOLATED DENUMERANTS AND LAMBERT SERIES.* 
By E. T. BELL. 


1. Notation, summary. The letters denoting integers, r > 1; a, 


are constant and >0; 2,,:--:,2, are variable and =O. The number 
of sets (2,,- - of solutions of 
(1.1) 
is (in Sylvester’s terminology) the denumerant of n with respect to a,,° °°, dy. 
To avoid trivial modifications in the notation, the case a,—=-:-*—a,;=a 
is excluded; D(n|a,---,a) is the coefficient of x" in (1—a2*)-", =1 if 
ai-r r-1 
n=0, =0 if afn, geet Ra II (n+ ja) if a|n. When there can be no 

ambiguity, we abbreviate D(n|a,,- - -,a,) to D(n), and, if necessary, state 
the associated equation (1.1). 

The L.C. M. of a,,- -,a, will be denoted by a. are any r 


different non-negative constant integers. Hence 
(1. 2) D(an; + b;), 0= <a, 
are constants. 
In Sylvester’s method,’ 
D(n) = A(n) +U 
where A(n) is a polynomial in n of degree r — 1, and U, the ‘ undulant’ part, 


contains roots of unity. Sylvester states * that 


r-1 
A(n) = 
j=0 
where A, is the coefficient of x* in the power series expansion of [[ (1 — 
3 1 
Cayley’s method leads to similar conclusions, the periodic terms in D(n) 
being collected in his ‘ circulators.’ Our object is to get rid of both undulants 


* Received January 14, 1942. 

1 For the history of denumerants to 1919, see L. E. Dickson, History of the Theory 
of Numbers, vol. 2 (1920), chap. 3, and for an exposition of Sylvester’s and Cayley’s 
methods, either the papers in vol. 2 of each of these authors’ works or P. Bachmann, 
Niedere Zahlentheorie, zweiter Teil (1910), Kap. 3. The results stated by Sylvester 
in 1857 (loc. cit., No. 17) are repeated, and most are proved, in several later papers. 

2 Loc. cit., 87, last footnote. 
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and circulators, as these are responsible for the most awkward part of the 
calculations, by suitably classifying the arguments n in D(n). 

It will be shown (2) that D(an +b), for any constant integer b such 
that 0 = 6b <a, is a polynomial in n of degree r —1, 


(1. 3) D(an+ b) 
with constant (independent of n) rational number coefficients which may be 
explicitly evaluated in terms of the constants (1.2) with bj =b. 


If in (1.3) we put n=7,,:--,n,, and from the resulting equations 
and (1.3) eliminate we get 
Dian-+-6) 1 
D(any +b) 1 
D(any +6) 1 n, 


The coefficients of the elements in the first column in the expanded form of 
the left are determinants of a familiar type (alternants), whose expansions 
as products of differences are immediate. Or we may apply Lagrange’s inter- 
polation formula to the polynomial in (1.3): 


P(e) (e—nj), (2) = 
j=1 
(1.4) D(an + b) —2 D(anj +b). 


The constants n; are at our disposal. For the choices nj =j—1, or 
nj (1.4) takes a simple form. If nj = j, 
(1.5) D(an +b) = (—1)"(n—1)(n—2)- - (n—r)d¢n(a, db), 
(—1)/D(ja + b) 
n(a,b) => — 


which may be written 


(1.6) (—1)'r! Dian +b) = (n—r) (—1)47 
j=l 

It follows that the r constants cj in (1.3) are linear homogeneous func- 
tions, with constant rational number coefficients, of the r constants D(an; + b) 
for (1.4), and of D(ja+ b) for (1.6), j7=1,---,r. The complete set of 
formulas D(n) for a particular equation (1.1) therefore requires the cal- 
culation of the numerical values of ra constants of the kind (1.2). This may 
be done directly or by the classical methods cited. 


= 
r 
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It will be seen in the proof of (1.3) that b may be any non-negative 
integer. Hence in (1.6) we may take n =O and in the result replace b by 


the arbitrary non-negative integer n, 


(1.7) D(n) + (—1)/(") D(n ja) 


For n = 0 this enables us to calculate D(ra) when D(ja), 7 =1,:--,r—1 
are known, 


(1.8) D(ja) = 0. 


From (1.6) we have 


(1. 9) D(an + b) = O(n"), 
the constant implied in the order function being 
fr 
— ( ) D(ja+b). 
J 


The enumeration of compound partitions being reducible to that of 
simple partitions, there are formulas corresponding to the preceding for com- 
pound denumerants. 

To illustrate (1.5), let r—=3. Then a is the L.C. M. of ay, ds, a3, and 
D(m) refers to 

+ Ao, + A3%, mM: 
2D(an + b) = (n—-2)(n —3)D(at+ 6) —2(n — 1) (n— 3) D(2a + dD) 
(n--1)(n—2)D(3a+ 5), 0=b<a. 


For large values of n, this D(an + b) is approximately 
D(a + b) —2D(2a+ 6) + D(38a+ b)]. 


2. Proof of (1.3). Everything will be proved if (1.3) is. With the 
notation as in 1, let a be a primitive a-th root of 1. Then if a; > 1, there 
is a least integer s; > 0 such that a*) is a primitive aj-th root of 1. In the 
set (¢t; = 1,- -,a; let y,8,- - -,7 be all the unequal 
elements of the set. and let these occur precisely g, d,- - - , ¢ times respectively. 


(If sjt; = siti, mod a, the corresponding elements are equal.) Then, | «| <1, 


> D(n)2" = Il (1 — 


n=0 


= (1—-2)-"(x — y)-9(a 


No 1—2z” has a repeated root, and there are precisely r such binomials. 


Hence none of g,d,- - -,¢ can exceed r. 
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The resolution of the above generating function into partial fractions 
may be written 
— G(x ¥)"* + (—1)’G,(a# y)7 
— D,(« — $)*+ 8)? +---+ 
+ 
in which the capital letters denote constants, and since g Sr, dSr,---,tSr, 
some of the constants G, D,- --,7' may be identically zero. The coefficient 
of y” in (1 y)~*, s an integer > 0, is (nm), where (mn), =n, and 
If A is any one of y, 8,- - -,7, 
(—1)4(a@— =AI(1— 
in which the coefficient of 2" is A“"-/(n);. Hence 
j=1 
Since A”"* = 1 for all integers m, 
D(an+b)=3> (an + b); T 
j=1 


The expression in square brackets is a constant, say kj. Hence 


D(an + b) = k;(an + b);, 
j=l 


which may be written in the form (1.3), since (an + 6); is a polynomial of 


degree j} —1 in n with coefficients independent of n. 


3. An _r-fold Lambert series. (ne possible generalization of a (1-fold) 
Lambert series 
J “n x 


n=1 din 
is suggested by (1.3). If in this series x is replaced by x°, where ¢ is a 
constant integer > 0, the effect is merely that of a change of variable with 
an obvious modification of the convergence criteria.* It will be assumed that 
f(n) is a numerical function of n such that all series in which it occurs here 
are convergent. 


5 As stated, for example, by K. Knopp, Journal fiir Mathematik, vol. 142 (1913), 
pp. 283-315, especially 287, Satz 1. 
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With a,,- --,a,, a, b, D(n) as in and we define 
an r-fold Lambert series L(x; ;/) =L (a) for the numerical func- 
tion f and the parameters a,,° a, 

F(n) 


To expand L(x) in a power series we multiply the generating function 
of D(n) by f(t)x*, replace x by zt in the result, and sum over all integers 


t>0O. Thus 
L(2) = + f(t) 
ler, 


t=1 m=s+1 


in which > refers to all pairs (d,8) of positive integers such that d’ = m, 
d>s. If there are none such for a particular m, S is vacuous and there is 


no term 
We now reduce D(d —s) by (1.3), noting that ¢j in (1.3) depends on 
b, say cj =c;(b), so that 


D(n) = > cj(b)n/, 


J 0 
For 0 fixed let (dy, 8) be a (d,8) such that 
d, =s +b moda, dy >s+b, = m. 


Then (d,—s— b)/a is an integer > 0, and by (1.3), 


r-1 


D (dy — 8s) = go(do) (b) (d, — s — 
j 0 
and hence 


Ef(8)D(d—s) 
Finally, then, the complete expansion is 
L(x) > [f (m/s) +S dy) f(8,) |x”, 
in which f(y) = 0 if y is not an integer > 0. 
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ON THE BLOCH-LANDAU CONSTANT.* 


By Hans RADEMACHER. 


1, In 1929 Landau, in a paper? on the Bloch constant, introduced 
besides the Bloch constant B two other constants 2 and & of which only the 
first will interest us here. The definition of the Bloch constant is the following: 

There exists a constant 8 such that to every « > 0 and to every function 
w= f(z) which is regular for | z| <1 with f’(0) —1 there exists a closed 
region @ = (;,. in the interior of the unit-circle which is mapped by the 
function w= f(z) in a one-to-one manner on a circle of radius 8 —e in the 
w-plane; whereas there exist functions of the same specifications by which no 
subregion of the unit-circle is mapped in a one-to-one manner on a Circle of 
radius B+ «. 

By dropping the condition of bi-uniqueness Landau defines another con- 
stant which we shall call the Bloch-Landau constant 2: There exists a constant 
2 such that every function w=f(z) which is regular in |z| <1 with 
(0) = 1 assumes ail values of a certain circle of radius £ —e in the w-plane, 
for any given positive «; whereas to every « > 0 there exists a function of the 
same specifications whose values do not completely fill any circle of radius 

Since every “ Bloch circle” in the w-plane is a “ Landau-circle,” but not 
conversely, it is clear that we have @®=&. Landau, in the paper cited, 
proves also 

0.39 <B, 043 <2 S 


where the latter upper bound is established by a certain example of f(z) 
which does not permit a greater Landau circle. Landau gave no separate 
upper estimate for 8. This was done some years ago by Ahlfors and Grunsky 


who proved * by a suitable example that 
B= (1 + 3%) = 04719- 


I shall show that their method leads also to an improved upper estimate ot 2%. 


* Received March 21, 1942. 

1E. Landau, “ ber die Blochsche Konstante und zwei verwandte Weltkonstanten,” 
Mathematische Zeitschrift, vol. 30 (1929), pp. 608-634. 

*L. V. Ahlfors and H. Grunsky, “ Uber die Blochsche Konstante,” Mathematische 
Zeitschrift, vol. 42 (1937), pp. 671-673. 
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2. The interior of the unit circle | | <1 is mapped on an equilateral 
equiangular circular triangle of angles ar (0=a< 1) by means of the 


function 


1 


2-a — t)~°/8+a/2 Ct) -1/6+a 


+ 34/2) , F(%- 2, 4 — 4/2, 


where F(a,6,c;z) is that branch of the hypergeometric function which is, 
for |z| <1, represented by the Gaussian hypergeometric series. In_ the 
mapping established by (1) the points 1, p,p? (p= e°™/*) of the ¢-plane 
correspond to the vertices of the circular triangle in the w-plane. 

The formula (1) is given by Ahlfors and Grunsky for «=1/hk, k >0. 
But it can be seen that the arguments which they set forth to establish its 
validity hold also for «0. The case «= 144 means the mapping of the 
unit circle on an equilateral straight triangle. We can, therefore, by com- 
bining two conformal mappings of the sort (1), map a zero-angled equilateral 
triangle on an equilateral straight triangle. If we introduce the inverse 
function 

f= 
of w= then 
(2) w=f(z) = C1173 (Po(C22) ) 


will produce that conformal mapping; the constants C; and C. can be used 
to adjust the sizes of the triangles. 

Now the conformal mapping (2) can be extended by reflection with 
respect to the sides of the triangles. The repeated images of the zero-angled 
triangle in the z-plane will just fill its orthogonal circle, which is its circum- 
circle, whereas the repeated images of: the equilateral triangle in the w-plane 
will build up a Riemann surface in which the ‘vertices of the triangle and of 
its reflections form a regular point lattice, and where all these vertices become 
branch-points of infinite order. The function w= f(z) therefore maps, by 
this process of analytic continuation, the interior of a whole circle K on the 
described Riemann surface. We can normalize the mapping 'so that K’ is the 
unit circle | z| <1 and so that f/(0)—1. 

Since, then, the image of any circle | z | = 1—e will not cover any of 
the branch-points in the Riemann surface over the w-plane a Landau circle 
cannot be greater than the circumcircle of the equilateral straight triangle 
which was the element out of which the Riemann surface was constructed. 
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ON THE BLOCH-LANDAU CONSTANT. 


The normalization just mentioned requires 


= 
or 
(3) Cz = do(1) 
and 


$'o(0) 


Equation (3) has to be understood as 


lam? 


lim = lim iz £5, 

Here we have hypergeometric functions F(a,b,c;é) with a+b—=c, in 


which case 


F(a,b,c;&)~ 
r(a 


Hence 


and therefore formula (4) furnishes 
$'(0) r(% 
3(0) r(%) 
The upper bound of the radii of all Landau circles, i.e. the cireumradius of 
the equilateral triangle in the w-plane, is 
f(1) = 3(1) = C1 = (% 
We have therefore our result: 
An upper bound of the Bloch-Landau constant x is given by 
(5) (14)T(%)” 
This result compares numerically with Landau’s estimate quoted above as 
follows : 
0.54325 < < 0.0. 54326 


0. 55488 < (74) * < 0.55489. 


“In a letter of May 30, 1942, Professor R. M. Robinson informs me that he found 
the same upper bound for @ in 1937. His result is briefly mentioned in a footnote to 
Ahlfor’s paper “ An extension of Schwarz’s Lemma,”* Transactions of the American 
Mathematical Society, vol. 43 (1938), p. 364 which I had overlooked: “In the other 
direction R. M. Robinson has proved Q < 544. This result has not been published.” 
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3. The configuration of the zero-angled circular triangle together with 
‘all its reflections is, of course, the well-known modular figure; and if it had 
not been for the convenience of the Ahlfors-Grunsky formula, I would have 
adhered to my first redaction of the present paper, which was based on the 
modular function k?(r) =A(r) into which ®)(w) can be transformed. 
Landau in his paper quotes previous results of Carathéodory and Hartogs, 
which are also obtained from the theory of modular functions. An analysis 
of Landau’s procedure shows that he has actually used the conformal mapping 
of the zero-angled triangle with the vertices 1, 1, —1 on the equilateral 
triangle. It is through a gain in symmetry of our example over Landau’s 
that we obtained the improved estimate (5). 

Ahlfors and Grunsky in their paper make it plausible that their appraisal 
of Bloch’s constant actually gives its true value. A similar reasoning would 
apply to the above example. The symmetry of the discussed conformal map- 
ping suggests strongly that (5) may contain the best value of the Bloch- 
Landau constant &. 

It is noteworthy that the Ahlfors-Grunsky example, as well as ours, 
exhibits functions which have the unit circle as natural boundary. This is in 
agreement with R. M. Robinson’s statements about “ Bloch functions,” * 
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‘Raphael M. Robinson, “ Bloch functions,’ Duke Mathematical Journal, vol. 2 
(1936), pp. 453-459. 


THE RECTILINEAR MOTION OF A GAS.* 


By Monroe H. Martin. 


Introduction. [n the rectilinear motion of a gas the particles of the gas 
move along straight lines. The lines are all parallel to each other and it is 
assumed that the velocity and density of the gas at any given instant is con- 
stant over any plane perpendicular to the straight lines along which the motion 
takes place. Letting « denote the distance along any one of these lines 
measured from some point as origin, and denoting the time by ¢, the physical 
problem we consider may be formulated as follows. Given the velocity u and 
the density p of the gas as a function of x for t=0, what can be said about 
the subsequent behavior of the gas, 1. e., how shall u,p be expressed as func- 
tions of x,t? In addition to assuming that the gas moves rectilinearly we 
shall suppose that it is perfect and that its pressure and density are connected 
by a characteristic equation p= p(p). 

This type of motion has been studied by Riemann in a celebrated memoir.’ 
In it he reduces the problem of determining the motion of the gas under 
certain conditions to finding a solution of a linear, hyperbolic partial differ- 
ential equation of second order of the type wrs—a(r,s)(wr—we) = 0, 
taking prescribed values for its first partial derivatives along a given arc.” 
The solution of this initial value problem given by Riemann has become 
classic and is known today as “ Riemann’s method.” LKssentially the method 
consists in determining a two parameter family of solutions of the adjoint 
partial differential equation subject to certain initial conditions. This two- 

* Received February 21, 1942; Revised August 22, 1942. 

‘“ ber die Fortpflanzung ebener Luftwellen von endlicher Schwingungsweite,” 
Géttingen Abhandlungen, vol. 8 (1858-59), pp. 245-264 or Werke, Leipzig (1892), 2nd 
Edition, pp. 156-175. See also H. Lamb, Hydrodynamics, 5th Ed., Cambridge 1924, 
pp. 456-459, where both the methods of Riemann and Earnshaw are discussed. For 
further work on the problem see Rayleigh’s Theory of Sound, vol. 2, p. 38 and vol. 3, 
pp. 214-216 of Durand’s Aerodynamic Theory. An excellent bibliography on the 
problem is given in the article on Compressible Fluids by H. Bateman in Bulletin 84 
of the National Research Council, pp. 543-551. An admirable treatment of the subject 
will be found in Die partiellen Differentialgleichungen der Mathematischen Physik 
prepared from Riemann‘s lectures by Heinrich Weber, Braunschweig 1919, vol. 2, 
pp. 503-562. 

An error in §4 of Riemann’s memoir is pointed out in the report of Christoffel, 
Fortschritte der Physik, vol. 15, p. 123. For comments on Lord Rayleigh’s objection 
to Riemann’s work see pp. 559-562 of Weber’s book cited above. 

* Werke, pp. 160-162 and pp. 169-170. 
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parameter family of solutions is the so-called “ Riemann’s function” and, 
once it has been obtained, the determination of the solution of the original 


initial value problem is reduced to a quadrature. Since Riemann’s time his 
method has been extended to the most general linear hyperbolic partial 


differential equation of the second order.* 


Sections 1 and 2 of the present paper contain no new results but are 
devoted to a resumé of certain results * due to Riemann which the author has 
felt advisable to include in order to properly prepare the ground for the new 
developments which follow. 

In section 3 the soiution of the initial value problem mentioned above 
is reduced to the determination of what is termed the resolvent of the partial 
differential equation. The resolvent is analogous to Riemann’s function and 
is a two parameter family of solutions v = O(r,s; 7,5) of the partial differ- 


ential equation v;s + a(r,s) (v,— vs) = 0 satisfying the initial conditions 


As in Riemann’s method, the solution of the initial value problem is reduced 
to a quadrature, once the resolvent is obtained. The new method offers two 
advantages. The initial conditions on the resolvent are simpler than those 
for Riemann’s function and the final expression for the solution as a quadra- 
ture is simpler in form. On the other hand the method does not appear to 
permit extension to the general linear, hyperbolic partial differential equation. 

The method developed in 8 finds an application in 4 to the solution of 
the problem of the rectilinear motion of a gas with which Riemann occupied 
himself in § 10 of his memoir. The culminating point of his work in this 
paragraph is contained in his equations (2). These equations involve a sune- 
tion w of four arguments r, s, 7’, s’.. However we do not find any formula in 
iiemann’s paper in which the nature of the dependence of » upon its four 
arguments is explicitly set forth. 

In the theory developed in the present paper the resolvment 7 plays the 
role of the function » of Riemann and in 5, 6 we calculate 7 explicitly for the 
important case in which the gas is moving adiabatically. In the determina- 
tion of 7 we are led to the generalized integral equation of Abel, the ‘solution 
of which yields a formula for 7 involving Appell’s hypergeometric function 
F,(«;B, B’sy; x,y) of two variables. 

We have selected in 7 what appears to be mathematically the simplest 


*Cf., for example, D. Hilbert und R. Courant, Methoden der Mathematischen 
Physik, II, Berlin (1937), pp. 311-315. 
* Werke, pp. 168-175. 
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problem to exemplify the application of the methods of this paper to the 
actual computation of the rectilinear motion of a gas. The application of 
the method to more complicated cases in which a shock wave develops is 
reserved for a later paper. 


1. The partial differential equations of the problem. If it be assumed 
that the particles of a gas move in straight lines parallel to the z-axis of a 
rectangular codrdinate system in three dimensional Euclidean space and that 
the velocity and density of the gas are constant over any plane 2 = const. at 
any given instant ¢, the partial differential equations for the velocity 

=u(x,¢) and density p= p(z,t) are® 


(1) p(t: +- + G*p, = 0, pt + pret + pur = 0, 
it being understood that there are no external forces. Here 
(2) G=G(p)=(p’)*, p=p(p), 


the latter being the characteristic equation of the gas,° so that G is the “local 


velocity of sound 


(3) f=f(p) = fGdp/p, 


is introduced in place of p, the partial differential equations take the more 


If a new dependent variable f, defined by 


symmetric form 
uutGfe=0, frt+feu+ Gus—0, 


which, in turn, may be replaced by 


(4) r+ (u+6)r.=0, St + (u— = 0, 
on taking 
(5) u=rt+s,f=—=r—s; 2r=u+ f, 2s—u—f. 


Thus u and p are functions of r+ s, r—s respectively. Therefore G, 
a function of p, is likewise a function of r—s and we have 


(6) ur = Us == i. G, = — Gs. 


Let us consider a solution 


~ 


( ) re=1(z,t), s=s(z,t), 


* Literal subscripts denote partiai differentiation. Equations (1) follow im- 
mediately from the hydrodynamical equations and the equation of continuity on setting 
Z 0, v= w=—0, = 0, = 0. 
*The prime denotes differentiation with respect to p. It is assumed that p(p) 
is a monotone increasing function of p possessing a continuous derivative with 
Tespect to p. 
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of (4) defined and of class C’ throughout a region X of the (z,¢)-plane and 
such that the Jacobian 


0(1, 8) 
J = | == 2G 
t) 
does not vanish in X. In addition we suppose the region X is mapped ina 
one-to-one manner upon a region R in the (r,s)-plane by (7) and denote by 
(7)’ (r,s), t = 8), 
the inverse transformation from R to X. Now 


JI Xe. => Jt,. st = J Zr. 8; >> 


and substituting for r+, rs, s:, sc in (4) we see that the functions (7)’ are 
solutions of two partial differential equations, linear in the unknown func- 
tions x, t 

(4)’ t,-— (u—G)t, = 0, (u+ G)t,=0. 


A solution (7) of (4) enables one to compute the state of the gas, that is, 
the velocity and density of the gas, on a plane & at time /; provided, of course, 
that the point (#,/) lies in XY. On the other hand, the solution (7)’ of (4)’ 
enables one to compute the location and time of occurrence of a preassigned 
state; provided, naturally, that the point (7,5) corresponding by (5) to the 
state in question lies in R. Either aspect is accordingly available for the 
description of the motion of the gas. 

Choosing the functions z, ¢ in (7)’ to be solutions of (4)’, the results in 


(6) may be used to infer that the line integral 
w= f{[x—(u+ G@)t]dr+ [a — (u— G@)t]ds}, 
is independent of the path. Moreover, if w=w/(r,s) is the function of r 
and s thereby defined, we have 
«—(u+G)t—v,, 
Subtracting the two equations, and differentiating the first partially with 


respect to s, we obtain 


(8) Wy -— Ws, = — 264, Wre = (G,—1)t, 


and w is accordingly a solution of the partial differential equation of the 


second order 
(9) (wy Ws) = (), 


Conversely, if w is a solution of (9), equations (7)”, when solved for 
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and a,t yield functions (7%)’ of r,s which are solutions of (4)’. This follows 
directly upon differentiating the first equation in (7)” partially with respect 
to s and the second equation partially with respect to r, the relations (6) 
being kept in mind. The integration of (4)’ is therefore reduced to the 
integration of (9). 
» ie Differentiating the first equation in (7)” partially with respect to r, the 
; second partially with respect to s, and using (4)’, one obtains 
Eliminating ¢ by means of (8), the Jacobian 
s) 
unc- . . 
is expressible in terms of w as follows, 
I 1 + G, 1 Gy Gr 
(w, — we) | [Wes — (wr — Ws 
t is, 
Tse, By hypothesis the Jacobian J is finite and not equal to zero in the region 
(4) X. It follows that the Jacobian J-* is not equal to zero and is finite in R. 
red Thus if w denotes the solution of (9) corresponding to the solution (7) of 
the (4), the region RF contains none of the loci 
the 
] + G; 
— Ws ) == ¢€, W as | (u == €, (e 0 or 
in 
2. The initial conditions. If the region X over which the solution (7) 
of (4) is defined is assumed to contain the segment a << x < B of the z-axis, 
this segment is transformed by (7) into an arc in the (r,s)-plane 
of (10) Co: r==1o(z) = r(2z,0), =s(2z,0), 
which, since J cannot vanish in X by hypothesis, is met by a line r = const. 
(or s const.) in at most a single point, or, in other words, Cy is a mono- 
with tone are. 
Thus the initial velocity vw) and initial density po of the gas are assumed 
to satisfy the inequalities 
Po Po 
along the segment « << « < B of the z-axis. Conversely, if the above inequali- 
ties are met, the arc Cy is a monotone arc. 
One implication of the assumption that the Jacobian J does not vanish 
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may now be pointed out. If the gas is initially at rest, the density must be a 
monotone function of z and the case of wave motion is ruled out.’ The case 
in which the Jacobian vanishes identically has been discussed by Riemann ® 
and various other writers. 

In view of the preceding results, in particular (7)”, the problem of con- 
structing a solution (7) of (4) subject to the initial conditions (10) is equiva. 
lent to the following boundary value problem: to construct a solution w of (9) 
such that the first partial derivatives take values prescribed by 
(11) Wy, = 2, We =X, 
on the monotone are C4. 

3. Conjugate partial differential equations. Two partial differential 
equations of the type 
L(u)= — — ty) = 0, M(v)=v.y,+ — vy) = 0, a—a(z,y), 
are said to be conjugate to each other. Taking u, v to be solutions of L(u)= 0, 
M(v) =0 respectively, it is readily verified that 


(UsV2)y = (Uyvy) 2, 
and the integral 
({u.vedr + uyrv dy}, 


taken around any closed path is accordingly equal to zero.° Applying this to 
the closed path formed by two straight lines e = Z, y= ¥ and a monotone are 
XY connecting a point X of z= to a point Y of y= ¥, one has 


P 
Uyvydy + f {urvedz + Uyvydy} + f = 0, 
P x 


7In case the density is not a monotone function of 2, one may apply the method 
to the intervals in which the density is monotone. However this will not suffice in 
general to define the motion of the gas for all a, t. An interesting problem of ap- 
parently some difficulty arises when one seeks to prolong the solution over the entire 


t)-plane. ; 

® Werke, pp. 162-164. Durand, Aerodynamic Theory, loc. cit. See in particular 
H. Weber, op. cit., pp. 516-522. 

® Naturally, this statement may require qualification if the closed path bounds a 
region containing points at which U,, Uy, Vz, Vy possess singularities. The referee 
points out that the relation of the method used in the text to classical theory is brought 
out by the identity 

(v,— v,)L (oe) + (6, — u,,) M(v) = P=4u,v,, Q=u vy,» 
analogous to the well-known identity 
(uw) —uM(v) = 

where M is adjoint to Z. Assuming that L(u) = M(v) =0, it is apparent that 
§(%e%2t” + u,,v,dy} taken around any closed path is zero. 
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bea where P is the point (#,%). Consequently if a solution 7 of M(v) =0 be 
case taken so that 
inn (12) on PX: on PY: — f, 
ois one finds the value u(P) of the solution u at the point P expressed in terms 
del of its values at X and Y, and the values of wz, uy along the monotone arc 
| (9) YY, by the formula ; 

(13) u(P) =4[u(X) + u(¥)]—4 + uydydy}. 
The solution 7 of M(v) 0 employed above obviously depends on the 
atial coordinates 7, Y of the point P, a fact which we signalize in the notation 
T=0(z,y; P). 
9), The function 7 is termed the resolvent of L(u) =O and the problem of 
= 0, integrating L(w) = 0 has been reduced to its determination. 
The resolvent 7 of Af(v) =0 is a solution of L(w) =0 equal to y* —y 
on 2* and to on y=y*. It is connected with the resolvent 
of L(u) = 0 by the relation 
(14) i(P; P*) +3(P*; P) =0, 
a analogous to the Principle of Reciprocity in the theory of the adjoint partial 
| differential equation. To prove this relation, consider the parallelogram with 
the four vertices 
P(z,9), P(2*,G), P*(x*,y*). 
| Replacing w by @ and the are TY in (13) by the broken line XP*Y, one finds 
sthod ,P* 
cea u(P; P*) =43[a(X; P*) + P*)] f ii2(P; P*)t2(P; P)dx 
fap Xx 
ntire 
and, since by definition 
ids a 
feree P*) =1 on P*X, iy(P; P*) on P*Y, 
vught 
the relation is established, inasmuch as 
a(P; P*)—y* —§, V(X; P)—9—y*, ; P) a* — 
aid Among the relations to be obtained from (14) by partial differentiation 


with respect to 7, 9, x*, y* we shall set down only two 
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(15) P*) + P) =0, 


as important for later application. 
Returning to the solution (13) and differentiating it partially with 

respect to ¢ and g, it is found that 


(16) uz = {u,Trzdx + u,tyzdy}, 


4 
17 = + u,dyzdy}. 


It will be recalled that (14) was derived from (13) by replacing wu by @ and 
the are YY by the broken line XP*Y. If the same procedure is applied to 
(16) and (17), it is found that 


P) = 2[as(P; P*) + P)], 
(¥; P)= 2[az(P; P*) + P)], 


+ 
inasmuch as 
(18) tz(Y;P) =—1, t7(X; P) —1. 
Referring to (15) it is apparent that 

=0, 
which, since may be taken arbitrarily on c= and Y may be taken 
arbitrarily on y=¥, leads to the following principle. The _ resolvent 
satisfies the relations 


(19) t+%—0, 
respectively along the lines x = and y= ¥. 

4. Rectilinear motion of a gas. In the first two sections of this paper 
the problem of determining the rectilinear motion of a gas under the hy- 
pothesis stated therein has been reduced to the boundary value problem 
formulated at the close of 2. For any equation of state p= p(p) the partial 
differential equation (9) takes the form 


(20) a(wy—w,) = 0. a=a(r,s), 


and, granting that the resolvent, 


T=T(r,s; 7,5), 


ith 
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of this partial differential equation has been determined, the value of the 
solution w at a point (7,5) is, from (13), computed by 


w(7,5) =4[w(k) + f {w,0,dr + 


in which R, S denote the points where '° the lines r=, s = §, respectively, 
intersect the monotone arc Cp. 

What is really needed, however, in (7”), is not the value of the function 
w at the point (7,5), but the values of its partial derivatives w,, w, at this 
point. For these we refer to (16), (17) and obtain 


8 

w7 = 4w,(R)[1+7,-(R)] {w,t,7dr + 
R 

wz = —4 f {w,t-gdr + 
R 


which, when the initial values (11) for w,, ws on Co are inserted, become 


R R 


where «(R), denote the abscissae of the points on the interval a <x < B 
of the z-axis which map into the points R, S under the transformation (10). 
After integrating by parts and applying (18) and (19) it is found that 


we = 4[2(R) +2(8)]+3 trdx, +2(8)] +34 f 


2(R) 


The motion of the gas is consequently to be determined from the fol- 


lowing equations 


(21) (a G) == v(R) «(8) ] +4 
x(S) 
G)t=}[a(R) + 2(8)] +3 
a(R) 


obtained by substituting for the partial derivatives in (7)” 


*° The values of w at the points R, S are indicated by w(R R), w (8) respectively. 
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The use of these formulae may best be understood with the aid of Figure 1, 
Suppose the initial velocity and density wo, po of the gas are known along the 
segment a <2 < f of the z-axis. To fix the ideas, let us suppose that along 
this segment Uo, po are chosen so that the are Cy of 2 is the monotone are 
ARSB of the figure. The basic idea in the use of these formulae is that they 
yield a mapping of the “triangular region” ABC of the (r,s)-plane upon 
the “triangular region” aBy of the (x,¢)-plane. Thus the mapping brings 
to each point (Z,/) inside «By a pair of values *,5 from which the velocity i 
and density p of the gas at the point (Z,7) is computed from (5) and (3). 

The details of the application of these formulae will be taken up in 7 


where they are applied to a specific problem in the rectilinear mution of a gas, 


u 


Figure 1. 


5. The adiabatic case. The results obtained so far are valid inde- 
pendent of the choice made for the function p(p) expressing the dependence 
of the pressure of a gas upon its density. From now on we shall suppose 
that the motion of the gas proceeds adiabatically. The characteristic equation 


of the gas takes the well-known form 
k*p7, 


where y is the ratio of the two specific heats of the gas and k is a constant. 


It is readily verified that 


G = (r—s), 


and that the partial differential equation (9) becomes 


3 y 


(22) (7 — 8) — m(w, — ws) = 0,7 == 
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This type of partial differential equation has received intensive study ** since 
the time of Euler. 
Equations (21) controlling the motion of the gas take the form 


and, once the resolvent T of (22) has been found, may be used to compute 


8 
si 
& 


+ 2(S)] +3 


8 
& 


the motion of the gas as explained in 4. 


According to the kinetic theory of gases ** 


where n = 0, 2 for monatomic and diatomic gases respectively, a result which 


is in good agreement with experiment. The corresponding value of m is 


m=4(n-+ 2), 


so that m= 1 for monatomic gases and m= 2 for diatomic gases. There 
appears '* to be no known law for n in the case of more complex gases, and 
we are accordingly led to consider the problem of integrating (22) for arbi- 
trary constant values of m. 
The general solutions of (22) for monatomic and diatomic gases are 
R—S R—S 


== 


r—s’ 


where R denotes an arbitrary function of r and S an arbitrary function of s. 
Although the general solution of (22) is available’! for arbitrary constant 
values of m, our problem is to find a solution subject to the initial conditions 
(11). Presumably one might achieve this by properly specializing the arbi- 
trary elements in the general solution. However it appears more feasible to 
proceed along the lines of 4 and we accordingly take up the problem of 
determining the resolvent of (22). 


6. The resolvent. The resolvent 7 of (22) is a two parameter family 


=T(r,s; 7,5) of solutions of the conjugate equation 


** See, for example, G. Darboux, Legons sur la théorie générale des surfaces, vol. 2, 
2nd Ed., Paris (1915), pp. 54-70. 

“J. H. Jeans, The Dynamical Theory of Gases, 4th Ed., Cambridge (1925), pp: 
191-193. Recent work by D. H. Andrews has enabled fairly accurate estimates of 7 
to be made for molecules containing up to 20 atoms. 


| 
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(24) (r— Urs + m (vr —v,) =0, 
such that 


It is easily verified ** that a one parameter family of solutions of (24) 


is given by 


v= (r—a)™(a—s)”, 
a being the parameter of the family. This one parameter family offers the 
starting point for our construction of the resolvent. 

To begin with we observe in (3) that f is a positive function of p and 
therefore it follows from (5) that r>s. We may accordingly confine our 
attention to the half plane lying under the line r=s. Taking 7 > 5, we 
shall construct the resolvent of (22) under the following restrictions 
(26) >? 35 m = A < 
Eventually we shall show how these restrictions may be removed. 

By forming the solution * 

v= f (r —a)*(a—s)>da, 
Js 
where (a) is an arbitrary function we gain a solution of (24) containing an 
arbitrary function. Taking for we arrive at a two- 


parameter family of solutions 


(27) v= + (r—a) (a s)~\da, 


? 
involving an arbitrary function. If this arbitrary function can now be de- 
termined so that the conditions laid down in (25) are met, the resolvent of 
(22) will have been determined. 

Thus we are to determine ¢(@) so that 
er 
(28) f (a— (r —a)*da = r o(a) (F 
r 8 
Let us consider the first of these integral equations. If we define 
y(t) = W(t) =0, y(t) V(t) =t—F, 


this integral equation is recognized as the generalized integral equation of Abel 


= 


18 Cf. G. Darboux, op. cit., p. 58. 
14 Cf. G. Darboux, op. cit., pp. 66-68. 
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The solution of this integral equation is *° 


ea a 
0 r 


from which it is readily ascertained that 


sin 7A 


(29) o(a) = (a—7)(a—S). 


The second integral equation goes over into the first when the variables 
in the pairs (7,s), (7,5) are interchanged. It follows that the solution is 
also given by (29). 

When ¢(a) is substituted in (27) and a new integration variable ¢ is 
introduced by setting 


a=fr+t(r—?), 


in the first and second integrals respectively, the resolvent 7 is obtained in 
the form 


t\[#F — § — t(s—- 


The resolvent will now be expressed in terms of Appell’s hypergeometric 


16 


function of two variables 


OO / 
(a, m + n)(B,m) nr) 
(99) BPG (y,m + n)(1,m) (1,7) 


Indeed if R(a) > 0, R(y— 2) > 0, Picard has shown ** that 


1 F , 
a-i(] tr)-8 (1 — ty) = (a)T(y 2) F(a; B, B’ sy 


If we put 1+A, B=—A, B’ =A, y= 2, in this formula and observe that 
sin 7A ’ 


we find on writing 


“The solution of this integral equation is well known. Cf., for example, G. 
Kowalewski, Integralgleichungen, Leipzig (1930), pp. 11-13. 

‘*P. Appell, “Sur les fonetions hypergéométriques de plusieurs variables,” 
Memorial des Sciences Mathématiques, Fase. III (1925), p. 4. 


*7 See Appell, op. cit., p. 9. 
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r—?7 s—$ s—$ 


that the resolvent may be expressed as follows 


3 


s—S§ 


It will be recalled that the resolvent was derived under the restrictions 
(26). These restrictions we now proceed to remove. Availing ourselves of 


the formula 


H 


we find alternative expressions for the hypergeometric functions in (31), 


namely, 


r—sS S—S s—s 
( yr (1. A3;—A,A;2;3 ), 
r—s s—f* s—r 


and thus we are led to consider three additional forms for T: 


eo 

zy 
> 

| 
> 

| 


(33) T= (=) F, 


r—f 
(34) (7 (=) € A; A,A, 
S—S s—S 


18 See, for example, W. N. Bailey, ‘‘ Generalized hypergeometric series,” Cambridge 
Tracts in Mathematics and Physics, No. 32, Cambridge (1935), p. 78. 


| 
r—s§\A 
r—s r—s 
PF, 
AY 
( ;) s—S§ s—S 
s—F 
r—s s—f’ s—r 
g 
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Appell’? has shown that the power series (30) converges for arbitrary 
real or complex values of a,B,y inside the unit square |2| <1, |y! <1, 
provided, of course, y is not a negative integer. The various formulae obtained 
for T accordingly define the resolvent in the following regions respectively : 


31)’ |r—#l <|#—S!,|r—F| <|s—F|,|s—s| < |F#—S|, |s—S| <|r—S]; 


2)’ |r—#| <|r—s|, |r—#F| <|r—s|, |s—3| < | |, | |; 


( 
( 
|r—#| <|#—S|, |r —#| <|s—#], |s—5| <|s—F|, 
(34)' | r—#| <|r—s|, <|r—s|, |s—8| <|s—#|, |2—5| 


The four forms (31), (32), (33), (34) obtained for @ are analytic con- 
tinuations of each other and taken together suffice to define the resolvent 
throughout the region 

Ss, 


for any real value of A. By hypothesis the are RS in (21) or (23) connecting 
a point of r= 7 to a point of s=S is a monotone arc lying in the region 
r>s and therefore cannot leave the above region. It follows that 7 is ade- 
quately determined for use in (21) or (23). 


0 


The formula ° 
F, (a; 32, 2) = (4,8+ 
reducing the hypergeometric function F, to an ordinary hypergeometric func- 
tion affords a partial check on our computation, inasmuch as, taking r=7 


and s== § in (31) yields 


= (§—s)F (1 + A, 0, —5—s, 


57,5) = (r—*)F (: + 2,0, 2; 


which agrees with (25). 


A further check is provided by the monatomic gases, the partial differ- 


ential equation (22) becoming 


(35) (7 —S)Wre — (Wr — We) = 0. 
Taking A= -— 1 in (31), the resolvent of (35) is calculated with the aid of 
(30) and found to be 
(36) (r —F¥)(#¥—s) — (s— S$) (r—S) 


° Op. cit., p. 4. 


20 P 


Appell, op. cit., pp. 3-4. 
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and a direct computation verifies that this function 7 is a solution of the 
conjugate equation to (35) and that 7 fulfills the conditions required in (25), 


7. The motion of a monatomic gas. Of all the gases the motion of a 
monatomic gas appears to be most susceptible to mathematical investigation. 
Thus the resolvent T has the particularly simple form (36). One readily 
verifies that placing this resolvent in (23) with y = 5/3 leads to the following 


explicit formulae: 


a(R) 
t == — $)-*{75[2(R) — x(S)]- f [rs —4(F + 5) (r+ 8) 
2(R) 


To illustrate the application of the method to a concrete problem let us 
investigate the rectilinear motion of a monatomic gas which initially is at 
rest everywhere, a vacuum (p = ()) being present for negative x and 
p= (2a/15k*)*/* for positive x, under the assumption that the gas moves 
adiabatically. The apparent artificiality in the choice of the initial density 
due to an attempt on our part to find an example which illustrates the method 
without burdening the reader with tedious and excessive computation. 

Referring to 1, one verifies that f= 15 kp’/* and that the parametric 
equations (10) of the are C, in 2 (the arc ARSB in Figure 1) are 


== (2/2), s= — (2/2). 
Substituting «(R) = 27", 2(S) = 25° and carrying out the integration in the 


above formulae, we find on dropping the bars from all letters that the region 
r = s of the (r,s)-plane is mapped upon a region in the (2,¢)-plane by 


== — (r? + 4rs + 8°), 3(r+s). 


These formulas take on greater physical significance if the variables u,f 
are employed in place of r,s. In terms of these variables, the mapping is 
given by 
(37) x= 4(f? —3u’), t = — 3u. 


It is apparent that the loci in the (z,¢)-plane of points where the gas has 
constant velocity comprise the lines ¢ = const., while the loci of points where 
the gas has constant density are the parabolas 7 = 4$(f? — ¢?/3), f = const. 
These curves are graphed below and form a network over a region in the 
(x, t)-plane, from which it is possible to read off the velocity and density of 


| 
(F + FS + 45°) a(R) — 5(4F + + 
2(S) 
[3(7 + 5)rs —2(F? + + +58) ]dz}, 


ric 
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the gas at any position and at any time. Outside this region both the density 
and velocity of the gas is zero, i. e., physically one has a vacuum. As a check 
on our computation, one finds on solving (37) for u,f in terms of z,t, that 


15k? 
t 
A 
u=-5 
u=-2 
\ 
Sunt” 
= 
u=6 


Figure 2. 


yield a solution of the hydrodynamical equations (1). It will be noted that 
a discontinuity in the velocity presents itself along the parabola p—0 of 
Figure 2. The discontinuity meets, however, the compatibility, impulse, and 


21 


energy conditions. 


UNIVERSITY OF MARYLAND, 
COLLEGE PARK, Mb. 


** For a formulation of these conditions, see Weber, op. cit., p. 511. 
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A CLASS OF SOLUTIONS FOR THE HEAT EQUATION AND 
ASSOCIATED BOUNDARY VALUE PROBLEMS.* 


By F. G. Dresset and E. R. Evurorr. 


1. Introduction. Associated with the heat equation 
(1) 0°u/dx? = du/dy, 


we have the region D bounded by the characteristics whose ordinates are h 


and g (h < g), and the curves 
t= x2(y), hSySq, Xi < Xz 


The problem of finding a solution of (1) in D, taking on preassigned con- 
tinuous values on the boundary of D* has been treated by several authors. 
Holmgren * solved the problem, using integral equations, under the assump- 
tion that the yi (7) had continuous derivatives. Later papers have sought to 
weaken the hypothesis on the yi(7)* or to use other methods for the solution.‘ 

In the first part of this paper we treat a problem similar to that of a 
former paper,’ where certain line integrals of solutions of (1) were made to 
take on boundary values of limited variation.® In Part 2 of the present paper, 


* Received June 19, 1941; Revised September 18, 1942; Presented, in part, to the 
Society, April 8, 1939. 

1 No boundary values are assigned on the line y = g. 

* E. Holmgren, “ Sur une application de l’équation intégral de M. Volterra,” Arkiv 
fér Matematik, vol. 3 (1906-07), pp. 1-4. 

3M. Gevrey, “Sur les équations aux dériveés partielles du type parabolique,” 
Journal de Mathématiques, (6), vol. 9 (1913), pp. 305-471. Gevrey showed that, 
if the x;(7) satisfy a Lipschitz condition of order a > 4 the problem has a solution. 

“W. Sternberg, “ Uber die Gleichung der Wiarmeleitung,” Mathematische Annalen, 
vol. 101 (1929), pp. 394-398. Here the problem was solved by Perron’s method under 
the assumption that the derivatives of the x;(7) were bounded. I. Petrowsky [“ Zur 
ersten Randwertaufgabe der Warmeleitunggleichung,” Compositio Mathematica, vol. 1 
(1934-35), pp. 383-419] improved Sternberg’s treatment and showed that a weaker 
condition than that of Gevrey was sufficient. 

5F. G. Dressel, “A boundary value problem for the heat equation,’ American 
Journal of Mathematics, vol. 55 (1933), pp. 641-653. This will be referred to as 
paper (A). 

° The results of Part 1 of this paper are taken from E. R. Elliott’s doctoral dis- 
sertation presented while in the graduate school of Duke University. The investigation 
was suggested by Professor E. R. C. Miles, to whom Mr, Elliott is indebted for sug- 
gestions and encouragement. 
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we show how to construct a solution of equation (1) in D which takes on 
preassigned functions of limited variation on the boundary. However, the 
approaches to the boundary of D are in a restricted sense that we call parabolic 
approaches. 


PART 1 


In this part of the paper we shall show how to construct a solution 
u(x, y) in the region D satisfying boundary conditions of the following type: 


(2) lim dt = B(x.) — B(2,), 5>0, 


(3) lim + dt = — (—1)*A,>0, i= 1,2, 


"1 
where B(a) and Gi(y) are preassigned functions of limited variation. The 
problem differs from that in paper (A) in that here the line integral of uz 
parallel to the curves x; is used in place of the line integral of the solution u 
itself. ‘To carry through some of the proofs we have found it necessary to 
replace the Riemann-Stieltjes integral used in our former paper by the Young- 
Stieltjes integral. 

Before stating our principal result we shall prove several introductory 
theorems about the class of functions we are using. 


2. Lemmas. We first prove a lemma on the change of order of 
integration. 
LemMA 1. Let G(a,y) be continuous in mSaSn, pSySq, with 
the possible exception of the point (a,b), m<a<n, p<b<q. Let 
f G(x, y)dz| <M, 
m 


If F(y) is of limited variation, then 


J de: | y)dF(y) = dF(y) - f. G (2, y) dz, 
m Jp p m 


when the above integrals are considered in the Young-Stieltjes sense. 
From the hypothesis, it follows that ; G(x,y)dx exists for each y, is 


bounded in y and continuous except for y= b; hence 


(4) dF (y) G(2,y)dx 
p m 


exists in the Young-Stieltjes sense. Further by noting that G(2,y) is con- 
tinuous for mS a2Sa—8, pSyXq, (8>0), we have 


4 
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a-5 q a-6 
m p p m 


As §— 0, the integral on the right is seen to have a bounded limit. A similar 
treatment of the interval a= 2=n shows that the integral 


exists. Since the two integrals (4) and (5) exist, they are equal by a theorem 
of W. H. Young.’ 
The fundamental solution of (1) will be denoted by 


U (x,y) = [4ry]+ exp (- =), ¥>0; =0, yS0. 
By V(z,y) we shall mean the function 
(6) =— 2U2(z,y), y>0; =(0, ¥=0, 


which is a solution of (1) for y > 0. 
We shall now prove several theorems about solutions of the heat equation 
constructed from the functions U and V. 


THEOREM 1. Let F(y) be of limited variation on hh Sy and cor 
tinuous at h. Let x(y) be continuous on the: same interval. Then 


h+d 
(8) = f U x(n), h+8—x)df(n) with 8—>0. 
ry h 


I(8) is the line integral of a solution of (1) along the line y=h+é 
Applying our first lemma, we have 


h+d Le 
f dein): ff x(n), h+8—n)dz | 
=M-T(h,h +8), 


where M is an upper bound for the integral of U with respect to « ani 


T(h, h + 8) denotes the total variation of F(y) over hSySy+6 
T (h,h +8) approaches zero with 8, since F'(y) is continuous at y = 


lim [(8) = 0. 


For completeness we state 


™W. H. Young, “Integration with respect to a function of bounded variation; 
Proceedings of the London Mathematical Society, Series 2, vol. 13 (1914), pp. 109-15 
See p. 148. 
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THEOREM 2. Let B(x) be of limited variation on aX 2 Sb, continuous 
at a and b, and have regular discontinuities fora< «<b. Then 


U (« — é,8)dB(é) = B(a2) — B(a,). 


6-0 1 


This follows from Theorem 2 of paper (A). 


THEOREM 3. Let y’(y) be continuous on h Sy Sg, and let B(a) be of 
limited variation on aS and continuous at x= y(h) (aS x(h) Sd). 
We then have 


Lim fa: VOW aB) 


The limit function is absolutely continuous onh= yg. (A is assumed to 
approach zero either through positive or through negative values.) 


We shall first show that s V | dt is uniformly bounded. From the 
identity ‘ 
(a—b)exp[— (a —b)*] (a—c)exp[— (a—4)*] + (c—b) 
exp[— (c —b)?]exp[—2(c — b) (a—c) Jexp[— (a—c)?] 
and the fact that (x(t) —x(h))/(t—A) is bounded, we see that 


(7) | V(x(t)+a—& | SM(t—h) (x(n) |, 


where M is a constant depending on x(t) and ,’(¢t). From (7) it follows 
that the integral 


(8) VQ) at 
7h 


is uniformly bounded for a given x(y). 
From Lemma 1 and a property of Young-Stieltjes integrals we may write 


a 
b 
dB(é) «lim V(x(t) -A—ét—h)dt, 
a Yh 


provided the limit of the integrand exists. For + y(h) and A sufficiently 
mall, V is continuous; this gives us the first of the following results: 


0, 
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lim V(x(t) +A—&t—h)dt —ét—h)di, x(h), 
h h 


The second result, when = y(h), is well known.® Thus the limit of 
the Stieltjes integrand exists ; in fact it is continuous in é except for € = x(/h), 
Since B(€) is continuous at x(h), the + 1 contributes nothing, and we obtain 


4 db 
lim f V(x(t) +A— é,t —h)dB(€&) 
h a 
f Bo "V (x(t) —&t—h)dt=4(y). 
a 


From (8) we may interchange the order of integration in ¢(y). Again using 
(8), we then see that d(y) is the indefinite integral of a summable function; 
hence ¢(y) is an absolutely continuous function. This completes the proof 
of the theorem. 
The function H(y(y),y) appearing in Theorem 4 of paper (A) is a 
Stieltjes integral; we now proceed to change it into an ordinary integral. 


y 
H(z, y) dG (n) : f (t)L (x(t) + d- x(n),¢— n) dt, 
h 
r—=A+y(y), AX. 


We may replace G(n) by G(y) = G@(y) — G(h) without changing the value 
of H(z,y). We assume this done, and we also assume y”(¢) continuous. 


After integrating by parts, we can write H in the form 


H (x,y) =- [2y'(t) dy. 
h On 


a 


On adding and subtracting = [2y’(m)U] to the integrand of the inside in- 


tegral, we obtain 


(9) H(z.y)=— G(n) J" — (x(t) 


y 
(x(t) +A— x(n), t —n)dtdy + 2y’(n)U(x(y) +A— x(n), y— Gl 


Since y” is continuous and since there exists a constant K such that 
!ajexp(—|a|) < Kexp(—|a|/2), 


the lim H(z, y) may.be obtained by putting A= 0 in (9). 


\-0 


8E. Goursat, Cours d@’Analyse, vol. 3, Paris (1927), p. 307. 
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To simplifv the notation in the next theorem, define 
pat) 


(10) R(do(y),¥(n)) = V(o(y) —¥(m),¥— 9) 


SOLUTIONS FOR THE HEAT EQUATION 


t—yn 


THEOREM 4. 
h = q. Then, 


following limit in two forms: 


lim 


A\>0+ 


is either 
(11) =+([G(y 
or 


(12) =+([G(y- 


The representation (11) follows from combining (9) for A= 0 with the 


results of Theorem 


A 0, we may change the order of integration and write 


7h n 


\>02 


f “dG (7) lim V (x(t) +A— x(n), — dt. 
\->0: 


From a theorem similar to one of Goursat’s ® 


Putting the above result in (13), we see that (12) is true. 


3. A boundary value problem. Consider the domain D defined in the 
introduction. We shall refer to the segments of the line yh +8, 8 >0, 
with end points (#,,h +8) and (a2,4-+ 8) as admissible 8-displacements, 
if for sufficiently small 8 the following set of points (z,y) lies in D plus 


its boundary 


* Loc, cit., p. 307. 


Let x” be continuous and G(y) be of limited variation on 


for each fixed y in h<ySq, we can represent the 


t 
f V (x(t) +A—x(),t—n)d@(n) 


t 


h 


0)—G(h)] + x(n), t 


t of paper (A). We now proceed to prove (12). For 
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The curves 
+ Ai, (—1)*A; > 0, 


will be called A-displacements if they lie in D except for their end points, 
which lie on the boundary of D. We shall show how to find a solution u(z, y) 
of the heat equation (1) in D which for admissible 8-displacements satisfies 
relation (2), where B(x) is a given function of bounded variation with 
regular discontinuities and continuous at ca and eb. Moreover, for 
A-displacements, u(x, y) satisfies (3), where the Gi(y) are preassigned func- 
tions of limited variation continuous from the left and continuous at yi (h), 
Set up the function 


b 2 1 
(14) u(z,y)= U (x—-&, y—) dB(é)— n) dF, 
a h 


= (©, y) — U2 y) — Us(2, y). 
We see that u(z,y) satisfies (1) in D if the Fi(m) are functions of limited 
variation. From Theorems 1 and 2, we have for admissible $-displacements 


lim fous + 8)dt = B(x.) — B(z,), 


6-0 1 
if the Fj(y) are continuous at y—h. For u(2z,y) to satisfy all our con- 
ditions, we have thus to show the existence of functions Fj(7) of limited 
variation, continuous from the left and continuous at »—/, and such that 
Ye 
(15) lim f Ux (xi(t) Ai. t)dt = Gi (ye) Gi(y1) 
"1 

for A-displacements. 

We shall assume the y”;(y) continuous, and since B(é) is continuous at 


é=a and é=b, we have from Theorem 3 


lim Ur(xi(t) + Ai, t)dt = di(y), 
b dx 


where the ¢i(y) are absolutely continuous functions of y. 


Putting y, =A and ys = y and making use of result (11) in Theorem 4, 


the pair of equations (15) can be put in the form 
2 

(16) fily) =gi(y) +2 Kis 0) dy, 


where the gi(y) are known functions of limited variation continuous from 


the left and continuous at y= 
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gi(y) = (—1)[Gi(y) — Gi(h) — gi(y)] 
and 
fi(y) = Pily) —Fi(h). 

(17) = (—1)*R(xi(y), xi (n)), 
the function FR being given by (10). 

Since the kernel Ki;(y,) is of the form of a continuous function of 
y and » divided by (y—~vy)!, the system (16) will have a unique solution 
continuous from the left and continuous at yh. It remains to show that 
the fi(y) are of limited variation. To show that the fi(y) are indeed of 
limited variation we proceed as we did in section 5 of paper (A). The only 


step we must supply here is to show that there exist functions Ti;(y, ) 
such that 


f° 
7) 


Comparing the results (11) and (12) of Theorem 4, we see that we can 
take the 
(18) (ys 9) = V(xi(y) —xi ¥—1)- 


We collect the principal result of Part 1 into the following theorem: 


THEOREM 5. Let B(x) be of limited variation for y,(h) =asSrsb 
=y2(h), continuous at a and b, and having regular discontinuities for 
a<a<b. Let Gi(y), Go(y) be of limited variation on h Sy Sq, con- 
tinuous at h, and continuous from the left forh<ySqg. Let x, and x2 have 
continuous second derivatives. Then in the class of functions (14), there 
exists a unique function u(ar,y) satisfying the heat equation (1) in D and 
satisfying the boundary conditions (2) and (3). 


4. A note on paper (A). In our first paper the function /(z, y) solving 
the boundary value problem was written 
i=1 a 


We wish to point out that by the use of Theorem 4 of this paper we could 
drop out the xi(7)U and write 


(19) ia) + vaB(e). 


i=1 h 


This simplifies 1(x, y). The kernels of the integral equations that (19) leads 
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to are the same as the functions Kj; of (17). These kernels are not simpler 
in form than the kernels of paper (A). However, the 7; then become those 
defined in (18) and are simpler than those appearing in our first paper. 


PART 2 


In this part of the paper we return to the classical problem of finding a 
solution of the heat equation (1) in the region D which takes on preassigned 
values on the boundary of D. Here, however, our boundary values are not 
continuous but are functions of limited variation. The approaches to the 
boundary are in a restricted sense that we describe in the next section. The 
Riemann-Stieltjes integral will suffice for the work of this part of the paper, 

5. Definitions and notation. A point P(z,y) lying in the region D 
will be said to approach a point M(a ,h) on the characteristic y = in the 
parabolic sense if 


(20) lim — 2%) /(y—h)* = 0. 


When a point P approaches a point Q(yi(Yo), Yo) through a set of points 
E of D, we divide F into two sets as follows: F, contains those points of F 
such that y > yo, and LF, = H—E#,. If Q is not a limit point of F,, or if 
Q is a limit point of F, and in addition 
(21) tim | — xs(yo)|/(y — yo)! = 


P 


we shall say that P approaches the point Q in the parabolic sense.'° 


Let W(z,y) denote the function 


(22) W(x, y) = f _exp(—u*)du, y > 0; =(, y= 0. 


y 
The following functions defined in terms of W are solutions of the heat 
equation (1) for points in D, and will play a fundamental réle in the rest 


of the paper: 


ah 
(23) u(a,y) [4—W(a—éy—h) 
(24) =2 f W((—1)**(a@ — xi(n)), 0) dG(n), eA xi ly); 
h 
h t] 


19 The above definitions of parabolic approaches are radically different, and it is 
to be noted that a corner point (x;(%),4) can not be reached simultaneously under a) 


both definitions. 
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The functions F(é) and G(n) are functions of bounded variation and 
a= xi(h), b= 2(h). 


6. Parabolic approaches. In this section we shall be concerned with 
the limits of w and v; under parabolic approaches to the boundary of D. 


THEOREM 6. If P(2,y) approaches the point M(a,h) (aSaSb) 
in the parabolic sense, then the function u(P) =u(a,y) defined by (23) 
has the following limit 

0) + F(a —0) F(b) + F(a) 


2 


lim u(P) 


where F(b+ 0) =F (b) and F(a—0) = F(a). 
Suppose first that a < a <b and write u(P) in the form 


Io db 
u(P) = f [4-— W(a—éy—h) + [4 —W(e@—éy—h) ]dF(€) 
=u,(P)+ u(P). 
Integrating u,(P?) by parts and replacing F'(é)dé by d f Feat, we may 
a 


write 
ui(P)=[3—W(a éy—h) (8) | + ( U(x —E,y—h)d& F(t)dt. 
a Ja a 
Under a parabolic approach the integrated part has the limit — F(a) /2. 


Since dt has a left-sided derivative at equal to F(x) 0), Theorem 
1 of paper -(A) shows that the above integral in u,(P) has the limit 


F(r,—0)/2. The limit of u,(P) is thus seen to be 


lim u,(P) = (a) 
P+M 2 2 


Noting that : F(t)dt has a right-sided derivative at x) equal to 


a 
F(z+ 0), the same type of argument shows that 


P(t) +0) F(d) 
5) 


lim us(P) = 


P+M 2 


~ 


From the above limits the theorem readily follows for a<a<b. If 2 is 
aor b then u,(P) or u.(P) is zero and no new considerations are needed in 
these cases, 

Turning now to the functions v;(P), we see that as the point P 
approaches the point M(2,h) with a < a < b these functions each have the 
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limit zero. Our next theorem is concerned with the limit of v;(P) as P 
approaches the point Q(yi(Yo), Yo). 


THEOREM 7. If P approaches Q in the parabolic sense, then the fune- 
tions vi(P) defined by (24) have the limit ij 


lim [4(Yo— 0) — G(h)] + (Yo), ¥0)- 


Let us consider the case of v,(P), and to simplify the notation write 
v(P) =v,(P) and x(7) =xi(y). It is then to be noted that >0 
for P in the region D. 

For a given §>0 we can select « >0 so that the total variation of 
G(n—0) + 0)] in y—eS gS yo is less than 
With « fixed, since y— yo, we assume y—e/2 SySyo+e/2, and write 
v(P) in the following form 


v(P) = 2W (2 — x(n), —7)dG(y—0) + 2WdG(y — 0) 


=I1,(z,y) + 12(2, y). 


Since G(h—0) = G(h) and the integrand W is continuous and equal to 


zero for y =7, the replacement of G() by G(y—0) is permissible. 
In J, the difference y— y is bounded away from zero, hence 
I, = Li(x(Yyo), yo). 
For y = yo, we have 
| y)| < 28, 


If y > yo, we write I.(z,y) as the sum of two integrals 


The first integral is seen to be less than 28, while the second can be put in 


the form 


y 

W(x— x(n), y—n)dG(yn + 90) + [E(yo + 0) — G(yo—0)] — 9- 
Here again the first term is less than 28. while the second has the limit zero 
when P approaches @ in the parabolic sense. Thus we see that 


lim | 7, | < 48. 


Since § is any positive number, and since 


= 


SOLUTIONS FOR THE HEAT EQUATION 


as P | 
(x(yo) — x(n), Yo— 0) 4G | < 28, 


we have 


‘une- 


1)dG(y—0) + — x(n)» Yo— 1) — W(0, 1) JdG(q — 0) 


write G(y,— 0) —G(h) + v(x(yo), Yo)- 


0 
. Theorem 7 is thus seen to hold for v(P) =v,(P) and a similar proof holds 
for v.(P). 
in 6, 7. The first boundary value problem. We shall prove the existence of 


vrite a function 1(z,y) which is a solution of (1) in D, and whose limit as the 
point P(x,y) approaches the point M(ao,h) in the parabolic sense is equal 
to a given function B(x.) of limited variation with regular discontinuities 
[here a= y,(h) < < xo(h) also the limit of 1(P) as P approaches 
in the parabolic sense a point Q(xi(Yo), Yo) With yo >h is equal to a given 
function H;(y.) of limited variation continuous from the left. 
1 to _ Let F(x) = B(x) fora<a2<b and let F(a) = F(b) =0, and set up 
the function 
(5) U(a,y)—S2 — — 0) dG (n) 


i=1 7h 


+ f 
a 
+ vo(z,y) + y). 
We see that in the domain D the function /(z,y) satisfies equation (1) 
if the Gi(m) are functions of limited variation. Since the v;(z,y) approach 


zero as P approaches M(2o,h) for a< x) <b, we see from Theorem 6 that 
under a parabolic approach 


(26) lim 1l(z, y) = B(x). 

»)s9- As P approaches Q(xi(Yo); Yo) in the parabolic sense with yo > h, we have 


from Theorem 7? 


zero (27) lim I(x, y) = [Gi(yo—0)— Gi(h)] + (xi (yo), Yo) + U(xi (Yo), 
j=l 


and we must show the existence of functions of limited variation Gi(y) so 
that the above limit is equal to the given functions: 


(28) Hi(Yyo). 
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Let us note that the vi(x;(y),y) approach zero with y approaching h 
when 1 = j and also when i ~j. Also as (xi(y) — xi(h))/(y — h)i 
approaches zero with y—h, we see from Theorem 6 that lim u(yi(y),y) a 
y approaches h exists and has a finite value. Taking this value for u(xi(h), h) 
the function u(yi(y),¥) is a continuous function of yin hSyS 4g. 

Set 
(29) fi(y)= Gi(y —0)— Gi(h) + u(yi(h), h)— Hi(h), y>h; =0, y=h, 
and on equating (28) to the right member of (27) with y= y, we have the 


following system of integral equations 
(30) f(y) — +E (i=i,2), 
j=l 
where 
gi(y) = [Hi(y) — Hi(h)] —Lu(xi(y), y) — u(xi(h), I, 
(31) Kis(ys9) — ], (0, 1)]. 


Here 8;' is 1 or 0 according asi—j or tj. The kernels Kj; can be put 


into the two integral forms 


—2x'i(t) | U (xi (t)— xi (9), 


We have seen ** that the system of integral equations (30) having kernels 
with the representations (32)-(33) and with gi(h) = 0, have unique solutions 
fi(y) of limited variation and continuous from the left provided that the 
gi(y) are such functions. The Hj(y) have these properties and _ the 
u(xi(y),y) are continuous; it remains to show that the u(yi(y),y) are of 
limited variation. We shall do this in the last section. Thus we state the 


(33) 


result: 


LEMMA 2. For a given set of functions gi(y) of limited variation, con- 
tinuous from the left, and vanishing at h, there exists a unique set of func 
tions fi(y) satisfying (30) and having these same properties. 


The equations (29) determine each Gj(y—0) to within an additive 
constant which does not affect the function I(x, y); moreover, /(z,y) is in- 


11 The constant u (x;(4),h) —H,(h) is put in to make 9; (hk) =0. 
12 Paper (A), pp. 651-653. 


420 
_ 


t—1)i 


rnels 
tions 
> the 

the 
of 
the 


SOLUTIONS FOR THE HEAT EQUATION 421 


variant under the replacement of Gi(y) by Gi(y—0). Considering then the 
gi(y), we see that if u(z,y) is fixed the gi(y) are unique. If we put 
on the restriction 


(34) F(a) = F(b) 


in forming the function u(z, y) defined in equation (23), it then follows from 
Theorem 6 that the only u(z,y) which, under a parabolic approach to 
M(a,h) (a < & < b), has the limiting value zero, is u(x, y) identically zero. 
We shall say that w(2,y) is regular if it is formed under the restriction (34). 


We now state the main theorem of Part 2: 


THEOREM 8. In the class of functions (25) with u(a,y) regular, there 
exists a unique function l(x,y) such that I(x, y) ts a solution of (1) in D 
and, under parabolic approaches to the boundary, takes on preassigned func- 
tions of limited variation. On the characteristic y=h the assigned function 
has regular discontinuities, elsewhere the assigned functions are continuous 


from the left. |Kach y’i(y) ts assumed continuous. | 


8. Functions of limited variation. In this section we shall show that 
the functions w(yi(y),y) appearing in the integral equations (30) are of 
bounded variation. 

Let x(y) have a continuous derivative and define 


*h X(y) 
(35) H(y) =  U(t—é,y—h)dtdF (é), h<y=qy, 


b X(y) 


yoh 


H(y) is seen to be equal to the function u(y(y).y) and Theorem 1 insures 


the existence of the above limit. 

LemMA 3. The function H(y) defined by (35) ts of limited variation 

We give the proof for x(i) =; the proof for other values of x(h) can 


be handled in a similar manner. 
Write H(y) in the form 


H(y) = U dt dF(é) +f j vaare, y>h 
Ja vé a b 


=T1,(y) + 12(y). 


ng h 
y) as 
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Under the transformation (¢— €)/2(y—-h)4 the function 


is seen to be monotonically decreasing with y for each €, and to have a total 
variation less than one on the set y > h. It then follows that the variation of 
I,(y) on this set is less than the total variation of F(é). 

With ¥.=y, >h, write 


b X (ye) 
T2(y2) U(t —&é,y2—h)dt dF(é) 


X 


exp | at dF (€) 
+ J, [er (— —exp( | aa 


Remembering that 6 = y(h) and that y’(7) is bounded, the first two integrals 
on the right are easily shown to be less than VW y2—h— Vyi—h times a 
constant C. The constant C depends only on the maximum of ,’(7) and the 
total variation of F(é). We obtain for the third integral a similar bound, 
if, after using the theorem of the mean, the following property of exp(—a) 
is used. If P,k are two positive constants with 0 < k < 1 and ¢ is a variable 
= 0, then there exists a constant K such that 


a? exp(— a) < K exp(— ka). 
Since the function Yy—h is of limited variation, the function /2(y) 
is seen to have a variation on h << yg less than a constant times the 


variation of Vy —h. 
As the function H(y) is bounded it then follows that it is of limited 


variation on the closed interval h = y= gq. 
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ANALYTIC CONTINUATION OF FUNCTIONS DEFINED BY 
FACTORIAL SERIES.* 


By GLen T. MILLER and Howarp K. HuGHEs. 


In this paper we shall be concerned with the determination of an analytic 

continuation of a function Q(z) defined by a factorial series of the form: 

oo m ) (— 1)m (e+) 

where & is a positive integer, and z is a complex variable. 

The results here obtained may be considered to be an extension of a 
theorem obtained by Hughes? for a function defined by a series of the type 
(1) when & =0. The methods used bear a close analogy to those used by 
Ford * and Newsom * in their work on functions defined by Maclaurin series. 

It is well known that the region of convergence of a factorial series is a 
half-plane bounded on the left by a vertical line R(z) =A, where R(z) denotes 
the real part of z and —o ZA. Furthermore the convergence is 
uniform in the half-plane R(z) —A-+ 8 where 6 is any positive number. 
Consequently the series (1) represents an analytic function in the region 
R(z) >A.4 


The results are given in the following 


THEOREM. (tiven the factorial series (1), whose abscissa of convergence 
is finite, and for which k is a positive integer; let it be assumed that the 
function g(m) occurring in the coefficient of (1), when considered as a func- 
tion g(w) of the complex variable w=ax-+ iy, satisfies the following 
conditions : 


(a) it is single-valued and analytic throughout the half-plane x= — 3; 


* Received January 15, 1942. 

1H. K. Hughes, “On the analytical extension of functions defined by factorial 
series,” American Journal of Mathematics, vol. 53 (1931), pp. 771-776. 

°W. B. Ford, “The asymptotic developments of functions defined by Maclaurin 
series,” University of Michigan Studies, Scientific Series, Voi. 11. 

°C. V. Newsom, “ On the character of certain entire functions in distant portions 
of the plane,’ American Journal of Mathematics, vol. 60 (1938), pp. 561-572. 

‘For an account of the properties of factorial series, see The Calculus of Finite 
Differences by Milne-Thomson (1933), Ch. X. 
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(b) for all c= —}, and for all values of | y | sufficiently large, we have: 


+ ty) 
(?) g(x) 


where « is an arbitrarily small positive number and M is a constant de- 


< Mexpl{r(k + 4) 


pendent on e. 
Moreover, let there exist a positive constant q, dependent on e, but inde- 


pendent of z, such that when z is confined to any finite region, we have, 
g(x) | N 
(3) < 
T(iz+2+1) | 
when x > q, and where N is a positive constant. 
Then the function Q(z), defined by (1) may be continued analytically 


throughout the whole finite plane, except in the neighborhoods of the points 
z= 0,—1,—2,---, and throughout this region will be defined by the 


gite 


equation 
© sin krax g(x) 
+00 F(y, z) 
+ T(z) (cosh ry )** 
where 


[sinh — y) ]*19(— 4 }+ is) 
(5) (9, 0 r(z+4-+1s) 


The proof of this theorem is based on the following theorem in the 


ls. 


Calculus of residues, due to Newsom: 

“Let P(w) be regular within a closed region bounded by the curve 
C in the w-complex plane, and enclosing the points 1, °°, 12, 
but no other vanishing points of sinww. Furthermore, let P(m) 40, 
Then 

P(w 1 m(k+1) k 

where m= 0,1,2,: --,n, and where (3a,%.° a) denotes the sum of the 
({;) products of ¢ factors each formed by taking the possible combinations 
of the k quantities + (k —1)i, + (k—3)1,-- - 0 ( t at a time. Here 
rs is +7or0 according as k is even or odd. The integration is performe 
in the positive sense.” ° 


5 Newsom, op. cit., pp. 562, 563. 
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Proof of the Theorem. 


Let us suppose that z has any fixed value which is real and positive and 
greater than A, the abscissa of convergence of series (1). We shall write the 
series (1) in the form: 

g(m)E(2)(— 1)" 


and shall seek to determine P(w) so that 


ot (S (k-t) g(m) 


This is a differential equation in m with constant coefficients. The 


roots of the characteristic equation are distinct, and are + (k—1)m, 


+ (k—3)mi,--- oh . lf we denote these roots by rjzt, where 
(9) rj=k+1— 2, is 


then a particular solution of (8) is given by 


Pp = >) exp rjz d 
(10) (m) exp rjrim f, m 


where 
1 1 
(11) bj= II 


(2m1)*? az=1 


In some of our later considerations we shall find the following form for 


b; useful : 


1 k—1 


If we replace m by the complex variable w= az -+ iy, then P(w) is an 
analytic function of w in the finite w-plane. Hence if P(w) has the form 
indicated by (10), then the equation (6) can be written in the form: 


T'(z)2*k! P(w) 
etm+i1) Jc (sin rw) 
where the contour C is so chosen as to enclose the points w = 0,1, 2,- - -,n, 


but no other vanishing points of sin w. 

If we substitute the form for P(w) from (10) in (13) we obtain a 
contour integral for the determination of the first n + 1 terms of the series 
(7) as follows: 


~ 


9) 
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)m(k+1) 9m) § exp rjyriw I(w, 2) 


(—1 
> T(z+m-+1) 2% bj (sin ww )** dw 
where 
(15) I(w,z) = dw, 


and where b; is given by (11). 

We shall choose as the path of integration the rectangle enclosed by 
the lines (1) w=a2-+ ip, (2) w=—}+ (3) w — ip, and 
(4) w=n+4-+4y, and shall denote the integrals along these sides by 
A, B, D, and EF, respectively. Then in order to evaluate the contour integral 


(sin ww )** 


we must consider & integrals of the form: 


(16) f exp I(w, 2) 
JC 


(sin ww 


where 1 =j =k, and I(w,z) is given by (15). 

In (15) let us take the path of integration along the z-axis from w =0 
to w =, and then along a line parallel to the y-axis to the point w = x + ty. 
The function 7(w,z) then can be written in the form: 


I(w,z) = R(x, z) + 18(2, y, 2), 
where 


(7 exp(—rjriz) g(z) 
R(2z,2) = f, i dx, 


and 


S(2. y,2 — f° ty) dy. 
(9.9) +1 +9) 
With regard to R(z,z) it is sufficient for our purpose to note that 
| R(x, z)| is bounded when z is bounded. 
It is known that for all z > 2, where 2, is arbitrary, and for all y 
sufficiently large, we have: ® 


‘ 


(17) | iy) 
From the above inequality, together with (2) and (9), it follows that 
the absolute value of the integrand in (15) is less than 


m, exp[ryry + {(k + 1) —te}ly J, 


< 


® Ford, op. cit., p. 61. 
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the constant m, being dependent only on e«. It follows that there exists a 
constant K. dependent only upon e, such that 


(18) | S(x,y, 2)j < Kexp[rjyry + +1) — y |]. 
We are now ready to consider in detail the various contributions to the 
integral (16). We shall first consider the contribution A. Here we have 
y= ip, and 
sin rw = sin rp(sin rz coth rp + i cos rz). 


Since lim (sin rx cothp + t cos rz) = 1, sin rw is of the order of exp rp 
px 


when p is large. 
In order to evaluate A we must consider & integrals of the form: 


rjri(x + ip)[R(a,z) + iS (2, p,z)] 


J exp mp(k + 1) 
where 1 =j=k. It is evident that we may write 
R(2,z “4 p, 2) 
(9) 1415 | exp mp(k-+ 1+ +| na exp 


Since R(z,z) is bounded if z is bounded, the first integral on the right 
of (19) approaches zero as p approaches infinity. 
With regard to the second integral on the right of (19), we have from 
(18) that 
| S(x, p,z)| K + & +1) — Help). 


The absolute value of the integrand in the integral in question is there- 
fore of the order of exp(— dep) for p large, and since « is positive we conclude 
that lim A; = 0. It follows that lim A = 0. 


In a similar manner it may be shown that lim D = 0. 
pw 


We next consider the contribution #. We shall choose n to be an even 
integer so that sin rw cosh zy. F consists of k integrals of the form: 
4 
exp (— rjr 
irs iR(n + 4,2) ay 
(cosh ry)* 


exp (— rjry)S(n + z) ay | e 


(cosh wy )** 


(20) 


We shall set 


*© exp(— rjry) 
f. (cosh wry )** 


Then by using the result of Newsom’ and (9) we can write 


"Newsom, op. cit., pp. 567, 568. 


= 
\ 
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ye jQk k 


Let us next consider Aa second integral on the right of (20). The 
integral can be written: 


*© exp(— ryry)S(n + 4,y,2) 
(31) (cosh wy) **? dy 
(2z+n+1) (cosh wry )**? “y 
where 


ty) 


We now make use of the following well known property of the gamma 


function : 
T(x + ty) = (27)4(x exp(— — yarc tan y/r) (1+ 8) 

where 6 approaches zero as either x or y becomes infinite. Then we have 
T(iz+n+1) | 
| 
By the use of (22) and (2), and taking into account the behavior of 
(cosh ry) for | y| large, it follows that the integrand of the right member 
of (21) is of the order of exp(—e|y|), for | y| large, and hence the integral 
converges. Furthermore, since the series (1) was assumed convergent, we 


(22) < K exp (4/2 | yl). 


have: 


g(n) 
lim 
Consequently the integral on the left of (21) tends to zero as n approaches 
infinity. The total contribution H#; is given by 


lim R(n + $,2 


OO 


We next consider the contribution B arising from the integration along 
the side where iy. Here sinww = cosh zy and we must con- 
sider & integrals of the form: 


(— cosh zy) *# 


exp rjri(— + ty) 3, 2) 
dy. 
(— cosh zy) **? 


co 


® Ford, op. cit., p. 61. 
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The first of the integrals above, when combined with the contribution E,, 
can be written in the form: 


iH; lim [R(n + 4,2) —R(—4, z)] 


exp (— rjriz) 


If we multiply the above integral by 6;, we obtain as the sum of the 


j integrals in question when substituted in (14) 


© sin g(x 
T(z+2+1) 
Finally let us investigate the second integral in (23). When substituted 
in (14), the sum of the & integrals in question may be written in the form: 


k 


iy P@+$+is) 
(cosh )**? dy. 
But 
(26) (— 1)?! exp rjr(s — y) = (sinh r(s — y) )*. 
ja 
Upon substituting (26) in (25), we find the value of this contribution 
to be: 


r(z+4-+ is) 


(27) 
(cosh wy) ** ay 


[sinh —y) 3 +18) 
0 


The value of the right member of (14) is therefore the sum of the 
expressions (24) and (27). We thus obtain equation (4). 

Up to this point we have restricted z to values which are real, positive 
and greater than A, the abscissa of convergence. ©(z) as defined by series 
(1) is a function which is analytic throughout the half-plane R(z) >A with 
the exception of the points z= 0,—1,—2,-- +. The right member of (4) 
furnishes an analytic continuation of Q(z) in any finite region which does not 
contain the points z= 0,—1,—2,- - -, since throughout such a region the 
member under discussion is an analytic function. To prove this last state- 
ment, it is sufficient to show that the integrals on the right of (4) converge 
uniformly for all values of z in the proposed region. From (3), we see that 
the first integral in the right member of (4) converges uniformly, since it is 
dominated by ("x «dr, From (2), (5), and (17), it is evident that the 


oo 
second integral converges uniformly, since it is dominated by f exp (— gey) dy. 


The theorem is therefore proved. 


| | — | 4 is) 
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Extensions of the Theorem. 


In the hypotheses of the theorem it was required that the function g(w) 
be analytic for all values of w=2x-+ iy for which r= — 4. This condition 
is unnecessarily restrictive, as we shall now show. Let us assume that g(w) 
has a singularity at a point ww, in this region, and that wy, is not an 
integer. We shall now show that. the theorem continues to hold provided 
that we subtract from the right member of (4) the loop integral 


T(z) T(z +w+1) 
g(w) 


Ab 


dw 


wherein the contour LZ, is described as follows: Draw a straight line from 
w = wy, to infinity in a direction of the fourth quadrant, and let this line be 
regarded as a cut in the w-plane. Let the contour ZL, consist of the two lines 
drawn on either side of and arbitrarily close to the cut, the ends near the 
point w= w, being joined by a circular arc of small radius about the same 
point. The integration is to be performed in the positive sense. We note that 
J, is analytic for all finite values of z except z=0,—1,—2,---. 

An examination of the analysis of the proof of the theorem shows that 
if g(w) has the singularity mentioned, then the theorem still holds if we 
subtract from the right member of (4) the loop integral 


L 


24 (Sin rw )**? 


where L, denotes the contour just described. We shall show that H, can be 
reduced to Jp. 

Noting the value of P(w) given by (10), it is evident that we may write 
P(w) in the following form: 


P(w) = ‘ exp rjri(w —t) 


By using the form of 6; given by (12) it is easily verified that 


[sin — t) 


m1(k—1)! 


> b; exp rjri(w —1) = 


Thus we have 


: g(t) 


— dw. 


Lp (sin aw 
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Let w be a point on the cut, and let L’, be that part of Z, drawn from 
w, to w, on both sides of the cut. Let 


(w) 


t) 
(w) (sin rw) dw. 
an 
del If we use the multiple angle formula and expand [sin x(w—t) ]** by 
the binomial theorem, we obtain: 
| (2) k—1 f (cos rt sin nt)s-19(t) 
i) J (— vot nw) esc? ru f, dt dw 
rom & After integrating by parts and making simple reductions, we obtain 
ines. 2 [cot nw sin rw (cos rw)*-i — (cos rw)* 
the sin + w+1) 
ame 
‘hat If we carry out the summation under the integral sign and make simple 
transformations, we obtain 
sin7w T(z+w-+1) dw. 
On setting 
g(w) 
— Lp Xp kriw T(z+w-+1) 
it is evident that lim J,(w) =Jd>. 
be 
te I We shall show that as # > © along the cut, we have lim [J,(w) — H,(w) ] 
Tite 
=0, and it will follow that lim H,(w#) =J»,. Now we have: 
w 
J,(w) — H,(w) = T(z) exp (— mw) 
) exp(— (k — m) exp(— (m—1) g(w) 
x sing7w T(z+w+1) dw. 


But if # —> o in the fourth quadrant we have for w large and for all w on L’y 


sin rw 


<1 and exp(— iw) < exp(— iw) and hence 


sin 


| g(w) 
(w) — ») | exp(— ri ok 


| 
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The integral exists because of the hypotheses of the theorem, and since 
the factor exp(— iw) approaches zero, the above quantity can be made arbi- 
trarily small. Hence if the function g(w) has a finite number of singularities 
at the wp, p= 1,2,- - -,n, the theorem continues to hold if we subtract from 
the right member of (4) the sum of the n integrals Jp. We note that if the 
singularity of g(w) at W, is a pole, then the integral J, reduces to the residue 
of the function 

P(w) 


P(2)*"sk! —— 
(sin ww )** 


where w= w, and P(w) is given by (10). 
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A KIND OF EVEN-DIMENSIONAL DIFFERENTIAL GEOMETRY 
AND ITS APPLICATION TO EXTERIOR CALCULUS.* 


By Hwa-Cuune LEE.’ 


Introduction. In this paper, we study the differential geometry of an 
even-dimensional space, into which a nonsingular skewsymmetric matrix is 
introduced as fundamental tensor. In this geometry we define a curvature 
tensor from which there arises the distinction between general, flat and con- 
formally flat spaces. Necessary and sufficient conditions for such spaces to 
be flat or conformally flat are found in 1, 2. The results are applied, in 3, 
to the exterior calculus to discover a necessary and sufficient condition that a 
nonsingular exterior differential form of degree two admit an integrating 
factor. The paper is concluded in 4 with a brief study of the infinitesimal 


automorphisms of the spaces. 


1. The space L.,. Consider a manifold of even dimensionality 2n 
whose points are described by a system of codrdinates z* (a, B, y, p, o, T 
=1,---+,2n), and let agg be a nonsingular skewsymmetric covariant tensor 
whose components are analytic functions of the z’s. The manifold equipped 
with such a tensor is denoted by Lgn, and the tensor is called its fundamental 
tensor. It is at once verified that the skewsymmetric quantity 


1.1) Kapy = 0,a Ogitya + = 0/0x% 
By By y“aB 


is a covariant tensor, which we call the curvature tensor of Lon. This is the 
simplest tensor involving the derivatives of the fundamental tensor. The 


curvature tensor evidently satisfies the identity 


(1. 2) Op K apy ppy OpK apy K agp 0 


which is a tensorial equation, holding for all codrdinate systems. 
In terms of Cartan’s well known notation in the exterior calculus,? there 
is associated with dag the erterior differential form of degree two 


(1.3) Q = 


which we call the fundamental form of Lon. As dag is nonsingular, Q is also 
called nonsingular. The exterior derivative of Q is 


* Received October 27, 1941. 
*The author wishes to express his thanks to Professor J. M. Thomas for reading 
the manuscript and giving valuable criticisms. 


* See e.g. E. Cartan, Legons sur les invariants intégraua, chap. VII. 
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(1. 4) = 


and (1.2) is equivalent to the well known fact 0” = 0. 
“An Ln is called flat if its curvature tensor vanishes. Evidently, every 
L, is flat. For a flat Le, we have 2’ —0 so that Q itself is the exterior 
derivative of a pfaffian form.* Taking, by a change of codrdinates, this 
pfaffian in its well known canonical form* we find that Q may be written as 
+ 


the last term being certainly [dr"dx"*'] because of the nonsingularity of Q. 
Hence in the new coordinate system the fundamental tensor assumes the form 


0 J 
== 
(1. 5) Aap 2 - 


where J denotes the matrix of order n with 1’s on the secondary diagonal 
and 0’s elsewhere and 0 denotes the zero matrix of order n. 
It is also clear that there are other codrdinate systems for which the com- 


ponents of dag are any constants. The converse is evident, and we have 


THEOREM 1. An Do» is flat if und only if there exists a coordinate system 
for which the components of the fundamental tensor are constants. 


2. Conformal spaces. An J», and an “Lo, over the same underlying 
co] 

manifold are said to be conformal to each other if their respective fundamental 

tensors dag(z) and ‘dag(x), referred to the same codrdinate system, are con- 

nected by a relation of the form 

(2. 1) ‘dap = dilap, 

where ¢ is a nonvanishing scalar function of the codrdinates in question. 

Evidently, every L. is conformal to every other Lz. We shall assume 2n >? 

throughout this section. By (2.1) we find 

4 

(2. 2) K apy = ¢(Kapy + log + log + dapi, log 

Since dag is nonsingular, it has an inverse a%* (skewsymmetric in @, 8) 
defined by a*¥a,g = 8g* where 8* is 1 or 0 according as @ is or is not equal 
to B. Define the covariant vector 

Ka = 
and call it the curvature vector of Lon. Since by (2.1) we have ’a% = $"a®, 
we find from (2.2) 
(2.3) ’Ka = Ky, — (2n-— 2) log ¢. 


* See e.g. E. Goursat, Lecons sur le probleme de Pfaff, p. 105. 


“Goursat, loc. ctt., p. 42. 
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Consequently the skewsymmetric covariant tensor 
(2. 4) bap = g — OpKa 
is an obsolute conformal invariant : 
(2. 5) ‘bap bap. 
We call bag the first conformal curvature tensor of Len. 
The result of eliminating the derivatives of log¢@ in (2.2) by means 
of (2.3) may be written 


(2. 6) ‘Capy = > Capy 
where 
(2.7) Capy = Kapy + (Kadpy + Kpaya + — 2). 


Thus the skewsymmetric covariant tensor Cag, is a relative conformal invariant 
which we call the second conformal curvature tensor of Lon. 

Solving (2.7) for Kag, and substituting into (1.2) give 
(2.8) dpabpy + ppbya + Apybap + bpattpy + bppaya + bpy4ap = Cpapy: 
where for the sake of brevity Cpagy indicates the aggregate of terms independent 
of the derivatives of K,. Because of (1.1) and (2.7) this expression may 
be written 


Cpapy = (2n — 2) (OpCapy — DaCppy — IpCpya — 9Cpag) 

+ (Kpcapy — Kappy — K pepya — KyCpag) 
(which is a skewsymmetric covariant tensor). Transvecting (2.8) by a°7 and 
writing Cag = 4" Cpagy we find 


(2n — 4) bag = baag — Cag 


where b ab. Transvecting the above equation by a% and writing 


c= we obtain 2(2n —2)b = —-c so that 
(2.9) — 2(2n — 2) (2n — 4) bag = cdag + 2(2n — 2) cag. 
We observe that (2.4) and (2.7) immediately imply 
THEOREM 2. The first conformal curvature tensor vanishes if, and only 


if, the curvature vector is a gradient and the second conformal curvature 


tensor vanishes if, and only if, the curvature tensor is of the form 
9 

(*. 10) apy = Jalpy Jptya + Jylap 

where Ja is a covariant vector. 


When 2n—4, (2.10) consists of four equations which are directly 
verified to be solvable for the four unknowns 4);, 92, 9s; ga. Hence 
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THEOREM 3. The second conformal curvature tensor of every L4 vanishes 
identically. 

An Lon is said to be conformally flat if a scalar ¢ exists such that the 
’Lon defined by (2.1) is flat. Suppose that Len is conformally flat. If ’L,, 
is the corresponding flat space, then ’Kag,—0 and therefore ’Kq = 0, ‘bag =0, 
‘Capy = 0 so that by (2.5) and (2.6) we have bag = 0, Cagy = 0. Suppose 
conversely that the latter hold. Then Kg is a gradient (by the first part of 
Theorem 2) and ’cagy—0 (by (2.6)). The former fact shows, by (2.3), 
that we can make ’K, vanish by properly choosing ¢. From ’cagy = 0 and 
’K, =0 we have ’Kagy =0 by (2.7), whence ‘Zon is flat and therefore Lo, 
is conformally flat. Thus we have proved 

THEOREM 4. A necessary and sufficient condition that an Lon (n >1) 
be conformally flat is that the first and second conformal curvature tensors 
vanish. 

When 2n > 4, the vanishing of the first conformal curvature tensor is a 
consequence of the vanishing of the second conformal curvature tensor, as is 
at once seen from (2.9). Hence, in view of Theorem 3 for 2n = 4, Theorem 4 


may be refined into 


THEOREM 5. A necessary and sufficient condition that an Len be con- 
formally flat is that the first conformal curvature tensor vanish for 2n =4 
and the second conformal curvature tensor vanish for 2n > 4. 

Given a conformally flat Z.,, the determination of the corresponding flat 
space is now immediate. For in this case we have, by Theorems 4, 2 equation 
(2.10) where gg =— Ka/(2n —2). The vector gq is a gradient and so may 
be written in the form 
(2.11) Ja = — 4g log ¢, 


where @ is different from zero and is determined save for a constant factor. 
Because of (1.1) and (2.11), (2.10) may be written 

(2. 12) + + = 0 

which shows that ¢aag is the required fundamental tensor of the corresponding 
flat space. (2.12) also implies the existence of a covariant vector dg, de 
termined save for an additive gradient, such that 

(2. 13) blag = — Ogda. 


3. Integrating factors of non-singular exterior differential forms of 
degree two. The theorems of the preceding section incidentally solve @ 
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particular problem in the exterior calculus. An exterior differential form 
(of any degree) is called an exact differential if its exterior derivative vanishes. 
An exterior differential form is said to admit an integrating factor if it can 
be made an exact differential by multiplication of a scalar factor.° Except 
for the case when the degree of the form is one (pfaffian form), a necessary 
and sufficient condition for a form to admit an integrating factor has not yet 
been found. 

It is clear that the form © defined by (1.3) admits an integrating factor 
if and only if the corresponding Lz, 1s conformally flat. When 2n> 4, 
Theorems 5 and 2 imply that LZ, is conformally flat if and only if (2.10) 
holds. Write the latter in the form 


(3.1) = [oO], 
where w denotes the pfaffian ggdz*. Thus Q is a factor of 9’,® and we have 


THEOREM 6. <A necessary and sufficient condition that a nonsingular 
exterior differential form of degree two in an even number (> 4) of variables 
admit an integrating factor ts that the form itself be a factor of its exterior 
derivative. 


Again by Theorems 5 and 2, an L, is conformally flat if and only if Kg is 
a gradient, i.e. the pfaffian form 6 = K,dzr* is an exact differential. Now 


= K,dx* = Kap = aK = a 


where 04g = Kag dx? is a system of (six) pfaftian forms, said to be associated 
with the form ’.? Calling 6=a@ag the contraction® of this associated 


system, we have 


THEOREM 7. A necessary and sufficient condition that a nonsingular 
exterior differential form of degree two in four variables admit an integrating 
factor is that the pfaffian form, which is the contraction of the system of 
plaffian forms associated with the exterior derivative of the given form, be an 
exact differential. 


4. Infinitesimal automorphisms of L.,. We conclude this paper by 
giving a brief study of the infinitesimal automorphisms of an Len, which are, 
by definition, infinitesimal (point) transformations leaving invariant the 
fundamental form ©. In the same way as Killing’s equation in Riemannian 


® Goursat, loc. cit., p. 118. 
*Goursat, loc. cit., p. 97. 
*Goursat, loc. cit., p. 126. 
*This operation is unfamiliar in the exterior calculus. 
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geometry is obtained, we find that a contravariant vector v* is the vector of 
an infinitesimal automorphism of Ls, if and only if it satisfies the equation 


(4.1) + + = 0. 


Introducing the covariant associate vg = dagv® of v* we may write this equation 
in the form 

(4. 2) 0ag — Opva = — VpK? ag, 

where K? ag = @7K yap. 

We shall only try to solve (4.2) in particular cases. If Lon is flat, we 
have K’ag—0 so that (4.2) implies that a covariant vector is the vector 
of an infinitesimal automorphism of a flat Lo, if and only if it is a gradient. 
More generally, if Lo, is conformally flat we may, using (2.10) and (2. 11), 
write (4.2) in the form 
(4. 3) — = Vpdag, 
where v = v’gp, which can be shown to be a constant. In fact, if we operate 
on (4.3) with @,, perform a cyclic permutation of a, 8,y and then add the 
three equations thus obtained, we find, on account of (2.12), ag,0av + dyadgu 
+ = 0, whence and v=const. Consequently, in view of 
(2.13), the general solution of (4.3) is of the form 
(4. 4) = + 

If we introduce the contravariant associate g* agp of gq and transvect 
(4.4) by g*, we find 

(4. 5) = —v(¢+a) 

where @ = g%dg. Hence 


THEOREM 8. A covariant vector vq is the vector of an infinitesimal auto- 
morphism of a conformally flat Len if and only if it is given by (4.4) where 
v is a constant and w is a solution of the linear equation (4.5). 


In this conformally flat case (n >1) it is,easily seen that v = 0 if and 


only if v, satisfies 


Vq(Ogvy -— Oyvg) Ova — Davy) + Vy (avg — Ogva) = 0, 


i.@., Ug is proportional to a gradient d,p, and this is the case if and only if 


the trajectories of the 1-parameter group of automorphisms generated by 
are at the same time trajectories of other 1-parameter groups of auto 
morphisms, the generators of the latter being then necessarily v* multiplied 
by functions of p. 


NATIONAL SZECHWAN UNIVERSITY, OMEI, CHINA. 
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OSCULATING CONICS OF THE PLANE SECTIONS THROUGH 
A POINT OF A SURFACE.* 


By BucuHIn Sv. 


1. Introduction. According to the classical theorem of Moutard all the 
osculating conics of oo° plane sections through a generic point O of an analytic 
surface immersed in a three-dimensional projective space are distributed on oo 
Moutard quadrics each of which corresponds to a non-asymptotic tangent of 
the surface at O. In connection with the definition of Darboux pencils of 
quadrics we have studied * Moutard pencils of quadrics (so called because each 
of these pencils contains a Moutard quadric). This notion has been generalized 
ingeniously by Louis Green * who shows, among other things, that two arbitrary 
non-asymptotic tangents at O suffice to determine a certain pencil of quadrics 
having a contact of the second order with the surface at O. 

In studying the residual osculating conic on a given Moutard quadric 
Chenkuo Pa has recently proved that * if a quadric through the two asymptotic 
tangents at. O possesses one and only one osculating conic of the surface at the 
same point, then it must necessarily belong to a Moutard pencil. Combining 
the investigation of L. Green with that of Pa we are naturally led to inquire 
whether the osculating conics belonging to different Moutard quadrics at O 
may be regulated in such a manner that they lie on a quadric of certain pencils. 

The purpose of this paper is to solve the problem of determining a quadric 
by means of an osculating conic of a given direction and another osculating 
conic of another given direction with a suitably chosen plane, showing that on 
the quadric thus obtained there exists a third osculating conic. Thus all the 
osculating conics of the plane sections through O are newly distributed, three 
by three, on a quadric having a contact of the second order with the surface 
at O and meeting it in a curve containing these tangents as two of its triple 
tangents. 


2. Analytic basis. If the four homogeneous projective coordinates z of 
a variable point on an analytic non-ruled surface S in ordinary space are given 


* Received July 17, 1941. 

1Buchin Su and Ichida, “ On certain cones connected with a surface in the affine 
space,” Japanese Journal of Mathematics, vol. 10 (1933), pp. 209-216. 

* Louis Green, “ Systems of quadrics associated with a point of a surface,” American 
Journal of Mathematics, vol. 60 (1938), pp. 649-666. 

*C. Pa, “On the quadrics of Moutard,” not yet published. 
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as analytic functions of two independent variables u, v, and if the parametric 
net on S§ is the asymptotic net, then the functions z are solutions of a com- 
pletely integrable system of differential equations in Fubini’s canonical form 


Luu = + Bry + pz, 
Lov = yLu + Xv + Gz (6 = log By). 


The coefficients of these equations are functions of wu, v and satisfy certain 


(1) 


conditions of integrability. 

The codrdinates of a point O’ on S near O can be represented by the form 
LyX LoLy + + For the non-homogeneous coordinates y, z of 0’ 
defined by putting 


we have the canonical expansion in the form 


(3) ay — + yy") + y) +- 


where 
(4) 9) = — y — — + yyy") 


and 


__ dlog By’ dlog B’y , 
(5) 


Consider the section of the surface produced by a plane through a non- 


asymptotic tangent ¢, whose equations are 
(6) z= 0, y— nz =0; 
the equation of this plane is evidently of the form 


(7) z= p(y—nz), 


so that the osculating conic of the plane section at O is given by the equations 


(7) and 


(8) —2°-+ YL{(p/3n*) (8 + yn*) — (1/n) }e 
— (p?/n?) {(1/3np) (B — 2yn*) + (1/n) 2) 
— (1/9n2) yn®)*}¥ + (p/n)] =0, 


where we have placed 
(9) Y =y— nz. 

‘Cfr. E. P. Lane, “ Power series expansions in the neighborhood of a point on 4 
surface,” Proceedings of the National Academy of Sciences, vol. 13 (1927), pp. 808-313. 
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Elimination of p from (7) and (8) gives the equation of the Moutard 
quadric belonging to the given tangent ¢,, namely, 


(10) (z— ry) + 12n? (28 — yn*) xz — 12n(B —2yn*) yz 
+ {4(B + yn*)? — 36n¢g,4(1, n) }2? = 0. 


3. Quadrics having a contact of the second order with a surface. The 
equation of the most general quadric through the two asymptotic tangents of 
the surface S at the point O is 


(11) ry + kyz + + + = 0, 


where k,, kz, ky, ky are arbitrary functions of wu and v. The conic intersection 


of this quadric and the plane (7) is expressed by (7) and 


(12) Y[{— (p/n) (ko + kgn) —1/n}a 
— (p?/n) ((1/p) hs + ks) ¥ — ky (p/n) ] = 0. 


In order that this conic be the osculating conic (8) of the plane section of the 
surface produced by the plane (7), it is necessary and sufficient that k.—=—1 
and 
( ko + hyn = — (1/3n) (B+ yn*), 
(13) 4 kz — (1/8n?) (B— 2yn*) = p{ (1/n?) (1, 
—- (1/9n*) (B+ yn*)? — ky}. 
Suppose that two arbitrary non-asymptotic tangents to S at O are given, 


respectively, by the equations 
(14) 0), nix = 0 (1 = 1, 2) 


and that through each of them there passes a plane 
(15) z= pi(y— nix) (4 = 1,2) 


on which lies the corresponding osculating conic K; of the section at O. If 
the quadric (11) passes through both conics K, and K, the conditions (18) 
must hold for n = ni, p= pi (t= 1,2) and the same set of and ky; thus 


ky + kan, = — (1/3m1) (B + yn’), 
hy -+ kgng = — (1/82) (B + 


kts — (1/3n,?) (8B — 2yn,') 
(17) = pi{ bs(1, m1) — (1/91?) (B+ yni?)? — ka}, 

ks — (1/8n2) (B — 2ynz*) 

| pof (1/ms?) b4(1, m2) — (1/9n2*) (B+ yne®)? — ky}. 


6 
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Solving (16) for k. and k; we obtain 


key = — (1/3) {(B/ninz) (ry + ne) — yryno}, 
ks = (1/8) {B/nynz— y(m1 + ne) 

The plane through one of the tangents ¢, and ¢, may be assigned arbitrarily 
and that through the other tangent is perfectly determined, since (17) gives 
k, for a given p, and therefore determines p2, except when the quadric thus 
obtained coincides with a non-Moutard quadric of the Moutard pencil belonging 
to the tangent ¢t, or t.. In the exceptional case the plane (14) coincides with 
the tangent plane of the surface at the point O, so that no osculating conic 
can be obtained. That is to say, there is one and only one osculating conic on 
a non-Moutard quadric of any Moutard pencil, as was shown by Pa. 

Let us consider the general case, confining ourselves to the determination 


(18) 


of the other osculating conics on the quadric in consideration. Writing the 


equation of the plane of such a conic, if any, in the form 


(19) z= p3(¥ — 32), 
we have 
(20) — (1/3n;*) (B 2yn;") 
= ps{ (1/ns?) (1, m3) — (1/9n3) (B + — ky}. 


The last condition evidently determines p; and consequently the plane (19) 
provided that the quadric is not the Moutard quadric belonging to the tan- 
gent ¢;. Substituting the values of /. and k, given by (18) into the first 
equation of (20) and reducing, we are led to the relation 


(21) Ni N2Ns = — B/y, 


which shows that the third tangent ¢, is completely determinate when the first 
two ¢, and ¢. are known. 
By virtue of (21) we can express (18) in the symmetrical form 
(22) { = (myn. + nons + 
ks hy (ny + Ne + Ns). 


Thus we arrive at the following conclusion: 


Given two arbitrary non-asymptotic tangents t,, t. to S at a genene 
point O and one plane m, through t, we can always determine a quadric Q 
through the two asymptotic tangents of S at O such that the osculating con 
of the section of S produced by 7m, as well as that of the section produced by 


wl 
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the corresponding plane m2 through t, hes on Q. This quadric, having a contact 
of the second order with S at O, contains a third osculating conic at O and 
forms a pencil when the plane m, turns about the tangent t,. 


4. A covariant line associated with two non-asymptotic tangents. Let 
us consider now the correspondence between the planes z, and zz passing 
through the non-asymptotic tangents ¢, and f2, respectively, such that the 
osculating conics of the sections produced by 7, z2 and the asymptotic tangents 
lie on the same quadric. By virtue of (17) it is easily seen that the planes 
(15) for 11, 2 are in the correspondence given by the equation 


(28) (1/p:) {ks — (1/3n,2) (8 — 2yn,*)} 
— (1/p2) {ks — (1/3n2*) (B — 2yns*)} 
= (1/n,?)$4(1, m1) — (1/2?) (1, ne) 
—(1/9n,*) (B + + (1/9n2*) (B + 


and are, therefore, in perspective position. 


In order to derive the plane, denoted by 2,2, on which lie all the lines of 
intersection of the corresponding pairs of planes x, and 72, we have to substitute 


1/pi = (1/2) (y— nz) (1 = 1,2) 


into (23) and then to reduce the resulting equation by means of the relations 
(21) and (22). A simple calculation suffices to obtain the equation of the 
plane namely, 


(24) B(2nz— p.)x + ynz(p2— 38nynz)y 


n No 
Y 
+ — — + + no”) } 
Ny No 
where we have placed 
(26) =n, + + Ng, Po = + + 


If we take the three tangents t; (i = 1, 2,3) at O given by 
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where 7;, M2, Ng satisfy (21), then there are three planes m2, 223, 731; Which 
correspond to the pairs of tangents ¢,, t2; t2, t; and ts, ¢; respectively. These 
planes are evidently concurrent in a line / associated with the tangents ¢, and 
t, (and consequently ¢;) of the surface S at O. 

From the equation (24), and that obtained by interchanging n, and n, 
the equations of / follow immediately. The result of carrying out the compn- 


tation is as follows: 


(28) y:2—=4(— + 2ypob — — + pe 
15B¢ — — ypip2) — Byps? — 9B + ypy 


Several particular cases of the above result are worth noting. If the three 
tangents are those of Darboux, then p,—p.,—0O and, consequently, the 
corresponding line / is the canonical line of the first kind for which k = — 5/12, 
which we shall designate as the first axis of Bompiani.® 

If two of the tangents lie in distinct Segre directions, the third is seen to 
lie in that Darboux direction which is conjugate to the residual Segre direction, 
Every Darboux direction thus determines a covariant line /. Putting 


= (B/y)*o"", m2 = = (B/y) 


where w® = 1, w ~ 1 and r= 1, 2, 3, so that 
2 (B/y) By", = 2 (B/y) 
we obtain the covariant line I, 


(11/20)y— (4/15) B*y%o" + (1/5) (y/B) 
(11/20) + (1/5) (B/y) — (4/15) 1. 


The harmonic polar of the tangent plane z = 0 with respect to the trihedron 


formed by 1,, 1., 1; is easily found to be the canonical line of the first kind fof 


which k = — 11/20. 
If two of the tangents are conjugate, namely, such that their direction f 
are n, —n, the third must lie in the direction B/yn*, so that 
Pix B/yn*, 


and the covariant line 1, corresponding to any non-asymptotic direction n (ani 

its conjugate —n) is given by the equations 

(29) {— (7/16) (77/88) n* — (1/16) (y/B) on? — (1/24) 
:{— (7/16) — (B*/8y) (1/n*) — (1/16) (B/y)y(1/n?) — (1/24) 


5 Cf. L. Green, loc. cit.’, p. 654. 
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In particular, for the Darboux directions 


n= (B/y) (r= 1, 2,3) 


(and the Segre directions) we have the lines: 


wry: {— (7/16)y — (1/6) B*y*%o" — (1/16) (y/B) 
{— (7/16) — (1/16) (B/y) 40" — (1/6) : 1, 


where 7 = 1,2,3. In consequence, the harmonic polar of the tangent plane 
z= 0 with respect to the trihedron formed the these lines is the canonical line 
of the first kind for which k = — 7/16. 

If two of the tangents coincide in a non-asymptotic direction , the third 
must lie in the direction — B/yn*, so that 


pi = (1/yn*) (B— 2yn'), Po = — (1/yn) (2B — yn’). 


By means of these values there is no difficulty in determining the line (28) in 
question. Without expressing the equations explicitly we shall only consider 
two interesting cases, namely, the lines corresponding to a Segre direction and 
a Darboux direction. In the latter case we have 


n =— (B/y) 
so that 
= — 3(B/y) %o", P2 = — 3(B/y) 


and the line becomes 


{(1/3)¥— (1/12) (7/8) — y(B/y) 
(1 + (1/12) (B/y) 1. 


These form a trihedron with respect to which the harmonic polar of the tangent 
plane z = 0 is obviously the first axis of Cech. 

In the similar way we obtain a covariant line corresponding to every 
tangent of Segre; the harmonic polar in question is the canonical line of the 
first kind for which k = — 7/16. 


5. The quadrics of Darboux. Let us confine ourselves to an important 
case which can be derived from the general result stated in Section 3, namely, 
the case where the tangents are in Darboux directions. The quadric in question 
must be a Darboux quadric, since (22) gives k, =k; = 0 for the roots n,, no, 


ns of the equation 


B+ yn* = 0, 
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so that the equation of the quadric becomes 
(30) ry—z2+ kz? =0. 


The planes of osculating conics lying on this Darboux quadric are 


(31) z= pi(y— nr) (1 = 1, 2,3), 
where 

(32) 1/pi = — (1/8) [os (1, ni) — nirky] 

and 

(33) ni = — (B/y) = 1, 2,3) 


In order that these planes be concurrent in a line it is necessary and 
sufficient that 
| 1| 
No 1/pe2 
Ns 1/ps ] 
or 
(34) ks = — $6ur. 


Hence we have the following theorem: 


Every plane through one tangent of Darboux uniquely determines a 
Darboux quadric Q such that the three osculating conics of sections produced 
by this plane and two others each passing through the remaining tangents of 
Darboux lie on Q. The planes of these conics are concurrent if, and only if, 


Q is the quadric of Wilczynski-Bompiant. 


The common line of these planes, in the last case, can be found by sub- 


stituting (34), (32) into (31), so that the equation of each plane is 
nila + (5/12)u2] —-[y + (5/12) 2] = 0. 


Hence the planes of the three osculating conics on the quadric of Wilczynsk- 
Bompiani are concurrent in the first axis of Bompiant. 

Since the expression (32) of 1/p; for a given nj, is linear in hy, we may 
conclude that if 7, rz, mw be the tangent plane and the planes of osculating 
conics of sections through a non-asymptotic tangent and lying, respectively, 
on the quadric of Lie and that of Wilczynski-Bompiani and if 7 denotes the 


plane through the same tangent such that the double ratio 


) = h, 


and 
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the quadric determined by x ts the Darboux quadric 


(35) ry —2z— {40u + (1—h) By}z? = 0. 


Moreover, by 


mine the first and the second axes of Bompiani. In fact, the three planes of 


means of the pencil of Darboux quadrics, we can also deter- 


osculating conics lying on the quadric (30) form a trihedron of edges 


nix —y + (1/pi)z2=0, njx—y + (1/pj)z=0, 


where j, 1,7 = 1, 2,3 and p; is given by (32). An easy calculation, which 
we shall omit here, gives the equations of the harmonic polar of the tangent 
plane with respect to the trihedron, namely, 


2=— (5/12)y:— (5/12)¢:1, 
which represent the first axis of Bompiani. 

In the similar way the three conics in consideration intersect each other 
at three points distinct from O; when the quadric varies in the pencil of 
Darboux the plane of these points of intersection describes a flat pencil whose 
axis is given by 

(5/12) + (5/12)yy +1=—0, 


that is, the second axis of Bompiant. 


6. The curve of intersection of a surface and a quadric. As was 
shown in Section 3, the planes 


z= pi(y — nix) (t= 1, 2,3) 
determine uniquely a quadric whose equation is 
(36) ry—2z+ + + =0, 


provided that the directions n,, m2, m3 are connected by the relation (21). 
From (36) it follows that the quadric thus determined must have a contact 
of the second order with the surface S at O. 

By means of the relation (22) we can easily demonstrate that the triple 
tangents of the curve of intersection of (36) and S at O coincide with the 
three tangents 
(37) (y — mx) (y— nox) (y— ngv) = 0. 


Hence the quadrics corresponding to the planes through a non-asymptotic 
tangent t, or t. form a pencil of quadrics having a contact of the second order 
with Sat O and meeting S in a one-parameter family of curves each containing 
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these tangents as two of its triple-point tangents. The residual tangent t; at 0 
is the same for all the curves of the family and is the common tangent of all 
the osculating conics each of which lies on a quadric of the pencil. 


Conversely, suppose that a quadric Q has a contact of the second order 
with the surface S at O; the curve of intersection has O for its triple point. 
From (3) and (36) we have the equations of the triple-point tangents 


—- 0 (i = 1, 2, 3), 
where 71, %2, nz are the roots of the equation 
4B + kon + + 4yn* = 0, 
and therefore 


MN = — B/y, he + nong + gn), ks = — 4y(m + + 


In order to obtain the osculating plane of each branch at O we have to assume 
that the expansions for the branch which touches the tangent ¢;, 


z= y— = 0, 
are given by the equations 
(38) 
and then to substitute (38) into the equation 
4 (Bx? + yy®) + ba(x,y) + + + hye? = 0; 
the left-hand side should be satisfied identically up to and including the term 
in «*, Thus the coefficient A in (38) is given by 


(39) i) — — (1/9n4*) (B 
= (A/n;) [ks — (1/3?) (B — 2yni*) |. 

Observing that the osculating plane of the branch (38) at O is 
(40) Az=ni(y— 
and comparing (39) with (17), we are led to the following result: 

On every quadric Q having a contact of the second order with the surface 
S at O there are, in general, three osculating conics of the sections of S at 0 
and their plaies coincide with the osculating planes to three branches of the 
curve of intersection of S and Q. 


From this we infer that the characterization of the quadric of Wilczynski- 


rm 
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Bompiani stated in Section 3 coincides with the one obtained by L. Green.® 
Similarly, we can put the interesting result of L. Green into a new form con- 
taining the osculating conics of sections of the surface. 

Several exceptions to the last theorem shall be mentioned. 

(i) If the quadric Q determined by an osculating conic C' coincides with 
the Moutard quadric belonging to the tangent of C at O, then one of the two 
other osculating conics must necessarily coincide with C and the remaining 
one is the residual conic of Q. 

(ii) If the quadric @ is a non-Moutard quadric in a Moutard pencil, 
then the curve of intersection of Q with S at O has two coincident triple-point 
tangents, the osculating plane to each being the tangent plane of S at O. 
Therefore there is only one osculating conic on a non-Moutard quadric in a 
Moutard pencil as Pa has shown. 

(iii) If the quadric Q is the Moutard quadric belonging to one .of the 
Darboux tangents at O, there is no residual osculating conic. 

(iv) If Q is a non-Moutard quadric in the Moutard pencil belonging 
to one of Darboux tangents, it contains no osculating conic of the surface at O. 


NATIONAL UNIVERSITY OF CHEKIANG, 
TsuNYI, KWEICHOW, CHINA. 


°L. Green, loc. cit.*, p. 654. 
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THE ASYMPTOTIC BEHAVIOR OF A CLASS OF ENTIRE 
FUNCTIONS.* 


By C. V. NEwsom. 


Introduction. This paper considers the behavior of the function, 


co 8 
(1) f(z) = (n) = 1/ T(n+ 95), 
n=0 j=1 
for large values of the complex variable z, if s is an integer > 1, and p; isa 
complex constant. In particular, if k= [s/2]+1 is chosen as an odd 
integer, it is shown that 
(k-1)/2 


n 
(2) f(z) = S / 4 + w( Zp, 2) / 
j=0 


p=-(k-1) /2 
for all values of z satisfying the condition —2/s < argz < w/s; wherein cy = 1; 


a 
arg z= arg z— 2pr/s; p= pm; and lim w(2zp,n) =0. 
m=1 


is even the same conclusion results except that the limits upon the first 
summation are —k/2 and (k—2)/2 when — 2zx/s < argz < 0, and are 
— (k—2)/2 and k/2 when 0 < argz < 2z/s. 

The asymptotic behavior of the same function was considered by 
Van Engen? in his dissertation, but the method employed here is sufficiently 


general that it appears to offer an approach to the study of a large class of 


entire functions related to (1). 


1. A preliminary lemma. ‘The present study employs an important 
lemma which is essentially a modification of a previous theorem due to the 
author.” The earlier theorem and the present lemma are both concerned with 
entire functions of the general type (1) in which the coefficient g(n) is such 
that the function g(w) of the complex variable w—2x-+ iy satisfies the 


following two conditions: 


* Received January 11, 1942. 

1 Henry Van Engen, “ On the asymptotic behavior of analytic functions,” Disserta- 
tion, University of Michigan, 1934. 

*C. V. Newsom, “ On the character of certain entire functions in distant portions 
of the plane,” American Journal of Mathematics, vol. 60 (1938), pp. 561-572. Note, 
also, C. V. Newsom and A. Franck, “ Upon the asymptotic representation of functions 
of the Bessel type,” Bol. Mat., vol. 13 (1940), pp. 11-14. 
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(3) g(w) is single valued and analytic throughout the finite w-plane, and 


(4) |g(w)| < K e*, where K is a constant independent of x and y, and 


k is a positive integer. 


The theorem then states that 


x 


(5) {g(x) [+ 2*]*(sin / (sin x) }de— g(m)2™ + &(0,2), 


-c-1/2 


for all values of z satisfying the condition —a < arg (+ z*) <; wherein ¢ 
is any arbitrary positive integer; & is the integer appearing in condition (4) 
above ; the upper or the lower of the signs + is to be taken according as k is 


odd or even; and lim 2°&(c,z) = 0, irrespective of the value chosen for c. 


To inaugurate the study of the representation (5), let k be odd. In view 
g I 
(k-1)/2 
of the fact that (sin krx)/(sinwr) =1+2 > cos 2pre, it follows that the 
p=1 


integral of (5) may be rewritten in the form, 


(k-1) /2 
(6) {g(x) + 2 cos 


p=1 


wherein — < argz < arg = arg z— 2pm/s; and |z| =| 


Since 2e?™!” cos 2prx = 1 + e*?™'*, the expression within the bracket of 
the integrand of (6) may be considerably simplified. In fact, under the 


conditions upon g(w) previously given, we have the conclusion, 


(k-1)/2 


g(x) — g(m)2™ + &(c, 2), 


p=-(k-1)/2 m=-C 


in which < argz < arg % and |z| =| zp |. 
If & is even, a similar analysis is possible. The conclusion obtained 

differs from (7). however. in the fact that the limits upon the summation 

are and (k —2)/2 when — < argz < 0, and are — (k — 2) /2 


and k}/2 when 0 < arg z < 


2. Proof of the theorem. To employ the result (7) in the study of 
the particular problem of this paper, it is merely necessary to specialize g(x) 


| 
is a 
odd 
( 2p)" 
If k' 
first 
by 
0 


(11) log F(y)=-4$(1—s)log 2m + logs + as/y+ + + [an + a(y,n)]/f 
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as indicated in (1). As a consequence of a lemma demonstrated by Ford,’ 
condition (4) is satisfied when k = [s/2]+1. Thus, it is apparent that a 
study should now be undertaken of integrals of the type, 


oo 


(8) + 0s) Mae, 


-c-1/2 
subject to the restrictions upon z listed after (7). 


& 
If «= (2y— 2p +s+1)/2s, where p=} pm, the integral under con- 
m=1 
sideration becomes 


(9) (y) + 1) 


d 
in which, d= p—sc—s—1/2, and F(y) is the ratio of the functions, 


8 

T(y+1) and ss [[T(y/s + 1/2 + 1/2s—p/s+p;). The path of integra- 
j=1 

tion extends from the lower limit of the integral to positive infinity in the 


direction of the real y-axis. 
The analysis of log F(y) follows immediately when it is recalled that 


(10) log = log 27 + (y— logy—y + o(y), 


in which 
w(y)= B,/1-2y — (—1)"* [Bn + e(y, n) ]/(2n — 1) (2n)y™ 


and where lim e(y,n) =0, and |argy|<-. So, for sufficiently Jarge 
values of y taken upon the path of integration, and after considerable algebraic 


simplification, it is determined that 


in which the a’s are complex constants, and lim a(y,n)—0. Of course, then, 
ly|->00 


(12) = exp{as/y + a2/y? [an + a(y, ]/y"}- 

As a consequence of a well known principle in the theory of asymptotic 
series, the exponential factor occurring in this equality may be represented 
asymptotically by a power series in 1/y, in which the constant term is unity. 


3 W. B. Ford, Asymptotic Developments, Scientific Series, University of Michigan, 
XI (1936), p. 61. 
‘Ford, op. cit., pp. 3-4. 
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This series, in turn, may be reformulated into a factorial series so that integral 
(9) now becomes : 
n 
(18) (An) + +1) ]dy 
j=0 


d 


4 f [B(y, n) + n+ 1) Jay}; 
d 


wherein Cy = 1, and the remaining c’s are complex constants; and lim B(y, n) 


=0(. In fact, each integral obtained from relation (7) will be of this type. 
Integrals similar to those appearing in (13) were studied by Ford® in 
his consideration of functions of exponential type. In his study, however, 


it was necessary to make the restriction, | argz| </2. This condition is 


satisfied in the present analysis by all the z’s except possibly 24-1) /2 and 2-1) /2 
when & is odd, and 2x-2)/2 and 2/2 when k is even. Thus with the exception 
of these particular z’s, typified temporarily by the variable wu, the part of (13) 
inclosed within the brace may be written immediately in the following form 
as the result of the study of Ford: 


(14) + o(z, n)/(s2)"]; 
j=0 


wherein lim w(z,n) = 0. 


For the case when the real part of wu 0, a slightly different analysis of 
(13) is necessary. Referring again to the study of Ford, it will be observed 
that the first integral of (13) gives a result similar to (14) except that the 
remainder term does not involve the exponential factor, and is, in fact, in 


the form aj(u,n)/(su)4, where lim a;(u,n) Thus, for any z having 
j=0 |u|—>o0 


the property that | arg z| < 7/2, it follows that the product of this remainder 
term and the factor e-**(sz)™ produces a function which has the limit zero 
as |z|==|w|-—> co, when m is any positive integer. Consequently it may be 
combined with any (z,n) occurring in (14). A similar method serves to 
dispose of that part of the right member of equality (7) which is outside 
the integral. 

The second integral of (13), when the complex variable w is 2%-1)/2, may 
be written thus: 


5 Op. cit., pp. 65-67. 


n 
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(15) 1/(su) f [8(2—n—1,n) 


d+n+1 


= 1/(s%)"{ (2/u)” f [ P(x) (x) ]dz}, 
d+n+1 
in which P(r) = B(x —n-—1, n)e'*-)-*)ri/s, Of course, a similar analysis 
is possible for any other value of w. 

As the result of the application of another theorem given by Ford,° the 
function inclosed within the brace of (15), when multiplied by e-*, has the 
limit zero as | %|—> oo. Thus (2% ,) may be augmented by such a product. 
Consequently, result (14)' is appropriate irrespective of the particular z 
selected. It is true, however, that the expansion ceases to be dominant for 
any z having the property that its real part S 0. 

In summary, the representation of f(z) given in (2) follows immediately. 
In practice, k should be chosen as [s/2] +1, or perhaps, [s/2] + 2, if it is 
desired to enlarge the sector for z over which the expansion is valid. 


THE UNIVERSITY OF NEW MExIco, 
ALBUQUERQUE, NEW MEXxIco. 


° Op. cit., p. 38. 
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CONFORMAL GEOMETRY OF HORN ANGLES OF 
HIGHER ORDER.* 


sy Mary ELizABetH LADUE. 


Introduction. Under a conformal transformation a single curve has no 
invariant, but two intersecting curves have as an invariant the size of the 
angle between their tangents at their point of intersection. Since the order 
of contact is also a conformal invariant, two curves with first or higher order 
contact will be transformed into a pair of curves having the same order contact 
and the invariant mentioned above will have the trivial value zero. ‘This con- 
figuration, an ordered pair of analytic curves with n-th order contact, is a 
horn angle of the n-th order with its vertex at their point of intersection. 
Kasner has shown that the lowest order non-trivial conformal invariant of a 
first order horn angle is a third order differential invariant. In fact he has 
shown that every horn angle has one and only one non-trivial conformal in- 
variant which for an n-th order horn angle will be of differential order 2n + 1. 
This 2x +- 1 order differential invariant of the n-th order horn angle he calls 
the measure of that horn angle, and he has found the expressions for the 
measures of first and second order horn angles (1). In this paper the measures 
of third, fourth and fifth order horn angles will be obtained. 

The measure of an n-th order horn angle is determined by the values 
assumed by the first 2n + 1 derivatives of the two curves forming the sides 
of the angle at their point of intersection. Consequently its formula can be 
expressed in terms of the first 2n + 2 coefficients of the power series expansions 
of these curves (actually the first n + 1 of these cancel out), or in terms of 
the curvatures and are-length derivatives of curvature evaluated at the vertex, 
or finally, in simplest form, in terms of n+ 1 certain simple functions of 
these curvatures and their arc-length derivatives. For any particular curve 
the n+ 1 values assumed by these functions at any particular point can be 
regarded as the codrdinates of a point in an n + 1 dimensional space. If the 
set of all curves (analytic ares) through a given point and having there n-th 
order contact with each other is called an n-th order horn set, then with such 
ahorn set is associated the set of points corresponding to the curves of that 
set. This set of points is called the K»,,:-space associated with the given horn 
set. Since the distance between any two points of this space is given by the 
measure of the horn angle giving rise to these points, the study of the geometry 
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of this space is equivalent to the study of the conformal invariants of this 
set of curves. The properties of the geometries of the K-spaces associated 
with first and second order horn angles have been studied in considerable 
detail by Kasner, de Cicco, and Comenetz (2), and in this paper the geometry 
of a third order horn set is investigated in order that the similarities and 
differences in the geometries corresponding to various ordered horn sets may 
be considered. 

Roughly dual to the group of conformal transformations is the group 
of equilong transformations of the plane which convert straight lines into 
straight lines in such a way that the distance between the points of contact 
of any two curves (defined as the envelope of an infinitude of lines) having 
a common tangent remains invariant (3). When this distance is zero and the 
two points of contact are the same, the configuration becomes a horn angle 
which Kasner has shown to have a non-trivial equilong invariant of higher 
order (la). As in the conformal case, it is possible to associate with each horn 
set an equilong geometry which will have as its metric the equilong measure 
of horn angles of that order of contact. Kasner has investigated the equilong 
geometry associated with first order horn sets and found it to be identical 
with that of the conformal geometry of first order horn sets (4). In this paper 
the equilong geometries corresponding to horn sets of second and third order 
are investigated and, as in the earlier case, are found to be identical with the 
corresponding conformal geometries. The comparison of the conformal and 
equilong geometries of several higher orders is to be published soon. 

The problem of a conformal measure of a horn angle is interesting not 
only as the study of the invariants of an important geometry—the geometry 
associated with the group of conformal transformations of the plane—but 
historically in that it answers by modern methods a question the ancient 
Greek mathematicians were unable to cope with satisfactorily. For the con- 
formal measure of a horn angle provides a means of distinguishing according 
to size between angles which the Greeks were unwillingly forced to call the 
same, a measure furthermore which is significant in that it is invariant under 
a group of transformations which includes the Euclidean group of rigid 


motions of Greek geometry. 


1. Conformal measure of Horn angles of higher order. In this section 
the measures of horn angles of several higher orders of contact will be de- 
veloped. Since these invariants can be expressed in terms of the metric 
invariants of curvature and the arc-length derivatives of curvature of the 
curves forming the sides of the angle, they are obviously invariant under the 
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group of rigid motions. Hence there is no loss in generality in considering 
the measure of the horn angle whose sides are the curves 


(1) y= au? + 

which has the origin as its common point, the positive direction of the x-axis 
as its common direction, and has contact of order n at the origin. Further- 
more, this study will be limited to the effect on this horn angle of those direct 
conformal transformations which carry the origin of the (zx, y)-plane into the 
origin of the (X, ¥)-plane and the positive direction on the z-axis into the 
positive direction on the X-axis. Such a transformation can, in general, be 


represented by a power series of the form 


(2) Z=cz+e2°+--- (with c, real and positive), 
or +1Y, and = + by 


= + a -— y?) 

= + — y’) + +--- (a, >0). 
Under the transformation (2’) the horn angle (1, 1’) is carried conformally 
into another horn angle of the same order of contact which may be given by 
the pair of curves 


(3) Y=A,X?4+ A,V3+--- 


The Measure of a Horn Angle of the Third Order. The measure of a 
horn angle of the third order will be derived by determining an expression 
in terms of the coefficients of the curves (1) and (1’) which is invariant under 
the transformation (2’). This will be done by first considering the expres- 
sions invariant under the modified form of the transformation which is ob- 
tained by setting a, = 1, 


(2”) A = - — 2Bovy + 
Y=y-+ B.(a*—-y’) + 


This will cause a simplification in the calculations and a slight modification 
i the results which can easily be corrected later by following the transforma- 
tion (2”) by the simple magnification 


Z = (a, >0). 


The condition that the horn angle (1, 1’) be equivalent to the horn angle 

(3,3’) under the transformation (2”) is expressed by the two identities in a 
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+ B2[x? +- +- °° -)?] + + °° 


+ Bz [ + a’ -)?] + (a’ ov? a’ | 


obtained by substituting (2”) into (3) and (3’) and eliminating y by means f 
of (1) and (1’). The following equations, representing the effect of the j 
transformation (2”) on the coefficients of the curve (1), are obtained by { 
equating to zero the coefficients of the various powers of x in the identity (4) 


and solving the resulting equations for -,A;: 


(5;) As = 2+ 
(52) As + Bs — 4 
53) Ag + +- (as — + ) + Bs — +) 


( 

(54) As = — + + + 482”) az + 
+ (3B; — 64282) + (2a, — +- — + 


(55)? Ag -+ 


(5,5)? A;=4; — + 1322,‘ Bo. 


The corresponding relations (5’,),- - -, (5’g) for the second side of the hom 
angle are obtained from these by replacing ds,: --,d;, by 

In order to find the invariants, it is necessary to eliminate the coefficients { 
a:, 8; of the transformation. Since Aa, = Aa, = 0,’ the coefficients a4, %5, %, ; 
Bs, Bs. Be, Bz can be eliminated by subtracting equations (5,),° -, from | 
the corresponding primed expressions. The resulting equations are ’ | 


+, +, A’; respectively in the above equations. 


(6) AA, = Aa; 


; ( 
(7) AA; Aa; — 
(8) AAg = Ad, — + 10B.d2Aa, — (a — 402° — 482”) Aa, 
1 The complete expressions for A, and A; in terms of a,,. - -,a; and the coefficient 
of the transformation which involve, respectively, forty-five and eighty-three terms may > 
be obtained from the author. 


? Throughout this section of the paper the A notation will be used to indicate th® 
difference between corresponding expressions connected with the two sides of a hong 


angle; e.g, Aa, = a’, — a,, ete. 


neans 
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(9) AA; = Aa; — 4a,Aag + 128,a,Aa; — (2a; — 8a,” — 48) Aas 

(44.4; + 12B.B; 324,82") Ady. 
Next 8., 8; are eliminated by substituting from equations (5,) and (5,) and 
%, @; by substituting from (7) and (8) into (9). Using (6) and rearranging, 


the expression 


(10) AA;*— 2AA,AA,AA, + (AA, + 44,2AA,) — 124,A,AA8 
— Aa;? — 2Aa;Aa;Aa, + (Aa; + Ady? 12a,a,Aa,° 

is obtained. 

Equations (6) and (10) exhibit two quantities which are invariant under 
the special transformation (2”). Since the effect of the magnification (2’”) 
on the curve (1) is to divide the coefficient a, by a,"-1, equations (6) and 
(10) still hold under the general transformation (2’) provided the factors 
1/a,° and 1/a,** are inserted in the right members of (6) and (10) respec- 
tively. (a,>40 for 2,* is the Jacobian of the transformation (2’) at the 


origin.) Consequently, the quantities 
(6’) 
(10’) Aa;* — 2AagAa;Aa, + (Aa;+ 4a2?Aa;) — 12a,a,Aa,° 
are relative invariants of the horn angle (1,1’) under the conformal trans- 
formation (2’) of index 1 and 4 respectively, and 
THEOREM I. The quantity 


(11) 
Aa;* — (Aa; + Aas” 12a,a;,Aa,? 


and the sign of Aa, are the invariants of a horn angle of the third order under 
the conformal transformation (2). Conversely, if two such horn angles agree 
in the value of (11) and the sign of Aa,, they can be made to coincide by a 
conformal transformation (at least to any finite order of contact). 


For the proof of the converse, see Kasner’s Cambridge paper (la). The 
quantity (11) will be called the conformal measure of a third order horn angle 
and will be denoted by M,.*. It is a differential invariant of the seventh order. 
The coefficients of the curve (1) are expressed in terms of the values, 
at the origin, of the curvature, y,, and its derivatives with respect to length 
of arc, by the standard formulae 
(12, 
12 
(12, 


202 


9 


) 
) 6a; = 
) 24a, 3y1° 
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(12.4) 120a; = 19y,?y,’ 
(125) = yil¥ + + 48y1917 + 99y1° 

(126) 50400, = + 55 yi? + 199 + + 1135y,4y1’ 

and the coefficients of the curve (1’) by similar formulae in ye, , 2", 
Consequently, the invariant measure (11) can be written as follows: 

(13) = 

875 


— + (25Ay" + Ay”? —- 225 


This might be regarded as the intrinsic form for the conformal measure of a 
third order horn angle since it is invariant under any conformal transforma- 
tion, but it is an unwieldly expression and it will be convenient to use a 
simpler form. 
Let the equations 

w—a, 
(14.) 
(14;) 
(144) + 44275 — 
define a correspondence between the number quadruples (w,2z,y,z) and the 
coefficients a,,- - -,a,; of the curve (1). Then the conformal measure of the 
third order horn angle can be written in the form 

Aw* 
(15) 


where w,2z,y,2z and w’,2’,y’,2’ can be expressed in terms of the metric in- 
variants y;, y:’,° and ys, y2’,* by combining the equations (12) 
and (14). 


The measure of a fourth order horn angle. The derivation of the con- 
formal invariant of a horn angle of the fourth order is similar to the procedure 
outlined above for the third order case, presenting no new difficulties other 
than the increasing mechanical difficulty in handling the computation. The 
work can be somewhat simplified, however, by omitting the determination of 
the set of equations analogous to the set (5) and, instead, forming im- 
mediately a set of difference equations obtained by equating to zero the coefii 
cient of x” in the identity (4) and in the identity (4’), then subtracting the 
first equation thus obtained from the second, and finally solving the resulting 


difference equation for AA, (n= 


t 
a 
Ar® — 2AvAyAw + AzAw* a 
h¢ 
of 
(1 
(1 
(1 
(1 
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Note: These three operations, if done in this order, can be performed 
almost entirely by inspection of the two identities (4) and (4’) because 
the difference eqnations contain considerably fewer terms than the equa- 
tions (5) due to the fact that in a fourth order horn angle Aa, = Aa; 
= Aa,=—0. However, the use of this method does entail the loss of the 
explicit statement of the effect of the transformation (2”) on the coeffi- 
cients of the curve (1) which is needed in obtaining the group of induced 
transformations of a fourth order horn set and the expressions for Bs and 
8, which would be needed in the elimination process in the derivation of 


the measures of horn angles of order greater than eight. 


From these difference equations the coeificients a;, B; of the transformation 
can be eliminated as before. As a result of this elimination the quantities 


(16) Aa; 
(17) Aag* — 3Aa;Aa,’Aa; + (Aa;* + 2AasAag + Aas? 
—— (Ady + 4a,"Aa; + 12a,a,Aa,) Aa;* + (9437 + 16aea,) Aa;s* 
are obtained as relative invariants of the fourth order horn angle formed by 
the curves (1,1’) under the conformal transformation (2’). They are of 


index 1 and 5 respectively. Consequently, 


THEOREM II. The quantity 
Aa;" 


(18) — 
Aagi —- + + (9a,7 16a2a,) Aa;* 


and the sign of Aa; are the invariants of a fourth order horn angle under the 
conformal transformation (2). Conversely, if two such horn angles agree in 
the value of (18) and in the sign of Aas, they can be made to coincide by a 


conformal transformations (at least to any finite order of contact). 


The quantity (18) will be called the conformal measure of a fourth order 
hornangle and denoted by V/,.°. It is a differential invariant of the ninth order. 

Were it desirable, the quantities (16), (17), (18) could now be expressed 
in terms of the curvatures and arec-length derivatives of curvature of the curves 
of the horn angle as was done for the horn angle of third order contact in 
expression (13), but us this form is not needed, it will be omitted. 

The measure, J/,.°, can be expressed in its simplified form analogous to 
(15) by means of the correspondence 
(19.) v—a; 
(19.) W = (lg 


(19;) t= il; 
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(194) + 
(19;) = ily + + — (9a3” + 160204) a5. 


Then 


x 
| 


Aw* — 3AvAw?Aa + (Ac? 2AwAy) Av? — AvAz 


(20) = 


is the measure of the fourth order horn angle in terms of certain simple func- 
tions of the coefficients of the curves of the horn angle and consequently in 
terms of simple functions of the curvatures and arc-length derivatives of 
curvature of these curves at the origin. Thus M,,.° is a metric invariant also. 
The measure of a fifth order horn angle. Analogous to the results in the 

two preceding sections, the quantities 
(21) Adg 
(22) Aa;>— 4Aa;*AasAdg + (3Ad;*Ady + 3Aa;Aag? + 4a22Aa;*) Ade? 

— (2Aa;Adyo + 2AdsAdy + 8as2Aa;Ads + 12a.a3Aa;7) 

+ [Ady, + 4a22Ady + + (9a3? + 16a204) Aa; ] Aag! 

— (20a.a; + 240,05) 


are relative invariants of the horn angle of fifth order contact formed by the 
curves (1, 1’) under the conformal transformation (2’). They are of index | 


and 6 respectively . Consequently, 


THEOREM III. The quantity 
(23) 
—: (20a.a5 + 24a,4,) 


and the sign of Ad, are the invariants of a fifth order horn angle under the 
conformal transformation (2). Conversely, if two such horn angles agree w 
the value of (23) and the sign of Adg, they can be made to coincide (at least 
to any finite order of contact) by a conformal transformation. 

The quantity (23) will be called the conformal measure of the fifth order 
horn angle and denoted by V,.,°. It is a differential invariant of the eleventh 
order. 

Again the expression of these invariants in terms of the intrinsic metric 
invariants will be omitted. The measure, M,.°, can be expressed in the sim- 
plified form 
(24) M,.° = 

Au® 


— 4AuAv'Aw + 3Au?(Av2Ar + AvAw?)— 2Au? (AvAy + AwAr) + 
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by means of the correspondence 


WwW = sg 


Y = Ayo + + 


) 
) 
25.) t= ily 
) 
) + + + + 164204) a7 —( + 240304) 


By combining equations (12;) and (25;), the functions u,v, wv, x, y,z can be 
expressed in terms of the metric invariants, curvature and its arc-length 
derivatives, indicating that M,.° is a metric invariant also. 

Since the existence of one and only one such conformal invariant has 
been established by Kasner for horn angles of any finite order of contact, and 
since the methods employed for each case are such simple extensions of the 
methods of the previous case, it seems probable that a general form for the 
invariant for the n-th order case exists. So far, however, it has not been 
possible to find a symbolism in which to express this conformal invariant of 
differential order 2n + 1 which is sufficiently compact to permit a proof of its 
existence to be carried through by producing the measure. 


2. Conformal geometry of third order Horn angles. In an earlier 
section of this paper it has been shown that two given horn angles of n-th 
order contact can be transformed into one another by a direct conformal 
transformation if and only if the two horn angles agree in measure (and in 
the sign of Aa,,,) and that the measure of a given horn angle is completely 
determined by the values of the first 2n + 1 derivatives of the sides at the 
vertex. Consequently, the study of the conformal geometry of horn angles 
may be reduced to the study of the conformal geometry: of the differential 
elements of the (2n + 1)-st order passing through a given point in a given 
direction. In the remainder of this section the term curve will mean the 
differential element defined by a point on a given curve and the derivatives 
of this curve at that point, and the term n-th order horn set will mean the 
totality of all curves | (2n + 1)-st order elements] through a given point and 
each having at this point n-th order contact with every curve of the set. 
Again the n-th order horn set for which the common point and common 
direction are the origin and horizontal axis may be considered without loss 
of generality. 

Associated with each third order horn set (y,y’) (i.e., the horn set of 
curves through the origin having there slope zero, curvature y, and arc-length 


{| 
( 
( 
( 
( 
in ( 
of 
30. 
he 
ler 
th 
ric 
m- 


464 MARY ELIZABETH LADUE. 


derivative of curvature y’) is a four dimensional space, K4(y, y’), which con- 
sists of the totality of number quadruples (w,2,y,z) where w, x,y,z are the 
simple functions of the coefficients of the curves of the horn set (y, y’) given 
by the equations (14;;1). The conformal transformation (2;1) induces a 
transformation between the curves of the third order horn set (y,y’) and 
those of another third order horn set (T, I”) and, consequently, between the 
elements (w, x,y,z) of the space K,4(y, 7’) and the elements (W, X,Y, Z) of 
the space K,(T,1). The geometry of the K,-space associated with a given 
horn set (i.e., the invariants of this space under the transformation thus 
induced) will be the subject of this section of the paper. 


Definitions of the elements of a k-space. A point of the K,4-space is a 
number quadruple (w.«#,y,z) as defined by the equations (14;; 1). 

A series consists of any single infinitude of points of the space. 

A line (linear series) consists of the single infinitude of points of the 
space which satisfy three linear equations of the form 
(1) r—pwtn, y=pwtr 
where pi, 7; (11, 2,3) are constants. When two lines are considered, double 
subscripts will be used, the second subscript indicating to which line the 
coefficient belongs. 

A plane (flat congruence) consists of the double infinitude of points of 
the space which satisfy two linear equations of the form 
(2) z=a.w + bor+ 
where aj, bi,¢; (11,2) are constants. Second subscripts will be used to 
distinguish between the elements when more than one is considered. 

A flat consists of the triple infinitude of points which satisfy the single 
linear equation 
(3) 
where d, f, g, k are constants, subscripts indicating the flat to which they 
belong when more than one flat is considered. 

Derivation of the fundamental group, G;. The transformation induced 
by the conformal transformation (2;1) and carrying the point (w, 2, y,z) of 
the space K,(y, y’) into the point (W, X, Y,Z) of the space K,(T, I’) is given 
by the equations 

W=m,'w+h, 
(4) X = + 2m,?m.w + he 
Y = m,*y + + + hs 
Z = + + (1g — w + hg 


he 


‘le 
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where m, 0, ms, mz, hi, he, hs, hg are constants. The verification of this 
set of equations lies in expressing the functions W, X, Y, Z, defined as func- 
tions of the coefficients of the transformed curve (3;1), as functions of the 
coefficients of the curve (1;1) and then expressing this last set of coefficients 
as functions of w, z, y, z The effect of the transformation (2”;1) on the 
coefficients of the curve (1;1) is given by the equations (5;1) from which 
p, and B; can be eliminated as before by the use of (5,;1) and (52;1) 


leaving 


A, =, + f (G2, dg, Bi) 
As = — + (de, ds, Bi) 
= — 3a.a; + (— a3 + 4a,” + + — a, 

+ ds, Bi) 
A; $4.0 +- 2(— a; + 4a,” + + 2A.a,— 4a2*)a; 

+ [— 4a,(— as + 2a.? + 24.0. 

+ 12(A.A, — Jay 4 (ds, %i, Bi) 
where f, g, h, are functions of do, ds, %:, Bi (t= 2,:--,7) and do not 
involve *,d;. The effect of using the modified transformation (2”; 1) 
in place of the transformation (2;1) is eliminated by dividing each of the 
coefficients a, of the right members of equations (5) by 2," as before and 
the resulting equations, 


A, = (1 41° )ay f’ 


= — 2(a2/a;*)a, + 9’ 
+ [—4(2, a5) (— + 20,20? + 2a,7A,? + 2a, At. — 
+ 12(1/2,°) — azas) Ja, + 


express the effect of a more general transformation like (2;1) on the coeffi- 
cients of the curve (1;1). From these equations (6), the equations (14;; 1) 
and the corresponding equations expressing W, X, Y, Z as functions of 
A,,- + -,A;, are obtained the equations 
(1/a,3)w +f" 
X = (1/a,*)2 — 2(a2/a,°)w + 9’ 
“== (1/a,°)y —3(a2/a,*)@ (— + 42,72, 
+ + — 6a.")w + h’ 
Z = (1/a,°)z —4(a2/a,5)y + 2(1/a,°) (— + 
+ 44,°A.? + — 
— 4(a2/a,°) (— + + — 6a.")w 


+ kh’ + 44,°9’ — 124,A,f’ 
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which become the equations of the induced transformation (4) when the 
substitutions 

= 1/4, Mz = — Ms =(1/a,*) (— + + + 
hg =h’, hg =k’ + —12A2A5f’ 


are made. Consequently, 


THEOREM IV. The group of conformal transformations of the plane 
induces a seven parameter group, G;, of transformations between the points 
of the two K,-spaces associated with the given and transformed horn sets. 


Furthermore, Theorem I may be restated as follows 

THEOREM V. An ordered pair of points of a space K,(y,y') has the 
unique absolute conformal invariant, the measure or distance between the 
points, M,.* (1531). 


Thus the differential of arc in this space is 


(9) _dw* 
dw*dz- 2dwdxdy + 


The general line. A null series shall mean any series of the K4-space for 
which the measure of any pair of distinct points of the series is either zero or 
infinity. Thus from (9) it is seen that the null series are solutions of the 


differential equation (prime indicating differentiation with respect to w) 
(10) 2’ — + 2% =0 
and the special lines for which 
(11) w=constant. 
The line (1) is said to be a general line if 
(12) ps—2pipe + pr® £0. 


Geometrically, a general line is a line which is not tangent to a null series 


(i. e., the line and series do not have two consecutive points in common). 
Any general line (1) of a space K4(y, y’) is converted into a general line 


of a space K,(I,I”) under transformations of the fundamental group. This 
statement is proved by obtaining the effect of the transformation (4) on the 


coefficients of the line (1) which is given by the equations 


lane 
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= MP; + 2m2/m, 
2 = pz + 3mMmeop, + m3/m,? 
(13) Pz = + + 2msp,/m, + (4m2/m,*) (m3 — 
R, -—— myhypy + he 2mehy/m, 
Ry = + 3m,>mery — + hs — m3h,/m," 
R,; + + 2m,?msr, — ps3 — 4m, Mehr po 
— + hy (4m2/m,*) (mz — Ay. 
Then it can be seen that 
(14) P, SPF. +- (ps 2 Pip» p*), 
and consequently the line (1’) is a general line if the line (1) is. 
Given two general lines, /; and J., it is possible to obtain from the equa- 


tions (13) the following expressions from which the coefficients m2, ms, 
of the transformation (4) have been eliminated 


AP, = m,Ap, 
(15) AP.AR, — AP, AR, = m,°(Ap.Ar, — Ap, Ar.) 
(AP;AR, AP, AR;) 4AP.,(AP.AR, AP, AR,) 
= m,*[3Ap, (Ap,Ar, — Ap, Ar;,) — 4Ap.(Ap.Ar, — Ap, Ar.) ]. 
The left sides of these equations (15) and of the two expressions (14) stated 


once for each line are then relative invariants, consequently, 


THEOREM VI. Two general lines possess the four independent absolute 


invariants 


(16;) 7 
P3i — 2 Pej Prj Pas 


Ap .Ar, — Ap, Are 
(Psi — + Pur" 
3Ap, (ApsAr, — Ap,Ar,) — 4Ap.(ApoAr, — Ap, 

— 


(165) 


(165) 


Note: Two intersecting general lines possess only the two independent 
invariants (16;). 

The direct dihorn angle, aij, between the lines 1; and 1; is defined to be 
the quantity 

* Throughout this section of the paper the second set of subscripts will be used to 
indicate to which of two elements the coefficient belongs and the A notation to indicate 
the difference between corresponding coefficients for the two elements. 

*A single line possesses no absolute invariants but has the expression (12) as a 
relative invariant. 
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(Pris — Pri)® 
Psi — + Pas? 

and the reverse dihorn angle, x;;, for this ordered pair of lines, is defined by 
formula (17) with the i and the 7 interchanged. From the definition of these 
two angles and definition (15;1) the trihornometry (2b) of the third order 


(17) 


can be developed. 
Similar to the proof of Theorem VI is the proof of 


THEOREM VII. A point (Wo, Xo. Yo,%0) and the general line (1), not 
containing the general point, have the three independent absolute invariants 
g 


Ly — PiWo — 
(Ps — 2pipe + pr*)*”* 
2 (Yo — — — (Lo — Pro — 1) 
(Ps — 2pip2 + pr*)*” 
(2 — PsWo — 13) — 2ps (Yo — P2Wo — 12) + (3p1? — (29 — piWo—t 
(Ps— 2pip>2 + p*)* 


(18,) 


(18.) lpi, = 


(18;) [pi,= 


The extremum distance between the point and the line is a function of these 


invariants 


ws 


(18,) Extremum D = 


where W is that function of Jpi,,/pi,, pi, obtained by solving the cubic 


equation 
(19) W* — — 0. 
The general plane. The plane (2) is a general plane if 
Geometrically a general plane is a plane which is not the osculating plane of 
any null series (i.e., the plane and series do not have three consecutive 
points in common). ( 
From equations (4) it is possible to show that under the fundamental 
group, G;, any general plane (2) of a space K,(y,y’) is converted into a — ( 
general plane of a space K,(T, I’) 
(2’) Z=AW+BX+ 0). 


The transformation between the planes of the two K,4-spaces is given by the 


equations 
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A, = m,*a, — + (mz — /m,? 
B, = m,b, + 3m./m, 
Cy = m,°C, + (2meh, -— db, — my*hya, + hg — 3mzh2/m, — 
(21) A, = — 2m, Mob. + 4m, mea, — 8ms"b,/m, — 
B, = + 4m.b, + 2m;/m,? 
Cp = M,°C2 + Mm, — mhz) — + 4m,*mee, 
+ (4m2/m,) (2meh, — mh.) b, —- + hy — 2mzh2/m,? + 


From equations (2!) it is possible to obtain the following invariants: 
THEOREM VIII. A general plane (2) possesses the absolute invariant 


[ 27a, 18), (2a, bo) — 16b,° |? 


4[3(2a,— b.) + 46,°]* 


THEOREM IX. Two general planes possess, in addition to the invariant 
(22) stated for each plane, the four independent absolute invariants 


(23 ) [3 (2a). — boo) + 
— 180,, (26,, — — — — ba) 


— — ba) — 16b,,° — 4[3 (2a —bei) + 4b,,°]° 
(23. ) (2Aa, Ab. )* 
[2%a21 — 18b,1(2a,; — — 166,15]? — 4[3 — bo1) + 


23,) (Aa.Ab, — Aa, Ab.) 
— 1801: — — 160,,°]? — 4[3(2a11 — + 4b,,7]§ 


Note: In general two planes have only a single point in common, but if 


they have a line in common, the invariant (23,) vanishes. 


THEOREM X. A point (Wo, Xo; Yo, %) and a general plane (2) not con- 
taining the given point possess, in addition to the invariant (22) of the plane, 


the two independent absolute invariants 


(24, ) (Yo— — b,% — c,)*/* 
‘lg — 18b, (2a, — b2) —16b,°]* — 4[3(2a, — b.) + 40,7]? 
(24 ) (Zo — Ce) 4b, (Yo —— —-b,2% C1) 
— 18b, (2a, — bs) — 16b,*]* — 4[3(2a, — bs) + 40,2]° 
THEOREM XI. A general line (1) and a general plane (2) possess, in 
addition to the invariant (22) of the plane, the four independent absolute 


invariants 


18), (2a, be) 160,° 
Ps— 2pip2 + 


(25,) 
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(2b, — 3p,)* 

Ps 
(a, + — 
(Ps 2pipe + pr*)” 
(25,) [ (a2 + — ps) (C1 + 1101 — —(ai + pibi — po) (C2 + 
(ps — 2pipe + 


Note: If the line (1) lies in the plane (2), the invariants (25;) and 


(254) vanish. 

The general flat. The flat (3) is said to be a general flat if 
(26) 8d + 4fg + 9° #0. 
Such a flat is not the osculating flat of a null series (i.e., the flat and the 
series do not have four consecutive points in common) at a point for which 
g/2. 

From equations (4) it is possible to show that under the fundamental 
group, G;, any general flat (3) of a space K,4(y, y’) is converted into a general 


| 


flat of a space K,(T, I’) 
(3’) Z=DW+AFX+GY+ K. 


The transformation between the flats of the two K-spaces is given by the § 


equations 
D = m,3d —- 2m,mof — (mz — g/m, — (4m2/m,*) — 
F = m,°f —3m.g + (2/m,?) (mz — 6m,*) 
(27) G = mg + 4m2/m, 
K = m,°k — + m,(2m2h, — mhz) f 
+ — 6m,")hi/m, + mhs]g + hy 
—(2/m,*) (mz — he + (4m2/m,') (ms — hy. 


From equations (27) it is possible to obtain the invariants listed in the F 


following theorems. 


THEOREM XII. Two general flats possess the two independent absolute 


invariants 


THEOREM XIII. A point (Wo, Xo, Yo; 20) and a general flat (3) not con 
taining the given point possess the absolute invariant 1 


5A single flat possesses no absolute invariant but has the expression (26) as*— 


relative invariant. 
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256 (Zo dwy {xo k) 
(8d + 4f9 + 9°)* 


which is the extremum distance between the point and the flat. 


(29) 


THEOREM XIV. A general line (1) and a general flat (3) possess the 
three independent absolute invariants 
8d + 4f9 + g° 
(2p, g)* 
(305) ps — 4 — fp — 
8d + 4f9 + 9° 


Note: In general a line and a flat have a single point in common, in 


(30, ) 


which case 


(31) ps -— d — fp, — 


and no extremum distance can be computed. If in particular the line and 
the flat do not intersect (i. e., the line is “ parallel ” to the flat), the expression 
(31) vanishes but 
and the extremum distance is given by the expression 
256 (rs — k — fr, — gre) 

(8d + 4f9 + g°)° 


which in this case only is an absolute invariant. If on the other hand the 


(33) Extremum D = 


line lies entirely within the flat, both (31) and (32) vanish and the expression 


(33) gives the extremum distance, which is then zero. 


THEOREM XV. A general plane (2) and a general flat (3) possess, in 
addition to the invariant (22) of the plane, the three independent absolute 
invariants 
27s 18), (2a, b.) 16b,° 

8d + + 
(4b, — 3g9)° 
8d + 4f9 + 9° 
[3(2a, —-f) + b,(4b, — 3¢) 
(8d + 4f9 + 9°)” 


Note: As in Theorem XIV, no extremum distance can be computed if 


(34, ) 


(342) 


(34; ) 


the plane and flat intersect either in a line or in a point; but if they do not 
intersect 
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(35) —d-— ag = b. —f —big = 0 
and the extremum distance is given by the expression 
256 (co— — k) 

(8d + 4f9 + 9°)° 

which in this case only is an absolute invariant, or if the flat contains the 
plane both (35) and (36) vanish, the extremum distance being zero. 


(36) Extremum D = 


3. Equilong geometry of second and third order Horn angles. 


The equilong measure of second and third order horn angles. In this 
section of this paper the equilong measure and the equilong geometry of second 
and third order horn angles will be developed in order that the rough duality 


existing between the conformal and the equilong theories may be studied 


further. 
An equilong transformation of the plane is any contact transformation 


of the plane which converts straight lines into straight lines in such a way that 
the distance, 8, between the points of contact of any two curves having a 
common tangent remains invariant. Since the element is the straight line, 
it will be convenient to use Hessian line codrdinates (u,v) to indicate that 


directed line having a normal through the origin which is v units long and | 
makes an angle of w radians with the prime direction. In this codrdinate 
system a curve is given as the envelope of a single infinitude of lines and 
therefore by an equation of the form t 


v= + au? 
Two such curves will have contact of exactly n-th order at the line v = @ if 
the first n + 1 coefficients in their power series expansions are the same and 


the (n + 2)-nd are distinct. If these two curves have contact of first or higher FB a 
order, the invariant distance, 8, has the trivial value zero and the pair of ( 
curves form a horn angle of the given order of contact. This horn angle will 
have a unique equilong invariant. The form of this invariant for second and 
third order horn angles is obtained below. As in the conformal case there Cc 
is no loss in generality in limiting the study to the invariants of the horn 
angle 
v= + 7 
(1) y= (with Oni (5 
Since the group of equilong transformations of the plane can be shown to 
consist of the set of all line transformations of the form 
U = $(u) (67 


V = (u)v+ 
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where ¢ and y are arbitrary analytic functions (3), the equations 


= (a, ~0) 
represent an adequately general equilong transformation under which the horn 
angle (1) is equivalent to another horn angle of the same order of contact 
given by the equations 


V =A,U?2 + A,U?+--- 
V = + A’,U* (with A; An A'ny Anes 


Again the two identities 
(4) (a, + 2a.u-+- - +) + agu®+---)+ Biu+ Bou? +: 
= Ao(au+ au? au?+---)+- 
(4) (a, + 
= A’,(a,u + au? +--+)? + 


are the conditions for the equivalence of the horn angles (1) and (3) under 
the transformation (2). As before, when for a horn angle of a given order 
the coefficients of like powers of wu are equated to zero in each identity, dif- 
ference equations formed from the corresponding pairs of equations, and the 
constants a; and £; of the transformation (2) eliminated, a pair of relative 
invariants are obtained for this horn angle. For the second order horn angle 
the relative invariants are 

a Ads of index % 

(5) 


Aa,* — AazAas of index 
and for the third order horn angle 


(6) 


of index 
2Aa,Aa;Ady + Aay*Aa; of index 


Consequently, 


THeoreM XVI. A horn angle of second order has the unique absolute 
equilong invariant 
2Aa,A, 
and a horn angle of third order has the unique absolute equilong invariant 


Aa;* — 2Aa,Aa;Aa, + Aa,?Aa; 
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These quantities (5’) and (6’) are called the equilong measures of second and 
third order horn angles respectively and are denoted by py2* and pis’. 

As in the conformal theory, the invariants can be expressed in terms of 
intrinsic metric invariants, in this case in terms of the radii of curvature and 
their derivatives with respect to the inclinations of the curves.° However, 
it is more desirable for comparing the equilong and conformal measures to 
leave the equilong measures expressed in terms of the power series coefficients, 
regarding these as functions of the metric invariants of radius of curvature 
and its inclination-derivatives. Then the comparison of the equilong measures f 
(5’) and (6’) with the corresponding conformal measures ° and 
M,.* (15;1) reveals that the two measures for second and third order hor 
angles are identical in algebraic structure. 

The “duality” of the conformal and equilong geometries of second and 
third order horn angles. In section 2 of this paper the transformation induced 
in the space of elements (w,2,y,z) associated with a third order horn set 
by the conformal transformation (2;1) has been discussed in detail. Thef 
conformal geometries associated with first and second order horn sets has also 
been discussed thoroughly by Kasner and de Cicco (2). Furthermore Kasnerf 
has shown that the equilong geometry associated with a first order horn set 
is identical with the conformal geometry associated with that horn set (4b). 
In this section of this paper the extension of this “ duality ” to the geometries 
of second and third order horn angles is established. 

The conformal K,-space associated with a given second order horn sé 


may be defined as the space of elements (2, y.z) given by the equations (2c, d) 


1 = ils 
(7) 
as 40273 


where the a; are the coefficients in the power series expansion (1;1) of the 
curves of the second order horn set. For this space the fundamental trans ; 
formation induced by the conformal transformation (2;1) of the Cartesian 
plane is given by the equations * 


. ar’ 
5Ar’?— 4Ar’Ar’” + 4Ar”?’ 
875 Ar’ 
4 — Ar’ Ariv + (25A0v + 45Ar’") Ay? — 


7 Equations (7) and (8) differ a little from the definition and transformation @ 
the K-space given by Kasner and de Cicco because the author found this definitio FF 
more convenient. The effect on the geometry is unimportant. 


d and 


ms of 
e and 
vever, 
res 0 
‘ents, 
rature 
and 


horn 


and 
cluced 
rn set 

The 
s also 
casner 


rn set 


(4b). 
retries 


rn set 


2¢, d) 


of thet 


trans 


‘tesial 


— Ark 


tion 
finition 


CONFORMAL GEOMETRY OF HORN ANGLES OF 


HIGHER ORDER. 


X=mr+h, 
(8) Y= + + 
Z = + 2m,?moy + + hz 


and the fundamental invariant distance under this transformation (i. e., 3/,2°, 


the conformal measure of a second order horn angle) by 


(9) 


If the equilong K’,-space associated with the second order horn set is 


defined as the space of elements (2, y,2) given by the equations 
(10) Yay, 2—As, 


where the a; are the coefficients in the power series expansion (1) of the curves 
of the horn set in Hessian line codrdinates, then the transformation induced 
in this equilong K’,-space by the equilong transformation (2) of the Hessian 
plane is also given by the equations (8) and the fundamental distance in- 
variant, the equilong measure of a second order horn angle, by equation (9). 
Thus the * duality ” extends to second order horn angles. 

In section 2 of this paper the conformal K’4-space associated with a third 
order horn set was defined as the space of elements (w,2, y,z) given by the 
equations (14;;1) with an induced transformation given by the equations 
(4;2) and a metric given by the expression (15;1). In this paragraph the 
corresponding equilong A’,-space will be defined and shown to be identical 
with this conformal A’y-space. The equilong K’,-space associated with a given 
third order horn set is defined to be the space of elements (w,.,y,z) given 
by the equations 


where the a; are the coefficients in the power series expansion (1) of the curves 
of the horn set in Hessian line codrdinates. The development of the equations 
of the transformation induced in this space by an equilong transformation (2) 
of the Hessian plane follows the same procedure as that outlined on pages 18 
and 19 for the conformal case. Given a third order horn angle expressed in 
Hessian line coérdinates in the form (1), it is possible to obtain from the 
identity (4) the set of equations 


A, =(1/2,*)a, + hy (do, ds, 

(12) As =(1/a,*) a5 — + he (de, ds, 
Ag =(1/a,°) — 3 (a2/a1%) as —(1/1) (a — + hs (de, ds, 
A; =(1/a,°) a; —- ag — 2(1/a1?) (43/015 — 


) ( — +- h 4 (a2, As, a: ) 
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(where h,,- - are functions of dz, ds, and do not involve -,a;), 
These equations, which express the effect of the equilong transformation (2) 
on the coefficients of the curve (1), become the equations of the desired trans- 
formation by means of (11) when the substitutions 


(13) m,=—1/a,, [a3/a,° — 4(a.?/a,°) | 


are made. It is only necessary to write out this transformation in order to 
see that it is given by the equations (4; 2) and consequently has the expression 


(15; 1) as its fundamental invariant. Consequently, 


THEOREM XVII. The equilong geometry of a second or third order horn 
set is identical with the conformal geometry of that horn set. 


BARNARD COLLEGE, 
COLUMBIA UNIVERSITY. 
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ON THE SINGULARITIES AND THE QUALITATIVE BEHAVIOR 
OF THE SOLUTIONS OF THE PROBLEM 
OF CENTRAL FORCES.* 


By E. K. HAvILANp. 


Introduction. The purpose of the present note is to investigate the 
motion of two bodies under the influence of their mutual attractions, when 
these latter satisfy certain specific conditions. By choosing the origin of the 
codrdinate system so that it coincides with one of the particles, the motion 
may be reduced essentially to that of a single particle under a static central 
force. The force function, U, will in addition be subjected to the requirement 
that it be homogeneous of degree A (40) in the Cartesian codrdinates x and 
y. In the particular case where 4 = 0 the latter requirement is replaced by 
the condition that the first partial derivatives Uz, Uy, be homogeneous in z, y 
of degree X—- 1 = —1. This leads in the former case to a force function U 
proportional to r* and in the latter to one proportional to logr, where the 
constant of proportionality is chosen so as to correspond to an attractive force.* 

While a number of the results concerning the nature of the paths in 
various particular cases have long been known, the present paper is concerned 
primarily with an aspect of the subject which appears not to have been treated 
in detail hitherto, viz., the analytic character of the collisions, a subject which 
appears elementary at first glance, but which turns out to be complicated 
enough to involve a rather deep analytical apparatus; in particular, the proof 
of the non-existence of an analytic continuation in certain cases will be made 
to depend on Bohr’s theory of almost-periodic functions of a complex variable. 

At the same time, it has seemed desirable to complete and systematize 
the description of all the various cases arising from force functions of the 
above type, and in this connection the most striking result is that in the range 
—2<A<0 the situation is the same as for the Newtonian exponent 
A=—1; in the sense that the non-vanishing of the constant angular mo- 
mentum, c, remains necessary and sufficient for the non-existence of collisions 


—. 


* Received March 25, 1942. 

*Cf. W. Whewell, On the Free Motion of Points, Cambridge (1836), especially 
pp. 39-79. E. T. Whittaker, Analytical Dynamics, Cambridge (1917), pp. 77-86. P. 
Stickel, Encyclopidie der mathematischen Wissenschaften, vol. 4, pp. 494-496. 
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and singularities, and that in addition the three subcases belonging to a fixed 
c 0 are qualitatively the same as in the classical case. 
1. Classification of cases. Turning now to the detailed treatment of 


the problem, we observe that, if polar codrdinates be introduced, the equations 


of motion possess the integrals 


(1) + —U(r) =h 
and 
(2) rd’ =C (’=d/dt), 


the former representing the conservation of energy and the latter the con- 
servation of angular momentum. Since ¢ is an ignorable codrdinate, (1) may 


be replaced by 


(1 bis) U(r) + =h 
i.e., by 37°? U* =h, where U* = U*(r;c) = U(r) — $c?/r’, so that 
dr’? == U(r) —4e*?r? +h = say. 


We proceed to deduce from (1 bis) and (2) the nature of the motions 
corresponding to various forms of U(r). Since the number of possibilities 
is rather large, it has been thought desirable to summarize the various cases 


and subcases in the following table: 
I c=0. (A) 0<A<+ 0 
(a) h>0O. 


1. A irrational. 2. A rational = p/q. 


1*. p odd, 2*. p even = 2 


~K, 


2,*. «+ q odd, 2.*. 


(b) h=0. 

(c) h< 0. 
(B) U=—-.logr 
(C) V= 


(a) O<A<1. (B) 1SA<+ 2. 
> 


0. (a) h>0. 


(b) h=0. 


‘2 
wl 
(3 
(b) h<0. 
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c#0. (c>0). 
(A) U=r°. 

(a) O<A<2. 

(a) h>0O. (b) h=O. (c) h<O. 
(i) F(r) > 0, (ii) =9, (iii) < 0. 

(8B) \=2. 
(aa) €>0,h >0; (ab) €>0,4 <0; (ac) =0. 
(ba) £=0,h > 0; (bb) £=,0,h =0, (be) £=0,h < 0. 
(c)f<0. 

(y) 2<A< 7 oO. 


(a) h>0. 


(i) <0. ) (iz) 
(ii) F (ro) =0. Ly. A irrational, 2. A 
(iii) > 0. J rational = p/q. 


1*. p odd. 

2*, p even = 2k. 
i*.« + odd. 
o*. «x + q even. 


(b) ASO. 
(B) U =—logr. 


(i) h < loge+ 4. (ii) h=loge +4, (iii) h > loge + 


(C) 
(i) <0. (ii) F(t) =0, (iii) > 0. 
2. The case of vanishing angular momentum. For convenience in 
treating the individual cases, we have written (1 bis) in the form 
(3) == (r) +h=—F(r,h) = F(r) 


whene 


(3 bis) — f + | F(r, h)|-4dr, (t) = arb. const.). 


r 


We shall first suppose that 
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(I) c=0; (A) U=r’', (0<A<+ 0); and (a) h>0O. 

Under the assumption (1), the motion is rectilinear, so that we may write 
x (> 0) for r, and (3) then becomes 

(3 bis) ho F(z), 

whence 


e 0 


0 
As «0, t-— 0, so that a collision occurs, corresponding to 0. To in- 
vestigate its nature, we consider only values of z for which 0 < ha < p, say, 
(ep <1). Then the parenthesis in the preceding integral may be expanded 
by the binomial theorem and the resulting series integrated term by term, 
leading to the equation 
(27(A + 2) —h/(3A+ 2)a% 
+ 3h*/[4(5A + 2) — + 2) +--+}. 


fegarding now. ¢ and as complex variables, we set Then 


Q-M (42) 9 + 2) —h/(8A + 2)y 


+ 3h2/[4(5A + 2) Jy? — 5h8/[8(7A + 2) 
wr, t?\/(+2)- + » denoting those branches which are real and positive for real 
positive values of the argument. Thus z= y%(y), where %(y) is a regular 


analytic function of y in some neighborhood of y=0 and &(0) #0. Con 


sequently, 2’(0) +0, so that y is a regular analytic function of z in some 
neighborhood of z 0 and may be developed into a power series there, and 
since the coefficients of the series for z are real, the same is true of the series 


On expressing y,z in terms of 2, ¢ respectively, this leads to 


for y. 
(4) «= (1+ Bt’ + +--+), where B, > 0 and v= 2A(A + 2), 
the Puiseux series in parenthesis converging uniformly for all ¢ such that 
lt | < pi, say, pr > 0. 

Two cases arise: 

1. A is irrational. Then 2/(A+ 2) is also irrational, and 2 has 3 
logarithmic branch point at ¢— 0, so that x(t) possesses no analytic con 


tinuation through ¢ = 0. 
To prove this, we observe that, if we start with that branch of « which 
is real and positive for real positive values of ¢ (which is seen to be possible 


from the fact that the coefficients in the expression for z are real and B, > 9); 
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then since 2/(A + 2) is irrational, it is clear that no number of circuits about 
{= 0 will restore the first factor on the right of (4) to a real value. If, then, 
z is to assume its original value, it must, in particular, be real and this is 
possible only if, after n circuits, the second factor is proportional to the con- 
jugate of the first, with a real proportionality factor k. Put t= pexp(wW). 
Then after n circuits 


exp (iden /(A + 2) ) exp (i2xa), 


where = (2n/(A + 2)) = 2n/(A+ 2) — [2n/(A-+ 2) ]. 

The second factor on the right-hand side of (4) would then have to be 
of the form k'p?/“*exp(—i2ma) for all p < p:;. But this is impossible, since 
inspection of the factor in question shows that, for p sufficiently small, it lies 
in an arbitrarily small neighborhood of t= 1 and hence in the interior of 
the angle bounded by the half-rays 0 = + 27a. From this contradiction, it 
follows that x has a logarithmic branch point at ¢ = 0. 


2. A is rational, say of form p/q, where p and q are relatively prime. 
Reasoning similar to that of case 1 shows that z is real if and only if ¢?/*?) 
isreal. Therefore, x will be continuable through ¢ = 0 under and only under 
conditions such that ¢?/*?) is real for negative ¢. 

This last factor may be written as ¢74/'»**@), where 2q/(p + 2q) is to be 
interpreted as reduced to its lowest terms, viz., 2¢/(p+ 2q), if p is odd; 


q/(4p + q). if pis even. Then 


(1*) If p is odd, p+ 2q is odd and real analytic continuation is always 
possible, g arbitrary. This includes the Newtonian case p=—gq —1. 

(2*) If p is even, put p= 2x. Then 2/(A-+ 2) becomes g/(«-+ q), and 
(2,*) If «+ q is odd, real analytic continuation is again possible; 
(2.*) If «+ q is even, g must be odd, since otherwise we should not 

have p and q relatively prime. Then « is necessarily odd. In 

this case, real analytic continuation is impossible, e.g. p= 2, 


If the positive square root be chosen in (3 bis), it is seen that x increases 
steadily with ¢. As and x~ (2h)4t. Hence the solution is 
not bounded. If there are in general to be equilibrium solutions or asymptotic 


solutions, we must have“ for some ry 


(5) F (ro) = 0, 


*Cf. equation (3) above and A. Wintner, Analytical Foundations of Celestial 
Mechanics, p. 132 and p. 153. 


— 
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or in the present case, F = F,—0. But from (3 bis) it is seen that neither 
F nor F, vanishes for 0 < « < + o, so that under the present assumptions 
neither equilibrium solutions nor asymptotic solutions exist. 

Suppose next that 


(b) h=0. Here the motion is the same as in (a), with respect both to 
boundedness and to the nature of the collision and to the non-existence of 
equilibrium or asymptotic solutions, but the results are more easily obtained, 


since 


t = 24 f = + 2) or const. 


Again, suppose that 


(c) h<0. Then the nature of the collision is as before and there are no 
equilibrium solutions or asymptotic paths, since PF, ~0, + 0), 
but x ~ (2h)#t implies that the solution is bounded. 

Let assumption (I) be retained, but (A) be replaced by 


(B) U =— logr, (r= 2), 
so that $2" = F(x) = — logr +h; Fx, =—2"; Fre Thus there are 
no finite positive zeros of either F, or F.2. For all h ~ 0, F(0+) =+ 0, 
F(+ 0) =——o,. The path is rectilinear and has a maximum at «=e! 
corresponding to a cusp.* 

Since 
(6) dt = 2-4(h —log x) “dz, 


it follows that, if ¢* be the time arbitrarily corresponding to the position 2%, 


(0 < < e*), then 
(7) — t* == 2-4 f. (h — log 


from which it is seen that the solution path, starting from 2*, reaches both 


x= 0 and « =e" in finite time. To =0 may be made to correspond = 0, 
these representing the place and the time of collision. 

A necessary condition that x(t) be continuable analytically through 
¢t = 0 is that it should have at most an algebraic singularity there, whereupon 
the same would be true of the inverse function ¢=¢(x) and hence of the 
derivative of the latter, since differentiation cannot introduce transcendental 
singularities. But (6) shows that this derivative has a transcendental (loga- 


3 Cf. A. Wintner, op. cit., §§ 185, 186. 
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rithmic) singularity at «0, from which it follows that the analytic con- 
tinuation of = through ¢ = 0 is impossible.‘ 

As in the preceding case, there are no equilibrium solutions and no 
asymptotic paths. 

Continuing the assumption c = 0, we now suppose 


(C) U=U(r) =—2, 

(2) < }, 

then F(z) =—a2+h; Py Frg=—A(A—1)2*, 80 that F; 
and F,, are ~0 for all 7, O0<4< + 0. Suppose further, that 

(a) h> 0. 

Then F(0) =h > 0, F(+ ©) =— o, so that there exists a bounded colli- 


sion path, the maximum value of x corresponding to a cusp. As before, there 
are no equilibrium solutions or asymptotic solutions. To examine the nature 
of the collision, we consider 

(8) x’? = 2(h 


Taking, as we may, the positive square root 2’ = (2h)#(1—2*‘/h)3, and 
setting a/h'/\ 1, we may expand the parenthesis by the binomial theorem 
and the resulting series converges uniformly for | r’ | < p, say, (p <1), and 


may be integrated there term by term, leading to the expression 


3 


(k = const. ~ 0). 


Regarding ¢ and r as complex variables, we put 7 = y, so that 


3 


z=? and yr denoting those branches which are real for real values of 


the argument. Then, as in case (Aa), z= yM(y), where M(y) is a regular 


‘This fact may also be verified by making the substitution log « = y, followed by 


y= h—2*, whereupon (7) becomes t — t* =— V2 exp(h ) exp (-—2°)dz. If z=2(t) 

represents the inverse of the foregoing, then @ is given as a function of t by 
*=exp(h—ez*(t)). Since the nature of the transcendental singularity of t = t(z) 
at z= 00 is quite different from that of y=loga at r=0 (or =), it follows 
that as a function of ¢ has a transcendental singularity at t= 0. For the nature 
of t(z) at z= ~™, ef. Polya and Szegi, Aufgaben und Lehrsdtze aus der Analysis, 
Berlin (1925), vol. IT. p. 33, no, 189. 
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analytic function of y in some neighborhood of y= 0 and (0) #0, so that 
y may be expressed as a regular analytic function of z in some neighborhood 
of z= 0, and this leads to an expression of the form 


where y: > 0 and y, 0, (k 22), whereupon reasoning similar to that of 
case (Aa), shows that 


1. If d is irrational, no analytic continuation is possible. 

2. If A is rational, = p/q say, where p and q are relatively prime, we 
write (k —1)p/q = m/s, where it is understood that m and s are 
relatively prime. Here real analytic continuation is possible, if s is 
odd, but not if s is even. 


If, on the other hand, 
(b) h=0, 


it is clear from (8) and from the expression for FY, that there is no motion 
at all. 

Finally, if 
(B) 1<A<+~o 


and (a): h > 0, there exists, as before, a bounded collision path, the situation 
with respect to the behavior at the maximum value of z and with respect to 
the analytic nature of the collision being precisely as in case (a). In the 
cases in which analytic continuation through the origin is possible and A >1, 
the motion is similar to the case A = 2 of the harmonic oscillator. Since F 
and F’, do not vanish simultaneously, there are no equilibrium or asymptotic 
solutions. 

If, however, (b) : A = 0, then there is no motion except, in the case h = 0, 
for the trivial equilibrium solution r= 0. 


3. The case of non-vanishing angular momentum. Having considered 
linear motions, we now turn to the case 


(II) c~0, 


where we may suppose without restriction that c > 0. As our first sub-case, 
we shall consider that in which 


(A) U(r) =r* and (a) 0<A< 2. 
Here (3) becomes 


(9) hr”? — +h Fir), 


that 
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whence 


(10) Fr(r) =—Ar** + and (11) Fy, (r) =A(A+ 1)r*? — 


As r—>0-+, 1”? ~-—c?*/r?, which is impossible for real r, so that there 


js no collision in any of the three cases h = 0. 

Asr—>-+ o, 1’ ~ 2h, so that if (a): h >0, then r~kt, k a constant, 
which may be supposed positive. In this case, a solution exists and is not 
bounded. dr/dt has just one zero, since F’, has but one positive real zero and 
F crosses the r-axis an odd number of times, therefore once. The zero of 
dr/dt corresponds to a minimum value of r, which we may suppose to occur 
at ¢==0. For the nature of the path at this point, cf. the reference under 
(i) below. By (2), ¢’=c/r’, > 0), so that increases steadily with ¢. 
Moreover, as r ~ ki, integration of (2) shows that ¢ approaches a finite limit 
as ¢ > + oo, and a similar result holds if > — o. 

In this case, = would imply =0 also. But 
the latter two equations imply A = 2, contrary to hypothesis, so that no equi- 
librium or asymptotic solutions exist. 

Again, if (b): h 2(dt/dr)? so that r~ const. t?/*?), Here 
r becomes infinite more slowly than in case (a), but the motion is otherwise 
the same. 

Thirdly, if (c): h <0, any possible solution will be bounded. To de- 
termine whether or not a solution exists at all, we observe from (11) that 


3c? 1/(2-A) 


Consequently, F'(7) is first concave and then convex from below. The ex- 
tremum of F(r) is located at where =0. For 
0<A< 2, it is seen that 7) <7, so that the extremum is a maximum. 


F(0+) =— F(+ 0) <0. Therefore, 


VIIA 


F,,(r) = 0, according as r 


(i) if P(r.) > 0, there will be a bounded (periodic) non-circular motion. 
for a more detailed discussion of which the reader is referred to A. Wintner, 
The Analytical Foundations of Celestial Mechanics, § 215, p. 153; 


(ii) if F (7) 0, there will be a circular motion and nothing else; 


(iii) if F(7)) < 0, no solution will exist. 


In case (i), 0 << Ryo Ry, say, so that from (2) | ¢—¢ | = | t |, 
and ¢ becomes positively or negatively infinite with ¢, and the same is true 
in case (ii). 
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F and F, vanish simultaneously in, and only in, case (ii). The circular 
path then represents an equilibrium solution, but there are no asymptotic paths, 


Under the assumptions (II): ¢c > 0 and (A): U =r”, we now consider 
case 
(8) A=2. 

We set 1—4c? =, and examine first the subcase corresponding to 


(aa): £>0;h>0. Then 
(12) ¢/r’), whence (13) = + ht)3, 


if t= 0 corresponds to r = 0, and conversely. In this case, there is a collision, 
but no real analytic continuation of r through ¢ = 0 is possible. Substitution 
of (13) into (2) and subsequent integration leads to 


— = c[2(2E)8]* log[t/( (26)* + ht) ], 


so that @ becomes logarithmically infinite as t—+0-+. As ¢ increases and 
becomes positively infinite, the same is true of r. @ also increases steadily, 
but approaches a finite limit as /—> + o. 

If, on the other hand, (ab): h < 0, the situation with respect to the 
collision as t—>0-+ is the same as before, but as ¢ increases, 7 remains 
bounded. ¢ increases steadily with ¢ and becomes infinite as t —(2£)4/h > 0, 
while as {> 0+, @ becomes negatively infinite. 


In both the foregoing cases, = r cos @ becomes 
= 243( (2¢)4 + ht)} + ht) ]} 
and the expression for y =r sin ¢ is similar. 


Thirdly, if (ac): h =0, we have 1’? = 2¢/r?, so r= (8¢)4#4. Here, as 
before, there is a collision, and real analytic continuation of r through t=? 


is impossible. = c/r? = c/((8Z)4t); ¢=c(8f)4log¢-+ const., so that 
at as ow, in which latter case 


r—+>-+ also. 
Again, if (ba): £=0, h >0, (12) reduces to 1”? = 2h, so r= (2h)¥ 
In this case, 7 can be continued analytically through t—0, ¢’ =c/2h#; 
= — so that as > 0 +. becomes negatively infinite (assuming 
c > 0), but more strongly than in the preceding three cases. As 
r—>+to and ¢6>0—. x= (2h)4t cos (c/2ht), y= — (2h) 4t sin (c/2ht). 
If (bb): £=0, (12) reduces to =0, r—p=const. and 
d = cpt em where the latter constant may be taken to be zero, and 
we have circular motion. 
If (be): £=0, h < 0, there is no real solution. 
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Finally if (ec): £ < 0, there is no collision for any value of h. IfhS0, 
there is no real solution at all, just as in case (be). If, on the other hand, 


n> 0, 
t 2-3 r(hr? + €) 4dr = (hr? + £)2/h23, 
R 
where is the minimum value of r, corresponding to 0. For the behavior 


of the path there, cf. the reference under case (II Aaci). r?=h-1(2h7t?7—£), 
so that upon substitution in (2). and integration, ¢ = cn/(2h) arctan (yt), 


where 4° = 2h*/(— £), from which it appears that @ approaches a finite limit 
as and as o, while in the latter case, r++ o. 

The usual tests applied to the case A = 2 show that there are no equi- 
librium solutions except in the case (bb), where F =0 and F,=0, in which 
ease all circular motions are equilibrium solutions. In no case do there exist 
paths asymptotic as t—> o. 


We next assume (II): c>0; (A): U=r'; 


(y) 2R<A<+ 0. 

In this case, 7, < ry) (cf. p. 485 for the notation), so that the extremum 
of F(*) is a minimum. 

We first suppose (a): h > 0, so that it is seen from (9) that F(r) tends 
to + « as r—>0-+ and is positive as r—>-+ o. Then according as 

(i) F(%) <0, (ii) F(t.) =9, (ii) >9, 
there will be 
(i,) a bounded collision path and (i,) an unbounded path without collision ; 
(ii) a circular equilibrium path and asymptotic paths corresponding to (i,) 

and (is) ; 
(iii) an unbounded collision path. 


We examine first case (i,). Here 


== — + h, (2<A<+ 0), 
so that 


(14) ; r/?{1 + (hr — $c?r-?) } adr, 


(15) 


0 


ar 
J yr A (hi +. 36 (hr* — ‘}dr. 


The series in brackets converges uniformly and absolutely for | -? | ($c? 
+|hr? |) <p, say, (0<p<1), i.e. for |r| < r**, for a suitable r** > 0 
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Consequently, (15) may be integrated term by term there and the resulting 
series will converge uniformly for r in the above range. Before carrying out 


this integration, it is convenient to write (15) as 


== 2-4 — (hr? $c?) + (hr? — $c?)? 
0 


2 


© 


so that 
A+ 2 24(3A + 2 


i.e., t is of the form 


or 
(18) t= fq, + + agr’s + +}, 


where ¢ and r are henceforth to be regarded as complex variables, the coeffi- 
cients and the exponents, however, being real and r” for any exponent » 
referring to that branch of the function which is real and positive for positive 
real r. 

From (16) and the remarks preceding it, it is seen that ¢ is a regular 
analytic function of r in the domain 0 < | r | < r**, Moreover, from the 
uniform convergence of the integrated series, it follows by the approximation 
theorem for almost-periodic functions of a complex variable’ that ¢ is an 
almost-periodic function of r in that domain. Furthermore, if d is irrational, 
t has an almost-periodic singularity (corresponding to a logarithmic branch 
point) at r= 0, as is seen by applying to (17), where a, > 0, the reasoning 
used in case (I A al). Also, since A > 2, (16) shows that the exponents 
of the series for ¢ are all positive and that amongst them there is a smallest, 
viz., $(A +2). Consequently, f(r) is an analytic almost-periodic func- 
tion in the domain 0 < | r | < r** and has a normal almost-periodic singu- 
larity at the point r 0. Then, by a theorem of Bohr,® for a sufficiently small 
r* > 0 the function {= f(r) will represent the neighborhood of r=0 
on a simply connected domain D of the Riemann surface S;, which domain D 


° Cf. H. Bohr, “On the inverse function of an analytic almost periodic function,” 
Annals of Mathematics, 2nd Ser., vol. 32 (1931), p. 249. 
°Cf. H. Bohr, ibid., p. 251. 
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will contain a complete neighborhood Hp. of t= 0. Hence in Eps there exists 
an inverse function r = g(t) and this function will (for sufficiently small p*) 
be analytic and almost-periodic in 0 < | t | < p* and have a normal almost- 
periodic singularity at 0. 

Furthermore, if A = $(A + 2) > 0 denotes the smallest exponent of f(7) 
and M > 0 the smallest exponent of g(t), then 


M = 1/A = 2/(A + 2), 


and if we write r4f,(r) and g(t) = ¢Mg,(t), the modules of f,(7) and 
of g,(t) will be connected by the equation 


Mg, = 


(where it is seen from (16) that Mf, is generated by A and 1). Consequently, 
the reasoning of case (I Aal), applied ubove to (17), now shows that r as a 
function of ¢ possesses, in the case A irrational, a logarithmic branch point 
at == 0, so that in this case analytic continuation of the motion through the 
instant of collision is impossible. 

If, on the other hand, A is rational, the right-hand member of (16) 
becomes a Puiseux series which may be inverted just as in case (IAa). As 
there, we conclude that if A= p/q, p and q relatively prime, then 


1* If p is odd, q arbitrary, real analytic continuation is possible ; 


2* If p is even, = 2x, then 
(2,*) If «+ q is odd, real analytic continuation is again possible. 


(2.*) If «+ q is even, no real analytic continuation is possible. 


The smaller zero of F'(r) is simple and corresponds to a maximum R, 
of r= r(t), which, by (3 bis), is reached in finite time.’ 7 decreases steadily 
with ¢ from r = RF, to r = 0, the point of collision. From (2) it follows that 
¢ increases steadily with ¢ as 7 increases from 0 to R, and then decreases to 0. 
If0 << r* < R, and if t*, t, be the times corresponding to r*, R, respectively, 


t 
then the finite value of ¢— ¢* (c/r*)dt for all shows 
t* 
that approaches a finite limit as r—R,. If t-0, then r—>0 and 


t 
k const., 4 > 2, 


80 that ¢ approaches a finite limit in this case also. 
With respect to the case (II A y ai,), it is seen that Rz, the larger zero 
of F(r), is also simple and corresponds to a minimum of r= r(t), which is 


9 
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reached in a finite time from a point on the path for which r= r** > R,, 
As ¢ varies in any finite range, the same is true of ¢, so that @ approaches g 
finite limit as t 0, if ¢ is made to correspond to r= As r—>-+ a, 
1” —~ (2h)4, so that (2h)4t and by (2) + ko, and 
stants), so that @ approaches a finite limit as {> -+ oo also. The motion js 
similar to that in cases (II A « a) and (II A ab). This completes the dis. 
cussion of the motion under the assumptions (II): c>0; (A): U=?r, 
(y):2<A<+ (a) h>0; (i) F(m) <0. 

If (ii): F(7)) =0, then, as already pointed out, there is a circular 
motion, corresponding to an equilibrium solution and asymptotic paths one 
lying inside the equilibrium solution and corresponding to case (i,) above, 
while the other lies outside and corresponds to (is). The former is a collision 
path behaving at the moment of collision as the path in case (i,) above, 
r increases steadily with ¢ and approaches ast likewise increase 
steadily with ¢ and becomes infinite as r— rp (i.e., as t—> ©) as may readily 
be seen from the integration of (2) in view of the fact that r remains bounded. 
On the outer asymptotic path, ¢ becomes infinite and ¢ approaches a finite 
limit as 7, starting from some r** > r,, becomes infinite, just as in case (i,), 
but as r decreases toward becomes negatively infinite and decreases 
indefinitely with ¢. 

If, finally (ii): > 0, then F(r) has no zeros in 0<r< +4, 
so that there exists an unbounded collision path. The nature of the collision, 
corresponding to t= 0, is exactly the same as in subcase (i,) above. As 
increases from 0 to + r increases steadily and ~ const. while ¢ 
approaches a finite limit, this portion of the path being similar to that of 
case (iz). 

Let assumptions (II), (A), (y) be retained, but let (a) be replaced by 
(b):kh=0. Here F(0 +) =+ o, but F(r) < 0 for all sufficiently large r. 
F(r) crosses the r-axis an odd number of times, hence just once, as F’(1) has 
but one positive real zero. The motion corresponds to a bounded collision 
path precisely similar to that of subcase (i,) under (a). This concludes the 
treatment of case (II A). The only equilibrium solutions or asymptotic paths 
are those treated under subcase (aii). 


tetaining the assumption (II): c > 0, we now suppose that 


(B) U =—logr. 
Then 
F(r) =—logr— $e?r? +h; (r) =— + er; (r) = 1? — 
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F,=0 implies r—=c(>0), while F,, = according as 3c’. For 
<= 7 . . 

all h = 0, we have F'(0 +-)== F (+ o)==— o. F(r) has one maximum, viz., 


5. 
at r= c, and F'(c) = 0 according as h loge + 3. 
Consequently, (i) if h < loge+ 4, there is no solution; 


ii) if h=loge+ 4, there is a circular motion and 
oy 


nothing else ; 


(iii) if h > loge + 4, there is no collision, but a bounded 
periodic motion like that of case (II A a ci). 


Finally let assumption (II) be retained but (B) be replaced by 


(C) U = — 1", 0<A<+0 

Then 
F(r)=—- pr — 2 + h: F,(r)= +- —A(A— 1) r\-? — 
As in the previous case, for all h, F(0 +) =F (+ 0) =—oo. F, has just 

one zero, viz., at Ty = (c?/r)/O), Then 


(a): if 0 << A <1, it follows that F,, = 0 according as r =r— 


= 

respectively, and < since 1 < 3/(1— A). 

(8): if A=1, direct calculation shows that F vanishes at ry = c?/*, 
while — 3c?rt < 0. 

(y) 1fA > 1, F, vanishes at 7) = (c?/A)”**) and F,, < 0 for all r > 0. 

Hence the extremum is a maximum, and according as F'(r,) = 0, there exists 

(i) no solution; (ii) a circular motion and nothing else; (iii) a bounded 
periodic motion like that of case (II A a ci). 
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RAMANUJAN IDENTITIES INVOLVING THE PARTITION 
FUNCTION FOR THE MODULI 11+.* 


By JosEPH LEHNER. 


1. Introduction.‘ Ramanujan conjectured that 
(1.1) p(n) =0(modq) if 24n=1(modq), 
¢ 


where p(n), with n >0, is the number of unrestricted partitions of n, and > 
q is one of the primes 5, 7, 11. Direct proofs of each of these three con-f ; 
gruences were given by Ramanujan.’ 

He also indicated an indirect method of proof; namely, as immediate f 
consequences of certain identities, which he stated without proof, and off 


which the following is an example: I 

(1.2) — + 49x — 
n=0 II (1 Il (1 am ) 8 ca 

1 1 

The identities for the moduli 5 and 7 have since been proved by several la 
writers.* Ramanujan did not give an identity for g = 11. & 
2amanujan stated more than (1.1). His most general conjecture—tha 
to: 


(1.3) p(n) =0 (mod q*) if 24n=1 (mod q*), «= integer > 0 


—is known to be false for q = 7%, a > 2, but G. N. Watson * [15] has proveif ,¢. 
the correctness of (1.3) for g=5, «21, and has proposed and proved sf },), 
substitute congruence for g=7 and all a. In the same paper he promiself 
a discussion of the moduli 11%, but to date it has not appeared. A proof af 
(1.3) for g* = 11? is contained in one of Ramanujan’s notebooks but has no— >: 
been published.® For higher powers of 11, nothing is known concerning th (9 


* Received March 16, 1942. ; 

1The author wishes to express his thanks to Professor H. Rademacher, wh 
suggested to him the problem treated in this paper, and to Professor J. B. Rosser fof , 
helpful advice in connection with the computations of 9. 

* [14], pp. 210-213, 232-238. (Square brackets refer to the bibliography at 
end of this paper.) 

* For references see footnotes 3, 4, 4a in [11]. 

‘See also [16] and the reference in footnote 3 of that paper. 6 

5 [4], p. 100. 


492 


— 


RAMANUJAN IDENTITIES FOR THE MODULI 11%. 493 


correctness of (1.3), except for an affirmative test by D. H. Lehmer ® of an 
jsolated value. 

The procedure of Watson and Zuckerman was to establish identities of 
type (1.2) containing a factor q* in the right member. Recently Rademacher 
| [12] developed a systematic method for the investigation of these identities, 
§ and applied it to the moduli 5%, 7, 7°, 13°. Rademacher’s method displays 
very clearly the underlying connection between the Ramanujan identities and 
7 the theory of modular functions. But it is not so powerful as Watson’s 

q process, Which employs modular equations: it does not establish, except for 
| the first few values of a where direct computation. is possible, the integrality 
| of the coefficients of the modular functions (which themselves have integral 
| coefficient expansions) in the identity, and so the congruence (1.3) does not 
follow. 
In the present paper we use Rademacher’s method to establish identities 
' for the moduli 11%. For «=1, the coefficients in the identity are ‘com- 
| puted—6—while in the case a 2—9—they are determined sufficiently to 
| prove the congruence (1.3). For higher values of « we cannot prove the 
| integrality of the coefficients in the identities and so cannot prove (1.3). 
| Thus, our goal, the proof or disproof of the congruences (1.3) in the unknown 
b cases (i.e., > 2), is not reached. 

It is felt that the identities possess some interest in themselves, particu- 
larly since they highlight the arithmetic character inherent in the larger part 
| of the theory of modular functions. From this point of view the Ramanujan 
congruences appear to be an arithmetic expression of certain facts pertaining 
+ to special subgroups of the modular group and the functions belonging to them. 
2. Outline of the method. We shall make use of a simple modification 


of an operator studied by Hecke.’ Let F(r) be an entire modular function 
| belonging to the subgroup T',(q) of the full modular group. To(q) is defined 


by the condition c= 0 (mod q) in the substitution( 2 . where g is a prime 
>3. We define the linear operator 
A=0 

In general, we cannot prove that U,F is entire modular on I'y(q). How- 

sever, if F exhibits a simple behavior under the (non-modular) substitution 


=—1/gr, such as 


(2.2) F(—1/qr) =cF(r), or F(—1/qr) =cF4(r), c= const., 


*See [9], p. 90. 
*See [5] for a readable sketch of this theory and [6] for a complete account. 
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then U,F will belong to [%(q). Moreover, in certain cases F may be a 
modular function of T'y(qg?), and not of T(q), but still the same conclusion 
will hold. This is true of the function 
n(7) 


where (7) is the well-known Dedekind function 
. . 
(2. 4) n(r) = exp(mir/12) (1—a"), 2 2zir. 


@(7;q°) is an entire modular function on but q°) belongs 
to the larger group T(q). 

We now take g = 11 and show in 5 that every entire modular function 
of T(11) is a rational function of certain basic functions A,,(7), C,(r). 
The latter have power series expansions in x = exp 2mir with integral coeff- 
cients. We find in 6 that 
(2.5) = 11(11AC? — 11°C + 2AC — 320 — 2), 

A=w=Ai:(r), C=C,,(r). 

On the other hand, we have as a special case (q—=11, «=—1) of 

Theorem 7, 


(2.6) 121) (1—2") > 6)a', x 


i=0 
This, together with (2.5), gives 
(2.7) + 11 J] (1 — 7 (11 AC? — 11°C + 240 — 320-1 
1=0 1 

From (2.7) we deduce, by comparison of coefficients in the two members, 
the correctness of Ramanujan’s conjecture (1.1) for g=— 11. 

The identity (2.5) is equivalent to (2.7) and is far more convenient, 
since it involves only modular functions. We shall obtain all identities in 
the form (2.5), the transition to the partition function then being made by 
analogues of (2.6). 

I. The Identities. 

3. Preliminaries. By the statement “F(r) is an entire modular func 

tion on a subgroup T of the full modular group” we mean that F satisfies 


the following conditions: 
(A) F(r) is an analytic function of + which is regular for &r > 0. 


(B) F(r) has at most polar singularities, measured in the uniformizing 


variable, at the parabolic vertices of I. 
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(C) For every substitution 
ar b 
Vr= —— 


 ertd 
belonging to T, we have 
(3. 2) F(Vr) = F(r). 


The subgroup of greatest interest to us is Ty(q), defined by c= 0 (mod q) 
in (3.1), where g is a prime > 3. Because of the primality of q, the funda- 
mental region of Ty)(q) has only two parabolic vertices, which we may take 


to be = 0, 10. 


The uniformizing variable for is evidently z= exp 2miz, since 
T,(q) contains the substitution St—7-+ 1. For r’ 0, we use the sub- 
stitution r’ = — 1/qr and employ the transformation equation of (2.2). This 
will give us an expansion in 2, so that x is a uniformizing variable also for 


7 


Notation. + is a complex variable with 2s >0. The point rico is 
the infinite point of the r-sphere considered as the limit of iy, y>+ o, 
y>0. gis a prime >3. 8S, 7, V, W,- - -, are linear sub- 
stitutions, 7’ = (ar + b)/(er +d), with integral a, b, c, d and determinant 
ad—bce=1. Y(1) is the full modular group, consisting of all such sub- 
stitutions. is the subgroup of (1) defined by c==0(modn). O(2*) 
denotes a power series in x with integral coefficients beginning with a term 


of order => s. 


Formulae. 


20 +X 

(3. 3) ace! ( ]) /? — 
1 x 

(3. 4) [I] p(n)a", p(0) =1. 

(3.5) n(r +1) = exp (mt 12) -n(7), n(r + 24) = 7(7). 

(3. 6) n( l/r) = ( it)4n(r), 


where the square-root is positive when 7 is purely imaginary. 
For every modular substitution Vr= (ar+ b)/(er +d) eT(1), with 
¢>0. we have 


a-+d 
(3.7) (Vr) = exp — (s(a.0) “| [—i(er+ d) |4n(r), 
2c 


*The substitutions of the group determine the form of the uniformizing variable. 
Thus, for the uniformizing variable for r= 0 is a’ =exp(—2mi/qr). It is 
only because the functions considered here have an additional property of the type 


(2.2) that x = exp 2ri7r can be used in place of 2’, and it is convenient to do this. 
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where the sign of the square-root is determined as in (3.6), and where s(a,c) 
is a “ Dedekind sum.” The Dedekind sums are fully discussed in [13]. 


(3.7) is the transformation formula of (7) and was discovered by 
Dedekind and Hermite. (3.6) is a special case of this. (3.5) is obvious 
from the definition (2.4) of (7), and the first two formulae are classical. 


4. Modular functions of [,(q). Our first object will be to develop a 
set of modular invariants of [)(11). We shall, however, speak more generally 
of T.(q), g being always a prime > 3. We consider, first, modular forms 
of negative dimension, from which modular functions (zero dimension) will 
be obtained by division. Since we need entire modular functions, we must 
have one modular form which is zero-free in the interior of the fundamental 


region. Such a form is provided by 


THEOREM 1. The functions 


(4.1) = }" 
where 
(4. 11) r(q +1) =0 (mod 24), 


are regular and zero-free in the upper half r-plane. They have the trans- 


formation equations 
(4. 12) Wqr(Vr) (a/q)"(er + d)"bar(r); VreTo(q). 
(4. 13) tar(—1/qr) = (— "Yar(7). 


(a/q) is the Legendre-Jacobi symbol. It will be indicated in this way when- 


ever confusion with the rational number a/q is not to be feared. Otherwise, 


a 
we shall use the symbol (a|q) or ~——. 


Theorem 1 is proved in Part II. 


A second set of modular forms, of dimension —2r, is given by the 
Eisenstein series 
Dy» +00 
(4. 21) gr (7) (mr+n)", r= 2 
(27) m,n=-OC 


with the well-known property ”? 


(4. 22) gr(Vr) = 


®>’ means in this case that the term m = n=O is to be omitted. 

19 These functions can be defined for r= 1, also, even though the series on the 
right is no longer absolutely convergent. See [7], p. 207 
11 [3], pp. 20-22. 
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They have the Fourier series *? 


(4. 23) — Br/2r + 2(— 1)" 
ai(n) => d', B,—1/6, B,—1/30,: - 
d\n 


The forms (4.21) belong to the full group ['(1) whereas we are interested 
in functions having the modular property on no group containing T'o(q) as a 
proper subgroup. But consider 


(4. 31) Gar(t) =9'9r(qr) —gr(t), 
Clearly 

(4. 32) Gqr(Vr) = (er + d)*Gor(r), VreTo(q), 


Still another class of modular forms can be obtained from the multiple 
theta-series. We assume from now on g==3(mod4). The function 


+00 
m,n=-O00 
is a modular form of dimension —1 with the properties '* 


(4.42) 6(Vr) = (er+d)(a/q)@(r),  (a/q) = Legendre-Jacobi symbol, 
(4. 43) 6(— 1/gqr) = — iq’/*78(r). 


For higher dimensions we use either powers of 6g or theta-series in more 
variables.** 

From these modular forms we construct entire modular functions by 
dividing each form by a y¢,r(t) of appropriate dimension. In particular, we 
cebtain the functions of the following theorem." 


‘2 A convenient method for treating the transformation of @, is developed in 
Hecke [8]. We first reduce Q to a sum of squares: Q = 4((2m-+n)? + qn’) 
=%(M*+qN*), where M=N(2). Hence Q=wyu'/qg, where uw belongs to the 
principal ideal (V—q) in K(V—q). We then have (7) = 0(r; 0, (1), V—gq) 
in Hecke’s notation (p. 222). The relations (4.42), (4.43) of the text are equivalent 
to Sitze 7 and 5, respectively, of Hecke’s paper. If qg=1(mod4), we must use 
quadratic forms of determinant —4q, rather than — q. 

13 The transformation of theta-series in more than two variables can be reduced 
to that of double theta-series. See [7], pp. 210-11. 

14The functions which can be formed in this way from modular forms are not 
necessarily linearly independent. The number of linearly independent modular func- 
tions having at most a given valence can be determined from the Riemann-Roch 
Theorem applied to the Riemann surface which is the map of the fundamental region 
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THEOREM 2. The functions 


6," (7) 


Ger (7) 
Wa,2r(7) 


are entire modular functions belonging to the subgroup T(q) : 


(4.52) Ag(Vr) =Ag(r), Ba(Vr) =Ba(r), VreTo(q), g=3(4). 


(4. 51) Aq(r) = By(t) = 


Furthermore,® 
(4. 53) Aq(— 1/qr) = Ag(r), Bq(— 1/qr) = — By(r). 
The valences of Aqand By are 2sq and 4s8q, respectively, where sqg=1(q +1) /24. 
We need to prove only the statements concerning valence. The valence of 
a modular function is defined to be the multiplicity of the root of the equation 
F(r) 
where c is an arbitrary complex number and +r ranges over a point set 
consisting of the interior of the fundamental region and a complete set of 
inequivalent boundary points. The multiplicity must be measured in terms 
of the uniformizing variable. We may take c=. Since A,g(r) and B,(r) 
are entire modular functions, they can have poles only at the parabolic points 
=0, 00. 
We can find the expansions of Ag and Bq from (4.1), (4.11), (2.4), 
(3.3), (3.4), (4.31), (4. 23), and 


oo 
(4. 61) = Ag(n) 2x", v= exp 


where Ag(n) is the number of representations '* of n by the form Y. Thus, 


for r=10, 


of (q). Since (11) is of genus one, it follows from the Riemann-Roch Theorem 
that there are m — 1 4+- 1 = m independent functions on r,( 11) of valence at most m 
if m > 1. In particular, there are 2 independent 2-valent functions (since univalent 
functions do not exist on Riemann surfaces of positive genus 1) which may be taken as 
A,,(7) and A,,(7) +1. Now, we can construct another 2-valent function, namely, 


where lg, (llr) — 9, (r) (see footnote 10). Hence 
there must be a linear relation between A,, and A’,,. Indeed, we find by comparing 


power series that 6A’, = 5A,, —§8. Similar relations hold among the m-valent fune- 
tions when m > 2. 


15 Expressed otherwise, A,(7) is an automorphic function of the 


r(q) =T,(¢) +W-T, (9), W=—1/qr. Of course, I'(q) is not a subgroup of the 


group 


modular group. 

16 Counted in the usual 
only in the order or in the sign of the unsquared summands. For special Ag (I) 
explicit formulae are known: for Q= m2? + mn + 3n’, det.Q=—11, we _ have 


way, i.e., two representations are distinct if they differ 
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(4.62) Ag(r) = exp {— (wir/12)r(q + 1)} pees, 

(4.63) Bq(r) = exp {— (wir/12) -2r(q+1)} 

At the other parabolic point, we use 7’ = —41/qr and (4.53), and have 
(4.64) Ag(r’) =Ag(r) =o +---, =— =— 


Hence the valences are 2sg for Aq and 4s, for Bg, as stated. 
We conclude this section with a proof of the fact, noted in 2, that 
U,®(r;q*) is an entire modular function of T(q). We recall the definition 


q-1 
4.7) U P(r; q’ 
) q 7") ] q q 


Condition (A) in the definition of an entire modular function — 3— is cer- 
tainly fulfilled since Ug® is simply a sum of functions, each of which is 
regular in the upper half-plane. ; 

In the sum of (4. 71) we can replace the range of summation of d by the 
equivalent condition: Amodulo q. This follows from the fact that ® is periodic 
with period 1; replacing A by A+ q leaves each summand unaltered. Hence 


24 

Using (2.4), (3.3), (3.4), we find 

(4. 73) (l+a¢-+ 2277+ ) 


Introducing (4.73) in (4.72), we have 
9°) dn exp + 24A/qg)), Amodg, 
dn exp exp (2rin(24A/q) ) 
X 


= $m exp 2xi(mr/q), q|m. 


x 
(4.74) Ug®(r; 97) = 3 dan exp = dont", v =[u/q] +1. 
n=v n=v 
Ay (1): 2 > d,— = d,}, where d, = quad. residue, d, = quad. non-residue, (mod 11), 
d,\l tly 
and 1 5£0(11). See Dickson, Jntroduction to the Theory of Numbers, 1929, p. 81, 
ex. 12. For l llk, we use the substitution m = M + 6N, n = — 2M —N, and deduce 


easily Ag (11k) {,(k). Hence for l= 118k, Ag(l) = Ag (k,). 
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Thus for r= 10, Ug has an expansion of the required type. The effect of 
Ug is simply to replace an by dgn. 

For =0, we write ——1/r and apply the transformation 
equations (3.6) and (3.7). First, 


where 7r = (7 + A)/q and where we have made use of (3.5). Hence, 
, 2 7 


q-1 
+ (T24T'r) }. 


Now 
(4. 82) ToT = VT 
where 
— + 1)/g 
(4. 83) Vp ( 
(4. 84) =— 1 (mod q), 


We see from (4. 84) that » runs with A over a reduced residue system modulo q. 
Using (4.82) in (4.81) and applying the transformation equation of (7), 


we get 
q-1 
(4. 85) 
q-1 
= + Cun (Tost) 
pal 
with 


Cu = exp {ri(s(24A, gq) —2(A— z)/q)}. 
A computation similar to [12], (5.4), shows that 
(4. 86) (6n/q) =alv/q), 


Hence, the sum in the last member of (4. 85) is 


9 9 
—« Dd exp (— ) (u/q) p(n) exp | , n=d, 
q-1 
=e, exp (— wtr/12q) p(n) exp (2zi(nr/q) ) (u/q) 


n=0 


7 


x exp {2mi- (24n —1)p/q}—since (0/q) = 0— 
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00 
= exp (— > 


n=0 
where we have made use of the Gaussian sum 


(4. 87) G(h,q) => (1/q) exp (2ri(hl/q)  Imodg, 


l 


) p(n) exp = 1 


We replace + by qr and obtain finally 
g(r’; = 1 (7; g?) + exp (— zir/12) (1+ O(2)) 
= = —1/qr, 


by (4. 73) and (2.4). Hence, Ug® has a pole of order u in x = exp 2air for 
7 =0, and condition (B) is satisfied. 

That U,® fulfills (C) with T'—T,(q) is seen as follows. From Theorem 
8, which is proved in Part II, we take the fact (a = 2) that 


(4. 91) ©(Vr; = q7), V T(q’). 
We require also the following lemma, proved in Part II. 


LemMa 1. To any VeT,(q) and 0SA < q, there exist 
and WneTo(q"), such that 


(4. 92) =WpTp, = (+ +A) 
Furthermore, » runs with d over a complete residue system modulo q. 
Therefore, 
(4, 93) U,b(Vr; + 9°); A mod q 


= 9? 9?) = (739), 


as required. This proves (C) and 


THEOREM 4. The function Ug®(r; q*) is an entire modular function of 
the group T,(q). 


5. A rational basis for [,)(q). The next step is to find a basis for the 
entire modular functions of T(q). For this purpose the following lemma 
will be useful.*? 


Lemma 2. Let fi(r), fo(r) be modular functions on a group T having 


7 For proof see [1], Th. 13, pp. 96-97. 
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the property that an arbitrary pair of values of the functions is assumed at 


but one point of each fundamental region of T, i.e., the equations 
fi(t) =f.(7’), fo(r) = f2(7’) 


shall imply that r and 7’ are congruent under T. Then every modular function 
of T is a rational function of f, and fo. 


We are now ready to prove 


THEOREM 5. If the function Ag(r) of (4.51) ts 2-valent on T)(q), 
then the functions Aq(t), Bg(r) constitute a rational basis for T)(q). 


Proof. By Lemma 2 we have only to show that 
(5. 1) Agq(r) == Aq(7’), Ba(r) == Ba(7’) 


is impossible if r+47’ and 7,7’ lie in the same fundamental region of T)(q). 
Since Ag(7r) is 2-valent it follows from (5.1), (4.53) that 7’ =— 1/qr. 
But then 

Ba(r) = By(— 1/97), 


which, in general, is inconsistent with (4.53), since By(r) is not identically 
zero.'* 

As a special case we take g = 11 =3 (mod 4), n = 2, s,, —1, satisfying 
the conditions of Theorem 2. Then A,,(7) is 2-valent on T)(11) and we have 


CoroLuary.’® The functions Ai:(t), form rational basis for 
the entire modular functions of T)(11). 


Note that 
(5. 2) val A;, = 2, val B,, = 4. 


6. An identity for modulus eleven. Since U,,6(7;11°) is an entire 
modular function of the group [)(11) (Theorem 4), it is a rational function 
of Aii(t), Bii(r): 

(6.1) = Bi). 


It will turn out, however, that P is a polynomial. That is, we assume P as a 
polynomial with undetermined coefficients, and then show how to determine 
these coefficients so that the difference U,,6 — P, which is an entire modular 


18 This argument shows, incidentally, that B,(r) must vanish at the fixed point of 
the substitution 7’ = — 1/gr, and indeed at the fixed points of all elliptic substitutions 
of the group I'(q) (see footnote 15). 
1° See [2], p. 406, (12). Fricke denotes A(r) by r(w) and B(r) by o(). 
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function, has no poles. Then it must be a constant, which can be evaluated at 
Say. 
The expansions of U,,® at the parabolic points were worked out in the 


proof of Theorem 4. We repeat them here. 

(6. 2) (7,,®(7r; 117) = O(z2) 

(6. 3) 11°U,,8(7'; 117) = a5 + O(1). = —1/1I1r, 
where we have used the expansion 


(1217) 
calculated from (2.4), (3.3), (3.4). 

We are now in a position to compare the principal parts of the two 
members of (6.1). At r=1t0, U,,® has a zero of the first order, whereas 
A,, =A and B,, = B both have poles, by (4. 62), (4.63). It will therefore 
be convenient to construct a new function, C,;,(7) = C;;, which has a zero of 
the first order at r—1i0. We wish C to have the lowest possible valence, 
i, e., 2, since the fundamental region of (11), considered as a closed Riemann 
surface, is of genus one. In other words, C must have a double pole at 7’ = 0. 
If we add the requirement that the expansion of C at rt shall begin 
with x, the function is then determined, by the Riemann-Roch Theorem.” 


From the expansions—see the table on p. 514— 


(6. 41) A(r) =2'+6+ liz + 467° 
(6. 42) B(r) = 2? + — 12 — 1162 — 59722 
we find 
(6. 43) C(r) = — 10A (+) — 22 — B(r)). 
Its expansion at tr = 10 is 

oO 
(6. 44) C(r) Benz". 

n=1 


From the last relation and (4.53), we see that 

2-11°C = A*(r’) —10A(7’) — 22 — B(7’), =—1/11r 
(6. 45) 
A?(r) —10A(r) —22 + B(r). 


*°The number of linearly independent functions having a pole of second order at 
t=0 is 2—1-+1=2. The most general such function contains two arbitrary con- 
stants, one of which is used to insure the zero at 7 =—1i%, the other to guarantee the 
coefficient 1 for the linear term. 


a 
e 
| 


504 JOSEPH LEHNER. 


This gives 
co 
(6. 46) 11°C =a? 4+ 4+: D cnt". 


n=-2 


LemMaA 3. The coefficients Cn, ¢'n of (6.44), (6.46) are rational integers. 


Proof. Since A,B both belong to T(11), they are connected by an 
algebraic relation, which is of degree 2 in B and of degree 4 in A, in accord- 
ance with the valences (5.2). By writing such a relation with undetermined 
coefficients and then using the expansions (6.41), (6.42) we find easily 


(6. 51) B? = A*— + 56A?— 444A. 
From ** (6, 43) 


(2-11°0 + B)?= (A2—10A — 22)? A*— 4+ 564? — 444 + 484(41 

== (mod 11’), 
2-11°C + B=B (mod 2-11?) 
or 2-117°C=0 (mod 2-117). 


Thus the c, are integers. 
Again, let C(—1/11r) = C’, etc. Then, from (6. 45), 


2-11°C’ —B = A*?— 10A — 22, 
(2- 11°C’ — B)? = (A? — 10A — 22)* = (mod 4), 
2-11°C’ — B = R (mod 2), 11°C’ = B (mod 1). 


This proves the lemma.’ 
For future reference we include another lemma. 


Lemma 4. AC —1=0 (mod11). 
Proof. First, 
224° — 55A* + 444° =0 (mod 11) 
or At(A? + 2A + 1) = A?(A*— 20A* + 56A? — 44A) (mod 11) 
== A*B* (mod 11). 


21 In what follows we employ the following simple theorem: Let F = + 
G=a-l+ Xba, »>—l; a, and b, are integers. Then if F? = G* (modk), k odd, 
we have F=G@(modk); if F? = (mod2\), then F=G@ (mod 2-1). In 
general, we would have F=+G. The correct sign can be determined from the first 
terms of the power series for F and G. 

22 The functions and = C(— 1/117) do not satisfy so simple a transformation 
equation as A and B do. However, this equation can be found easily by eliminating 
A and B from the three relations (6.43), (6.45), and (6.51). The result is 
— 12.1100’ —11(C + C’) —1=0. 


Wi 
(6 
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A*(A + 1) = AB (mod 11) 
22-11242(A + 1) = 2?- 11°AB (mod 2?- 11°) 
A? (B? + 27-117(A +:1)) = 2-114 + 2?-11°AB + A?B? (mod 2?- 
A?(A? — 10A — 22)? == (2-112 + AB)? (mod 2?- 11°) 


if we use the computation in the lines following (6.51). Hence, taking 


of footnote 72, we obtain 


A(A* — 10A — 22) == 2-11* + AB (mod 2 - 11°) 
A(A*? — 10A — 22 — B) = 2-11? (mod 2- 11°), 


and by (6. 43) 
2-11? (mod 2-11*), 


which completes the proof. 

Following the procedure outlined at the beginning of this section, we set 
U,,6 equal to a polynomial P in A and C with undetermined coefficients. 
The poles of even order in the expansion (6.3) of U,,® at 7’ =0 can be 
| eliminated by powers of C, those of odd order by A times a power of C. Since 
both U,,® and P are bounded at +r i, it is clear that the term in 2? 
(at 77 0) must drop out automatically. Otherwise, the difference U,,# — P 
| would be an entire modular function having a single pole of the first order on 
T,(11), a subgroup of genus one. 

We therefore write 


11°U,,®(— 1/1l1r; 117) = O(1) 
= 2A C0"? + + + + doo, = C(—1/11r) 


= O(2) = + + a, AC + + doo. 


This gives us the following five equations, which are solved seriatim beginning 


with the uppermost. 


11-* 
11-* + aoe 11-4 
21-11-* + + a,,°11° 
— 20-11-74 + =—0 
+ doo = 0 
The values of the a’s are: ==— 115, ay, = 2-11, ao: — 32-11, 
= — 2-11, 


This constitutes the proof of 


(6.52) 


10 
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THEOREM 6. For g=—11, there is a “ Ramanujan identity” 


2). 


(6.6) = 11(11A0? — 11°C + 2A0 — 320 


Ramanujan’s congruence (1.3) with q* = 11 follows from (6.6) and (2.6), 
which together yteld (2.7). 


7. An induction. We now wish to establish identities of the same sort 

as (6.6) but involving higher powers of 11. But it is not yet clear what : 
functions, analogous to U;,;®(r; 121), should appear in the left member of 
such an identity. Our guide in this respect will be formula (2.6), which 
connects the partition function with U,,. The first step, therefore, will be 


un appropriate generalization of this formula. 


THEOREM 7. Let 


(7.11) L(r;1) 


L(r; == UgL(r; B=0. 


Then fora =1, 


(7.2) L(r3q*) =x (1—2™) t+ p)z', = exp 2mir, 
1 1=0 


where 

(7. 21) t= t(q*) for q < 24, 
29 1, even 
: q, « odd 

and 


(7. 23) 24p = 1 (mod q*), p< q. 


(2.6), is the case «1 of (7.2). 


We replace (7.2) by two statements, 


240 
(7. 31) = (+>), 


Proof. 


A mod q*, 


00 + 24A 
(7. 32) TT (1 — 2) + ( ), 


and proceed to prove the first. Suppose a1; then BO, and (7.3I— 


reduces to 


24d 
(7. 4) Ug®?(r; (=+**), A mod q. 


23 Actually, the value of ¢ is 1 if a is even and [q/24] + 1 if a is odd. 
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Now (see the lines following (4. 71) ) 


U,®(r; 9°) ‘2 ® : A mod q, 


24r 24r 
q 
where we have used (3.5). This is (7.4). 
Suppose now that (7.31) has been proved for «28. From (7.12) 


we have 


{ (7) 


= q PU, ) } 


‘ 9 : 
= Sy (° Sy 24a) ) 
q q 
24 


which establishes the inductive step from 28 to 28+ 1. The step from 
28+ 1 to 28 + 2 is accomplished in the same way. Thus (7.31) is proved. 
To demonstrate (7.32) we employ the formula 


n (7) 1=0 


which follows immediately from (2.4) and (3.4). Introducing this in the 
right member of (7.32) we see that it becomes 


2 2. 
(er) exp ( — )3 p(n) exp 4) , Amod q* 
r 


= (er) exp (- = p(n) exp (2-122) exp (24n —1)) 


vmod 


= exp (— + (i —2") > m) exp ( 


24m = 1 (mod 
from which (7. 32) follows readily with 


24q% 


(7.6) 


Clearly ¢ is an integer. For, eq* = q?8 = (q®)?==1 (mod 24) for «= 2£, 
and eq? q%*! = (qSt!)?==1 (mod 24) for 1; moreover, 


. 6), 
sort 
what 
rhich 
ll be 

mr 
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q* divides 24p—1 by (7.23). Since 0 < p < q* and since we assumed in 
(7.21) that q < 24, we find that 


0<t1<1+9/24< 2; 


therefore, as an integer, £1. This completes the proof of the theorem. 

Theorem 7 shows us that we must attempt to construct identities in which 
one member is L(7;11*) while the other is a polynomial in A and (, 
Obviously we shall proceed by induction, in view of the recursive definition 
(7.12), (7.13) of the L’s. The first step of the induction (#1, B = 0) 
has already been proved; it is Theorem 6. 

The nature of the inductive step from a to a+ 1 will depend upon the f 
parity of a, as the equations (7.12), (7.13) clearly show.** However, we can | 
bring these two cases together. Suppose that we have an identity 


(7.7) L(r;11*) = P(A(r), C(r)), 


where P is a polynomial in A and C with rational coefficients, in which 4 
appears at most linearly and always multiplied by a power of C. In order to 
obtain L(r; 11%), either we apply U,, directly to P or we first multiply by 
@(7;11°) and then operate with U;,. Since U,; is linear we may apply it to § 
each term of P or ®P separately. Thus the most general function which we § 
need to consider is 
(7. 8) = (7; gq?) (7) C*(r), 

j—Oorl, 
The induction will be completely established if we prove that G(r;1,),k) 
can be expressed as a polynomial of type P(A(r), C(r)) m (7.7). 

The procedure will parallel that used in the case ¢ = 1. The whole prooi 
consists in showing that G@ is an entire modular function of T(11). The 
expansion of Gat tr =i will not contain any terms with negative exponents. 
By examining the proof of Theorem 6 we see readily that any such function 
can be expressed in the form of a polynomial of the desired type. 


8. Identities for powers of eleven. 


THEOREM 8. Let G(r3;i,j,/) be defined by (7.8). Then G(r; 1, j,k) 
is an entire modular function of T,(11). Moreover, we have 


G(r; i, j,k) = P(A,C) 


where P is of the same type as in (7.7). 


24 See [12], pp. 624-625. 
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Proof. We must show that conditions (A)-(C) of 3 are fulfilled. This 
is obvious for (A) since all the functions involved in G are regular in the 
upper t-half-plane and U,, is simply a sum of such functions. 

We now consider the expansions of @ at the parabolic vertices. At 
7=10, we have 
(8.1) 11°) = O(a), Ai(r) = O(a), = O(a), 
by (4. 73), (6.41), (6.44). Then 


Di AIC* (7) = O(a”) = daz", w= 0. 
n=w 
Hence, as we saw in (4. 74), 
ao 
(8. 2) G(r; 1, j,k) == == diz", = 1. 
naw’ 
For 7 = 0 we write = Tr = — 1/r, and have 


10 
G(r’ 3%, j,k) = 117 7) 


A\=0 


10 
A=0 


10 
N=1 
10 
— 1/117) + 114 oaur), 
if we recall (4.82). The functions A and C are invariant under Vy, which 
belongs to T)(11), but ®(T sur) picks up the factor 
= (w/11) ; where (/11) is the Legendre-Jacobi symbol, 


as is seen if we put q = 11 in (4. 81)-(4. 86). This gives 


10 
(8.31) 11G(r’) = 1/11) + 11-776 (p/11) 
We write 


(8. 32) Pi (r; 11°) = > ; Ai(r) Anz", C*(r) Cnr", 
n=54 n=k 
where the coefficients are integers. Then 
11G(r’) = 1/111) 
fr + 24 
(p/11)dmancp exp 2rt(m + n+ p) ( 


Lem, 


The sum in the right member is 


exp (2xi(m + n+ p)(r/11)) (u/11) exp {2ri(m + n + p) (24n/11)} 


9 
m,n 


hich 
0; 
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0) 
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where we have evaluated a Gaussian sum, (4.87). Using (8.32) and replacing 


t by 117, we now have 7’ = — 1/117, 


(8.2) and (8.4) show that G@ has only polar singularities in 
x = exp 2xtr at the parabolic vertices. 
It remains to prove that 


(8. 5) G(Vr) = G(r), VeT)(11). 


Here we require (4.92) of Lemma 1. 


10 
G(Vr) V eT(11) 
A 


=0 


10 
= pid ( ) 9 Wa € T, (11°) 


10 
= 11-1 (Tyr) = G(r), 

since all three functions are invariant under T)(11*). This completes thef 
proof of the first statement of Theorem 8. 
That G can be expressed as a polynomial in A and C now follows in the 
same way as in the proof of Theorem 6. We set @ equal to a polynomial 


G = 3ajxAIC*, j= 0 or i. k = 4 


with undetermined aj. Then, equating coefficients of like powers of x in thi 
expansions at the parabolic points, we obtain a set of linear equations analogow 
to (6.52). The first equation will involve only one coefficient, and each sue 
ceeding equation will add one new coefficient. Hence, the coefficients arg 
determined uniquely. They are obviously rational numbers, since the expar 
sions of all functions involved have rational coefficients. This establishegg 
Theorem 8. 

Making use of the remark in the lines following (7.8), we see that wa 


have proved 


THEOREM 9. For every integer «= 1, there is a “ Ramanujan identity’® 
of the form 
(8. 6) L(r;11*) = P(A,C), 
where L is defined in (7.12), (7.18), A= An im (4.51), and C in (6.43) 


P is a polynomial in A and C with rational coefficients. 


By using Theorem 7 we can put (8.6) in the form 
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> p(q*l + P(A, 0) (1— 
=0 1 


with the p, « of (7.23), (7. 22). Unfortunately, we have not been able to prove 
the integrality of the coefficients of P except when ¢ = 1,2. Thus, for higher 
values of a, (8.7) gives no information concerning Ramanujan’s congruence. 
The discussion for « = 2 is given in the next section. 

The following remark may be of interest. If P(r) is an entire modular 
function of T)(q) and we subject U,F to the transformation +r’ =— 1/qr, 
we find, following the indications already given, 


(8.81) qU,F (— 1/qr) — (qr) = F(— 1/q’*r) — F(z). 

If, in addition, (+) has expansions in the same uniformizing variable (say 
= exp 2zir) at both parabolic points, then U,F will be entire modular on 
T,(q). Now suppose that Ty(q) is a subgroup of genus zero. This is true for 
5, 7,13. Then there exists a univalent function 


( (qr 
“4 Parle) (7) 


where r(qg — 1) ==0 (24) and r is minimal. Moreover, 


(8. 83) 1/qr) = (7). 


F(r) and U,F(r) are polynomials in ®g,(7), and this is true also of F(7’) 
and U,F(7’), by (8.83). Substituting these expressions in (8.81) we obtain 
an algebraic equation connecting g(r) and #,,-(qr). If this is irreducible, 
it is the modular equation of (7). For the primes mentioned, these equa- 
tions have been discussed by Rademacher [12]. 

Another approach to the modular equation is given in a forthcoming 
paper of the author [10]. 


9. The Ramanujan congruence for « = 2. For « = 2 we can determine 
the coefficients of P in (8.6) to an extent sufficient to establish the Ramanujan 
congruence. (6.6) can be written 

L(r311) = 112(AC?-— 110?) + 11(2(AC — 1) — 32C). 
Also, by (7. 13), 


The application of U,, to a power series simply replaces each coefficient an by 
(iin. It is clear, then, that 


L(r; 11?) =2-11U,,(AC —1) + 110,,C (mod 11’). 
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(9.1) 


(9. 2) 


(9. 21) 


(9. 31) 


(9. 32) 


(9. 4) 


(9. 41) 


have from (6. 45) 


JOSEPH LEHNER. 


L(r; =11U,,C (mod 11’). 


L(r; 117) =0 ( 


it is sufficient to show that 


mod 11?). 


= 0 (mod 11). 


0 (mo 


occurrence, and determine g last. 
The expansion of U,,C at r’ =0 is, by (8.81) and (6. 46), 


11°U,,C =a + 22+ O(1), 
and this we set equal to 


g + (7) = g (7) C*(7’), 


dil). 


or 1, 


Putting 11°C = FE (r) and bjx = cjxl1™**, we get 


Tt is therefore sufficient to determine the bj, modulo 11. 

The coefficients bj will be determined by considering the expansions at 
the parabolic point r’ 0, except for boo, which is to be determined at 
t=10. We therefore replace boo by g, a constant which may vary with each 


k>j, 


11°U,,C (7’) = > cjxAi(r) E*(r). 


Cj, = 0 (mod 11%), 


j=0,1. 


has an expansion in z with integral coefficients. 
Clearly the term of highest degree in FZ in (9.4) will be Z". We shall 
therefore need the power series of #7 and A to 22 terms. To compute F we 


But by Lemma 4, each coefficient of AC —1 is a multiple of 11; hence, 
By (7.2) Ramanujan’s congruence for g* = 11? is equivalent to 


Therefore, in order to prove Ramanujan’s congruence for the modulus 11°, 


This was done by a direct computation of the coefficients of the polynomial 
= Q(A(r), C(r)) = bind 


Since U,,C, A, and C have expansions (in 2 ~exp2zir) with integral 
3 5 


coefficients, the statement (9.1) reduces to 


7 1/1l1r 


We now see that (9.21) will be completely proved if we can show that 


(The constant, g, will be determined separately.) This is true since E(r) 


ial 
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(9. 5) A? —10A — 22 + B. 


Since we must compute cj, modulo 117, we shall need A and B modulo 2: 11?. 
The calculation may be carried out as follows (see the table on p. 514). 
The coefficients of A are computed from its definition (4. 51) 


Am Ay = (112). 
First, 7*(7)7*(11r) is obtained from (3.4); then this is squared. Next, 
the coefficients of @ are calculated from (4.61) and footnote **, and the 
resulting series squared. From these two series we obtain A and A? by 
multiplication and squaring. 
The definition (4.51) of B is 


Gi1,2(7) 


B= B,, 
Wi1,4(T) 


Yas 18 simply the square of ¥ 71, Gii,2 is defined by (4.31) and (4. 23). 
By multiplying the last two series we get the coefficients of B. 

The computation for each series (except that for 6? which has very 
small coefficients) was checked by recomputing the coefficients modulo 2- 11°. 

Now, using (9.5), we obtain the coefficients of #, reduced modulo 11?. 
We raise to-the eleventh power by the formula E- E?- 
working modulo 11*. In order to have a check on the coefficients of E™, we 
observe that, by an elementary theorem on binomial coefficients, if 


fo 
E(r) > eve’, 
y=0 


then 
le @) 
EX (7) (mod 11). 
y=0 
Therefore, 
(9.6) (7) = E(11r) = (— 1/1217) 


= + 29711 + O(1) (mod 11), 


where we have taken the last expansion from (6.46) with + replaced by 11+. 
(9.6) shows that all the computed coefficients e’, of H*t must be multiples 
of 11 except e’..==1(mod11) and e’_.;==2 (modi11). As a further check 
we computed £** (modulo 11) in two ways: 
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Residues of the Coefficients of Certain Power Series 
Power Function 
modulo 2-11? mod 11° 
0 1 1 1] 1 1 ] 1 1 
1 2 4 1 12 4 —2 2 2 22 
2 5 4 6 70 14 —18 —12 —12 —833 
3 10 8 17 54 40 —56 —116 5 —55 
4 20 20 46 105 105 96 —113 8 —bd5 
5 36 16 116 G2 10 —10 —4120 1 0 
6 65 32 10 —14 90 —20 —114 7 0 
7 110 16 19 116 30 38 110 —11] 22 
8 —b5bB7 36 66 —d6 — 82 40 10 10 22 
9 58 28 25 -—4 98 — 62 —12 —12 dd 
10 — 3 40 —90 50 — 56 —90 —18 —18 11 
11 28 4 —39 72 94 —2 —22 —22 2 
12 — 41 64 6 —35 —107 26 — 95 26 —ll 
13 86 40 —- Fi 0 — 82 —40 110 —11 —44 
14 23 64 2 — 60 90 100 — 80 4] 11 
15 —118 56 75 68 — 86 — 38 44 44 — 22 
16 86 68 90 — 83 — 21 76 —15 —15 0 
17 —112 40 30 —38 —114 94 — 96 25 — 22 
18 108 100 76 — 32 —6 —7 — 66 55 0) 
19 — 48 48 —%6 20 —120 7 54 54 33 
20 32 104 28 —11% 23 —4 87 — 34 5d 


21 — 94 80 — 36 — 4 104 96 36 —46 —ll1 


Setting cjx = 11fjx, we can now rewrite (9.4) as 
(9. 7) 11°U,,C =11"- + E*(7) + 9 (mod 11). 
We need to compute the fj, only modulo 11, prove 
(9. 71) fj: =0 (mod 11), 


and determine the constant g. 


This is done exactly as in the proof of Theorem 6. We obtain a set of & 


linear equations, analogous to (6.52), which can be solved stepwise for the fiz. 
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The coefficients of the various powers of A and EF are found quite readily 


since the computations are carried out modulo 11. The result is: 


(9.8)  112U,,C (r’) = 11 - BM — + + 6E® + 4E°A + 
+ B'A + 9E' + YE°A + + + 28°) + (mod 11). 


The computation was checked as follows: 


start with: 2HA + 8A+6 

multiply by 6H 

add 44 + 8: 8FA+6F+4A+8 
multiply by + + 6H? + + 82, 


etc. This procedure mixes the numbers quite thoroughly, and it is very 
unlikely that an error in the check would duplicate one in the original 
computation. 

Recalling that L(+) =11°C(7’), we see that (9.8) implies (9.71). 
Now we have only to compute g. Let r—>7’ in (9.8). Since A(r) = A(7’), 
we get 

=g + 1190(4,C) (mod 11) 
= g (mod 11), 


since Y(A,C) is a polynomial with integer coefficients. Since 
U,,0(r) 
we have 
g=0(11), 

which proves (9.1). By the remark in the lines preceding (9.1), this 
establishes 

THEOREM 10. p(1211-+ 16) =0 (mod 121), == 0,1, 2,° - -). 

Il. Proofs of certain theorems. 


10. Proof of Theorem 1. That the functions 


War (7) = 
are regular and zero-free in the interior of the upper half-r-plane is clear 
from the definition (2.4) of »(r). We consider their behavior under 
V=(ar+b)/(cer +d) eQTo(q). First, if c=0, Vr—7r+ 5), and using 
(3.5) we obtain 


Wqr( Vr) = exp { (wt/12) br(qg = War(r), 


because of condition (4.11). Thus (4.12) is proved for this case. 
If c > 0, we apply (3.7) and find 


22 
1 
22 
33 
| 
4 
11 
44 
1] 
22 
0) 
22 
0 
55 
ll 
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War(Vr) = N(a, d) (— + d) 
with 
(10. 11) (—1)"N (a, b, c,d) = exp (— mrA) 
(10. 12) A = (s(a,c) — (a+ d)/12c) 
+ (s(a,¢/q) — (a+ d)/(12¢/q)) +4. 


Therefore, in order to establish (4.12), we must prove 
(10. 13) rA =r-4(1— (a/q)) (mod 2). 


a/q is the Legendre-Jacobi symbol. 


We shall need the following properties of s(h, k): 
For (h,k) =1, we have 


(10. 21) 12hks(h, k) =h? + 1 (mod 6-k) 
(10. 22) 12ks(h, k) =0 (mod 3) if 3fk. 


Here for 3{k and = 3 for 3|k. 
For odd k& we have 


(10. 23) 12ks(h,k) =k + 1—2(h/k) (mod 8). 

Let & = 241, X= 0, and let / and h be odd. Then 
(10.24) 12hks(h,k) =h? + + 3k +1-+ 2k(k/h) (mod 
These properties appear as Theorems 17, 18, 19 in [13]. 


From (10.21) with k —c, c/q, we get 


12ac{s(a,c) — (a+ d)/12c} =a? + 1—a(a+ d) =— be (mod 6c) 
12ac{s(a,c/q) — (a + d)/(12c/q)} 
+1) —qa(a+ d) =— (mod 6c). 


Hence, by addition, 


12acrA =—r(q + 1)bc=0 (mod 6c). 
Also, by (10. 22) 
12crA =0 (mod 3) if 3fc. 
Therefore, we have 
(10. 31) 12crA =0 (mod 3c). 


Now suppose that c is odd. Then from (10.23) we obtain 
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12c{s(a, c) — (a+ d)/12c} =c + 1-— 2(a/c) — (a+ d) (mod 8), 
12c{s(a,c) — (a+ d)/(12c/q)} 
=c+q—2q(alc/q) —q(a+ (mod 8). 


If we note that the Legendre-Jacobi symbols satisfy 
(a/(¢/q)) = (a/(¢/q)) (a/q*) = (a/cq) = (4/c) (4/9), 
we obtain 
12crA = 2cr + r(q +1) —2r(a/c) (q(a/q) +1) 
—r(g+1)(a+ d) + 6re (mod 8) 


== — 2r(a/c)(q(a/q) +1) (mod 8) 
== — 2r(a/c) (1 — (a/q)) (mod 8) 


| 
| 


by virtue of (4. 11) 
== — 2rc(a/c) (1 — (a/q)) (mod 8) 


since c is odd and 4|2(1— (a/q)) 


(10.32)  12crA =— 2re(1— (a/q)) =6re(1— (a/q)) (mod 8). 
By (10. 31) 
(10. 33) 12crA = 6re(1— (a/q)) (mod 3c). 


From (10.32) and (10.33) we infer 


12crA = 6re(1— (a/q)) (mod 24c), 
which is equivalent to (10.13). 
Now let c= A > 0, yodd. Then, (10. 24), 
l2ac{s(a,c) — (a + d)/12c} 
=a? +c? + 8c +1-+ 2c(c/a) —a(a+ d) (mod 2***) 
12ac{s(a, c/q) — (a + d) | (12¢/q)} 
== ga? + c?/¢ + 3c 4- + 2c(¢/q|a) — qa(a+ da) 
(mod 2**), 


from which it follows that 


12acrA = re*?/q(q +1) + 6er—rbe(¢q +1) 
+ 2er(c/a) ((q/a) +1) + 6aer (mod 
== -+ 2er((q/a) +1) + 6er(a +1) (mod 2) 
(14 — (@+)} 2), 
(10. 41) 12acrA = (q/a) — a} (mod 2*°), 
This shows that 


(10. 42) 12crA = 0(2**) if ris even, 
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since a must be odd if ¢ is even, because of (a,c) = 1. 
If r is odd, we must have g ==— 1(4) and therefore 


(q/a) = (a/q) (— 1) =a(a/q) (mod 4). 
Hence, from (10. 41), 


12acrA = 2*rya((a/q) —1) (mod 
so that 
0 (mod if risevenorif (a/q) 


10. 43 12crA = 
( | 2\*2 (mod if (a/q) =—1 and ris odd. 


The congruence (10.43) proves (10.13) in case c is even. Thus the 
transformation equation (4.12) is established in all cases. 
To prove (4.13) we use the formula (3.6). 


Wa.r(—1/r) = g/t) 1/7) = (— 1)" (7/9) 
= (— 19-7) 


(4.13) now follows upon replacing + by gr. This completes the proof of 
Theorem 1. 


THEOREM 3. Let q bea prime > 3, « an even integer, and 


Then 


= ®(r) 
for every V= (‘ 3 belonging to To(q*), i. e., with c==0 (mod q*). 
Proof. We shall make use of the transformation equations (3.5) and 
(3.7). Ifc=0, Vr=7r-+), and 
(10. 5) (Vr) —exp (i/12) = ®(r), 


since 24| (q*—-1) by virtue of the conditions placed on q and a. 
For c > 0, we apply (3.7) and find 


(10. 61) (Vr) = M(a, b, c,d) ®(r) 
with 

(10. 62) M = exp (mH) 

and 


(10.63) H =(s(a,c) —(a + d)/12c) —(s(a, ¢/q*) —(a + d)/(12¢/q*)). 


E 


the 


and 


RAMANUJAN IDENTITIES FOR THE MODULI 11%. 


We have therefore to show that 
(10. 64) H =0 (mod 2). 


The proof of (10.64) is word for word the same as that of Lemma 5 in 
Rademacher [12] if we merely replace 2 by a. 
Proof of Lemma 1. With 
a b a, b 


we must have 


qc qd Cy Ci qd; 
or 
10.7) a, 
= + qd, = qd. 


Equations (10.7) determine c, and a, uniquely with c=0 (mod q*) and 
(a,,¢:) =1. From the equations on the right we find 
(10. 8) pa = b + Ad (mod q). 
Since q|c, gfa, and (10.8) together with the condition 0 < » < q determines 
p uniquely, Furthermore, =. implies (A, —A.)d=0 (mod q), which, 
since gf d, is equivalent to A; =A, (mod q), i.e., Ay = Az. 

Now we find d; from (10. 7): 


d, = d — pe,/q = d — pe. 
If c > 0, we have 
(10.9) (a+Ac)(d—pe)—1 _ c(b + Ad — pa) — prc? 


since ad —be = 1. b, is an integer, as we verify from (10.8) and (10.9). 
If c= 0, we necessarily have a=d=—1, a, =1, ¢, =0, p=b+A(Qq). 
From (10.7), = (b +A—p)/q, an integer, and a,d;— bic; =1. Thus 


the lemma is proved. 


CORNELL UNIVERSITY, 
ItHaca, N. Y. 
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